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ABSTRACT 

This  thesis  investigates  stable,  asymptotically  stable 
and  unstable  dynamical  systems  from  a  topological  point  of 
view  with  direct  application  and  interpretation  to  systems 
of  differential  equations.   Knowledge  of  topology  is  not 
a  prerequisite.   Specifically,  such  concepts  as  Poisson 
and  Liapunov  stability  as  well  as  parallelizable  and 
dispersive  systems  are  investigated. 
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INTRODUCTION 

As  it  is  currently  available,  stability  theory  of 
dynamical  systems  requires  an  extensive  background  in  higher 
mathematics.   The  purpose  of  this  thesis  then  is  to  make 
this  theory  accessible  to  the  reader  who  does  not  have  this 
background.   It  is  assumed  that  the  reader  has  taken  courses 
in  elementary  differential  equations,  intermediate  mathe- 
matical analysis,  and  linear  algebra.   All  the  topological 
properties  needed  are  described  in  Section  1. 

Sections  1-3  deal  with  developing  basic  concepts  of 
dynamical  systems.   Sections  4-6  use  these  concepts  in 
studying  various  kinds  of  stable  or  unstable  systems.   Topics 
of  recursiveness  and  dispersiveness  are  also  discussed. 

Among  the  stable  systems,  Poisson  and  Liapunov  stability 
will  be  discussed  as  well  as  non-wandering  systems.   Poisson 
Instability,  Lagrange  instability,  and  complete  instability 
will  also  be  investigated.   Topics  in  dispersiveness  include 
parallelizable  and  dispersive  systems.   The  relationship 
between  several  of  these  systems  will  be  shown.   One  such 
relationship  to  be  shown  is  that  a  parallelizable  system  is 
dispersive,  completely  unstable,  Lagrange  unstable  and 
Poisson  unstable. 

Several  important  and  revealing  theorems  will  be  pre- 
sented.  Among  these  are  the  theorems  involving  the  Liapunov 
function  which  allows  analysis  of  a  system  as  to  stability 
or  instability  without  actually  solving  the  system. 


The  book  Stability  Theory  of  Dynamical  Systems  by 
Bhatla  and  Szego  will  be  used  extensively  throughout  this 
thesis  and  its  bibliography  contains  an  excellent  list  of 
references. 


I.   TOPOLOGICAL  PROPERTIES 

The  purpose  of  this  section  is  to  expose  the  reader  to 
all  of  the  topological  properties  needed  to  understand  the 
material  presented  in  this  thesis.   It  is  important  that 
the  reader  develop  a  firm  understanding  and  familiarity 
with  these  concepts  to  facilitate  in  understanding  of  the 
material  to  follow  in  later  sections. 

We  begin  with  the  concept  of  distance. 

If  X  is  a  non-empty  set,  a  metric  or  distance  function, 
on  X  is  a  real-valued  function  d  of  ordered  pairs  of  elements 
of  X  which  satisfy  three  conditions.   First,  d(x,y)  _>  0, 
and  d(x,y)  =  0  if  and  only  if  x  =  y.   Second,  d(x,y)  =  d(y,x). 
This  is  called  the  symmetric  property.   Finally, 
d(x,y)  <_  d(x,z)  +  d(z,y).   This  property  is  called  the 
triangle  inequality.   A  metric  space  X  is  then  a  set  X  with 
a  metric  defined  on  it. 

If  Xq  is  any  point  of  a  metric  space  and  r  is  a  positive 
real  number,  the  open  sphere  S(Xp^,r)  with  center  Xq  and 
radius  r  is  the  subset  of  X  defined  by  S(x^,r)  =  {x :d(x,XQ) <r} 
Now  a  subset  G  of  X  is  said  to  be  open  if  and  only  if  it  is 
the  union  of  open  spheres.   Since  x^^  is  clearly  an  element 
of  S(x^,r),  we  say  that  S(Xp^,r)  is  a  (spherical) 
neighborhood  of  x^. 

If  A  is  a  subset  of  a  metric  space  X,  a  point  x  in  X  is 
called  an  accumulation  point  of  A  if  each  open  sphere 


centered  on  x  contains  at  least  one  point  of  A  different 
from  X.   The  essential  idea  here  is  that  the  points  of  A 
different  from  x  get  arbitrarily  close  to  x.   Then  a  subset 
F  of  a  metric  space  X  is  called  a  closed  set  if  it  contains 
all  its  accumulation  points.   Another  characterization  of 
a  closed  set  is  that  a  set  F  C  X  is_  closed  if  and  only  if 
its  complement  in  X  is_  an  open  set .   Then  by  definition, 
the  whole  space  X  and  the  empty  set  0  are  open  sets.   Then 
since  the  empty  set  is  the  complement  of  the  whole  space, 
we  have  by  our  previous  observation  that  the  whole  space  X 
and  the  empty  set  0  are  also  closed  sets. 

Now,  in  any  metric  space,  the  union  of  any  collection 
of  open  sets  is  itself  an  open  set  and  the  finite  inter- 
section of  open  sets  is  an  open  set.   Taking  set  complements 
it  follows  that  any  intersection  of  closed  sets  is  a  closed 
set  and  the  finite  union  of  closed  sets  is  a  closed  set. 

If  x-^  is  a  point  in  the  metric  space  X  andr  is  a  non- 
negative  real  number,  then  the  closed  sphere  S[X|-.,r]  with 
center  x^  and  radius  r  is  the  subset  of  X  defined  by 
S[xQ,r]  =  {x:d(x,XQ)  <.r}. 

Let  X  be  an  arbitrary  metric  space  and  let  A  be  a  subset 
of  X.   A  point  in  A  is  called  an  interior  point  of  A  if  it 
is  the  center  of  some  open  sphere  contained  in  A;  and  the 
interior  of  A,  denoted  by  int„A,  is  the  set  of  all  its 
interior  points.   Symbolically,  we  have 

intyA  =  {x:x  e  A   and   S(x,r)  C  A  for  some  r} .   Since  the 
interior  of  a  set  is  a  union  of  open  spheres,  it  follows 


that  int^A  is  an  open  set.   Hence  we  have  the  characterization 
that  a  set  A  is_  open  if  and  only  if  A  =  int„A.   We  often 
write  int  A  in  place  of  int„A  whenever  no  confusion  will 
result . 

If  X  is  a  metric  space  and  A  is  a  subset  of  X,  the 
closure  of  A  in  X,  denoted  by  cl^A,  is  the  union  of  A  and 
the  set  of  all  its  accumulation  points.   Intuitively,  cl„A 
is  A  itself  together  with  all  other  points  in  X  arbitrarily 
close  to  A.   Since  it  is  clear  that  cl^A  contains  all  its 
accumulation  points,  cl„A  is  a  closed  set.   Thus  we  have 
an  additional  characterization  of  a  closed  set :   A  C  X  is 
closed  if  and  only  if.  A  =  cl„A.   We  often  write  cl  A  in 
place  of  cl„A. 

Again  let  X  be  a  metric  space  and  let  A  be  a  subset 
of  X.   We  will  use  the  notation  X-A  to  denote  the  complement 
of  A  in  X.   A  point  in  X  is  called  a  boundary  point  of  A 
if  each  open  sphere  centered  on  this  point  intersects  both 
A  and  its  complement  X-A.   Now  the  boundary  of  A,  denoted 
bdyA,  is  the  set  of  all  the  boundary  points  of  A.   It  is 
easily  shown  that  bdy  A  =  cl^A  Hcl^CX-A).   Moreover,  the 
boundary  of  A  is  a  closed  set,  and  a  set  A  C  X  is_  closed 
if  and  only  if  it  contains  all  of  its  boundary  points . 

Every  metric  space  has,  in  addition  to  the  properties 

already  stated,  two  important  "separation"  properties. 

First  the  Hausdorff  property:   if  x  and  y  are  two  distinct 

points  in  X,  then  there  exist  disjoint  open  sets  U   and  U 

X      y 
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such  that  X  e  U  and  y  e  U  .   Second,  the  normal  space 
property:   if  A  and  B  are  disjoint  closed  sets  of  X,  then 
there  exist  disjoint  open  sets  U.  and  U^  such  that  A  C  U. 
and  B  C  Ug. 

We  now  use  these  basic  definitions  to  describe  more 
advanced  ideas  in  topology.   We  are  interested  in  character- 
izing "compactness",  "completeness"  and  "connectedness". 

Let  X  be  any  metric  space.   A  class  (G. }  of  open  subsets 
of  X  is  said  to  be  an  open  cover  of  X  if  each  point  in  X 
belongs  to  at  least  one  G. ,  that  is  Vj  G.  =  X.   A  subclass 
of  an  open  cover  which  is  itself  an  open  cover  is  called  a 
subcover.   Then  a  compact  space  is  defined  to  be  a  metric 
space  in  which  every  open  cover  has  a  finite  subcover.   This 
becomes  a  familiar  definition  v;hen  one  recalls  the  Heine- 
Borel  Theorem  which  states  that  i_f  [a,b]  ls_  a  closed  interval 
of  real  numbers  and  if  (G  },  x  e  [a,b]  is_  an  open  cover  of 
[ajb],  then  there  exists  a  finite  number  of  the  G  's  whose 
union  contains  [a,b].   This  is  just  the  statement  that 
closed  intervals  on  the  real  line  are  compact  sets. 

A  much  used  property  of  a  compact  space  is  that  any 
closed  subspace  of  a  compact  space  is  itself  compact .  Another 
useful  property  is  that  any  infinite  sequence  in  a  compact 
set  contains  a  subsequence  which  converges  to  some  point 
in  the  compact  set ,   This  also  becomes  a  familiar  property 
when  one  recalls  the  Bolzano-Weierstrass  Theorem:   Every 
bounded  infinite  set  A  of  real  numbers  has  at  least  one  point 
of  accumulation.   This  theorem  is  just  an  application  of 


this  property  applied  to  the  real  line.   These  last  tv7o 

properties  are  very  important  and  will  be  used  extensively 

throughout  this  thesis. 

To  define  the  notion  of  a  complete  metric  space  v;e  need 

the  concept  of  convergence.   If  {x  }  is  a  sequence  in  X, 

we  say  that  (x  }  is  convergent  to  the  point  x  in  X  if  given 

any  e  >  0,  there  exists  a  positive  integer  N  such  that 

d(x  ,x)  <  e   whenever  n  >  N.   Then  in  particular  we  have 
n'  ^ 

d(x  ,x  )  <  G  whenever  n,m  >  N.   A  sequence  with  this  latter 
property  is  called  a  Cauchy  sequence .   Now,  not  every  Cauchy 
sequence  is  a  convergent  sequence.   For  example,  consider 
the  following  sequence  in  the  metric  space  of  rational 
numbers,  (3.,  3.1,  3.1^,  3-1^1,  3-1^15,  ...).   Clearly,  this 
is  a  Cauchy  sequence  but  it  is  not  convergent  because  it 
converges  to  it;  but  it  is  not  a  rational  number.   Motivated 
by  this  example  we  say  that  a  metric  space  is  complete  if 
every  Cauchy  sequence  in  X  is  convergent.   Completeness  is 
related  to  compactness  in  that  a  set  is  compact  if  and  only 
if  it  is  complete -and  totally  bounded.   By  totally  bounded 
we  mean  that  the  space  X  is  the  union  of  a  finite  number  of 
open  spheres  of  radius  less  than  r  for  each  positive  r. 
Finally,  we  define  the  concept  of  connectedness.   A 
connected  space  is  a  metric  space  X  which  cannot  be  repre- 
sented as  the  union  of  two  disjoint  non-empty  open  sets. 
If  X  =  A  U  B  where  A  and  B  are  disjoint  non-empty  open  sets, 
then  X-A  =  B  and  X-B  =  A  so  A  and  B  are  also  closed  sets. 
Thus  a  space  X  ls_  connected  if  and  only  if  the  only  sets 
which  are  both  open  and  closed  are  the  whole  space  X  and  the 
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empty  set  0.   This  property  is  very  important  and  it  can  be 
shown  that  the  intermediate  value  theorem  in  calculus  is 
equivalent  to  this  property  for  intervals  on  the  real  line. 

We  now  turn  our  attention  to  the  concept  of  a  continuous 
function.   This  concept  has  several  different  characteriza- 
tions and  we  begin  with  the  basic  definition  that  looks  like 
the  definition  used  in  elementary  calculus. 

Let  f  be  a  function  from  a  metric  space  X  with  metric  d 
to  a  metric  space  Y  with  metric  d'.   Then  f  is  said  to  be 
continuous  at_  x^  e  X,  if  for  each  e  >  0,  there  exists  a 
6  >  0  such  that  for  x  e  X,  we  have  d' (f (x) ,f (Xq) )  <  c 
whenever  d(x,Xp.)  <  6.   The  function  f  is  said  to  be  continuous 
if  it  is  continuous  at  each  point  of  X.   The  definition  we 
saw  in  calculus  was  for  real-valued  functions  so  the  metrics 
just  involved  the  usual  absolute  value  for  the  real  line. 

Several  other  characterizations  of  a  continuous  function 
will  be  needed  in  this  thesis  and  are  stated  here.   A  function 
is  continuous  if  and  only  if  the  inverse  image  of  each  open 
(closed)  set  is  an  open  (closed)  set .   Also  a  function  f  is 
continuous  if  and  only  if  for  each  sequence  {x  }  which 
converges  to  x,  the  sequence  {f(x  )}  converges  to  f(x). 

We  now  state  some  important  properties  of  continuous 
functions, 

A  continuous  real-valued  function  f 
defined  on  a  compact  set  A  attains  its 
maximum  and  minimum  values  on  A  in  the 
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following  sense :   if^  a  =  inf{f(x):x  e  A} 
and  b  =  sup{f(x):x  e  A},  then  there  exists 
x«,Xp  e  A  such  that  f(x-.)  =  a,  and 
fCx^)  =  b. 

The  continuous  image  of  a  compact  (connected) 
set  is  compact  ( connected) . 

Next,  we  define  the  concept  of  a  homeomorphism.   A 
homeomorphism  is  a  one-to-one  continuous  function  from  one 
metric  space  onto  another  metric  space  such  that  the  inverse 
of  this  function  is  also  continuous.   (A  one-to-one  function 
is  a  function  which  assigns  different  images  to  different 
elements  in  its  domain.   A  function  f  from  a  metric  space 
X  to  a  metric  space  Y  is  onto  if  the  image  of  X  under  f  is 
the  entire  space  Y.   If  a  function  is  one-to-one  and  onto, 
then  its  inverse  is  well  defined.) 

We  have  already  stated  the  definition  of  a  compact  space. 
A  closely  related  idea  is  that  of  a  locally  compact  space. 
As  its  name  suggests,  a  metric  space  X  is  locally  compact 
if  and  only  if  each  point  of  X  has  at  least  one  compact 
neighborhood.   Hence,  a  compact  space  is  automatically 
locally  compact. 

One  of  the  most  important  properties  of  a  locally  compact 
space  is  that  through  the  addition  of  one  point,  usually 
referred  to  as  the  ideal  point,  the  space  can  be  extended 
to  a  compact  space.   This  extension  is  known  as  the  one-point 
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compactification.   The  one-point  compact If Icat Ion  of  a 
metric  space  X  is  the  space  X'  =  X  U  (w)  where  the  open  sets 
in  X'  are  all  the  open  sets  in  X  together  with  all  subsets 
U  of  X'  such  that  X'-U  is  a  compact  subset  of  X.   In  other 
words,  a  set  U  is  open  in  X'  if  and  only  if  U  r\  X  is  open 
and  whenever  o)  £  U,  X-U  is  compact.   As  this  is  an  important 
fact,  we  give  a  formal  proof.   V/e  must  show  that  X'  is  compact 
and  that  unions  and  finite  intersections  of  open  sets  are 
open. 

THEOREM:   The  one-point  compactification  X'  of  a  metric 
space  X  is_  compact . 

PROOF:   Finite  intersections  and  arbitrary  unions  of  open 
sets  in  X'  intersect  X  in  open  sets.   If  w  is  a  member  of 
the  intersection  of  two  open  subsets  of  X',  then  the  comple- 
ment of  the  intersection  is  the  union  of  two  closed  compact 
subsets  of  X  and  is  therefore  closed  and  compact.   Hence 
the  intersection  of  two  open  subsets  is  open  in  X'  .   If  to 
belongs  to  the  union  of  the  members  of  a  family  of  open 
subsets  of  X',  then  to  belongs  to  some  member  U  of  the  family. 
Moreover  the  complement  of  the  union  is  a  closed  subset  of 
the  compact  set  X-U  and  is  therefore  closed  and  compact. 
Hence  the  union  of  open  sets  is  open  in  X' .   Consequently, 
arbitrary  unions  and  finite  intersections  of  open  sets  in 
X'  are  open.   Next,  let 'U.  be  an  open  covering  of  X'.   Then 
to  is  a  member  of  some  U  in  "U  and  X-U  is  compact.   Hence 
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there  is  a  finite  subcover  of  I*-  which  covers  X-U  and  thus 
this  finite  subcover  together  with  U  is  a  finite  subcover 
of  X'.   Thus  X'  is  compact.  // 

Our  final  topic  in  this  section  of  topological  prereq- 
uisites is  the  concept  of  product  spaces.   A  finite  product 
space  is  the  cartesian  product  of  two  metric  spaces.   In 
this  thesis,  one  of  the  metric  spaces  will  usually  be  the 
real  line  R;  consequently  if  X  is  a  metric  space,  we  have 
the  product  space  to  be  XxR.   Hence  any  point  in  the  product 
space  will  be  written  (x,t)  where  x  e  X  and  t  e  R.   Now  the 
open  sets  in  the  product  space  are  just  the  cartesian  products 
of  open  sets  in  the  metric  spaces.   That  is  if  U  is  an  open 
set  in  X  and  V  is  an  open  set  in  R,  then  UxV  is  an  open  set 
in  XxR,  and  every  open  set  in  XxR  is  of  that  form.   If  d  is 
the  metric  for  X,  we  define  the  metric  d'  on  XxR  in  the 
following  manner.   If  x,y  e  XxR,  then  there  are  elements 
a,b  e  X  and  r,s  e  R  such  that  x  =  (a,r)  and  y  =  (b,s).   Then 
d'(x,y)  =  d(a,b)  +  |r  -  s|.   The  spaces  X  and  R  are  called 
the  coordinate  spaces  of  the  product  and  the  functions  P.. 
and  P„,  which  carry  the  point  (x,t)  of  XxR  into  x  and  t 
respectively,  are  called  the  projections  into  the  coordinate 
spaces.   It  is  easily  shown  that  the  projection  mappings 
are  continuous  functions  and  that  their  inverses  are  also 
continuous  functions.   Moreover  f:Y  -►  X  is  continuous  if 
and  only  if  each  of  the  compositions  Pyof  and  PyOf  is 
continuous,  for  any  metric  space  Y. 
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A  function  f  =  f(x,t)  defined  on  the  product  space  is 
continuous  in  x  if  and  only  if  for  each  t  the  function  f(  ,t) 
whose  value  at  x  is  f(x,t),  is  continuous.   Similarly, 
f(x,t)  is  continuous  in  t  if  and  only  if  for  each  x  e  X, 
the  function  f(x,  )  such  that  f(x,  )(t)  =  f(x,t),  is  con- 
tinuous.  Finally  if  f  is  continuous  on  the  product  space 
then  f  is  continuous  in  both  x  and  t  by  our  previous  remarks. 
However,  the  converse  is  not,  in  general,  true. 
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II.   DYNAMICAL  SYSTEMS 

We  seek  to  motivate  the  concept  of  dynamical  systems 
by  some  examples  of  differential  equations.   But  first  we 
need  the  existence  and  uniqueness  theorem  associated  with 
such  systems. 

We  consider  systems  of  first-order  differential  equations 
written  in  the  form 


dx, 

dt"  "  f*i(t;x^,...,x^) 


dx 


n  _ 


IvL^JX-j  ft   •  •  )  X  )   • 


dt     ^^v^,^^-j^,.  ..  ,^^^ 

Such  systems  frequently  arise  out  of  physical  problems  where 
the  independent  variable  t  represents  time.   Then  a  solution 
of  the  above  system  describes  the  "trajectory"  of  a  particle 
moving  in  n-space.   To  see  this  more  clearly,  we  look  at 
the  following  two  examples. 

Consider  the  following  system  defined  in  R  : 


dXp 

-rr—   =   cost 
dt 
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dF"  "  ^• 


Then  the  general  solution  to  this  system  is  given  by 


x^  =  cost  +  c. 


Xp  =  sint  +  Cp 


X-.  =  t  +  Ct 


where  the  c's  are  all  constants  of  integration. 
Suppose  we  specify  that  this  solution  passes  through  the 
point  X,  =  1,  Xp  =  0,  and  x^  =  0  at  time  t  =  0.   Then  the 
solution  is  given  by 


X,  =  cost 


Xp  =  sint 


X3  =  t. 


Such  a  solution  describes  the  "trajectory"  sketched  below, 
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It  is  clear  that  at  any  time,  we  may  find  the  exact 
position  of  the  particle.   The  same  is  true  of  the  following 
example . 

Consider  the  system  defined  in  R   given  by 


dx. 


dt 
dx. 


dt 
dx 


=  2 


=  2t 


dt 


3  =  t2 


where  we  specify  that  its  solution  passes  through  the  point 
x^  =  0,  x^  =  0,  and  x^  =  0  at  time  t  =  0,  Then  this  system 
has  the  solution  given  by 


x^  =  2t 
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Xp  =  t' 


^3    3 


Its  "trajectory"  is  sketched  below, 


Thus  an  Initial  value  problem  for  such  a  system  consists 
of  selecting  an  initial  time  t  =  t-^  and  a  point  X_  in  R  , 
and  then  seeking  a  solution  whose  graph  passes  through  X,. 
at  time  t  =  t^^. 

By  setting  X(t)  =  (x^  (t ) , .  .  .  ,Xj^(t ) )  , 

f(t,X)  =  (f^(t,X),...,f^(t,X)), 
this  system  may  be  rewritten  as 


11=  f(t,X),   where   ||=  (x^ ' (t ) , . . . ,x^' (t ) ) . 


We  next  want  to  examine  conditions  under  which  unique 
solutions  exist  for  an  initial  value  problem  of  the  system. 
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Suppose  that  f  =  f(t,X)  is  bounded  and  continuous  in  a  region 
D  of  R    .   Furthermore,  assume  that  f  satisfies  a  Lipschitz 
condition  in  D;  that  is  there  exists  a  non-negative  real 
number  b  such  that  whenever  (t,X-,)  and  (t,Xp)  are  any  two 
points  in  R,  d(f(t,X-j^),  fCtjX^))  <  bdCX^jX^)  where  d  is 
the  ordinary  Euclidean  metric.   Then  we  have  the  following 
theorem. 


Let  f  =  f(t,X)  be_  continuous,  bounded,  and 
satisfy  a  Lipschitz  condition  in  a  region  D 
of  R    ,  and  let  (t«,Xj-^)  be_  any  point  in  D. 
Then  the  initial  value  problem 

xctg)  =  Xq 

has  a  unique  solution  on  an  interval  I  of  R 
and  the  graph  of  the  solution  lies  in  D. 

We  will  be  Interested  mainly  in  systems  defined  on  all 
of  R.   Then  by  ensuring  that  f  is  continuous  and  bounded  on 
all  of  R    and  satisfies  a  global  Lipschitz  condition,  we 
will  have  existence  and  uniqueness  of  solutions  to  the  above 
system. 

One  elegant  and  simple  proof  of  the  above  theorem  requires 

a  knowledge  of  contraction  mappings  and  some  theory  of 

fixed  points  and  will  not  be  presented  here  [Kreider,  Kuller, 

Ostberg,  p.  385]. 
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A  classical  proof  may  be  found  in  Theory  of  Ordinary 
Differential  Equations  by  Coddington  and  Levinson. 

Of  critical  importance  is  the  question:   what  happens 
to  the  solution  if  small  changes  are  made  in  the  function  f 
or  in  the  initial  conditions  t^  and  X^?   In  most  physical 
problems  whose  solutions  are  obtained  from  such  a  system, 
experimental  error  and  estimations  of  physical  constants 
may  affect  the  initial  conditions,  and  the  presence  of 
variable  parameters,  such  as  temperature  or  density,  may 
affect  f.   We  would  hope  that  such  small  variations  would 
not  cause  a  drastic  variation  in  the  solution.   It  can  be 
shown  that  these  unique  solutions  guaranteed  by  the  theorem 
depend  continuously  on  f  and  on  t^  and  X^  for  a  certain 
amount  of  time;  that  is  small  changes  in  f  or  in  t^  and  X^ 
produce  small  changes  in  the  solution.   More  precisely. 

Let  f  =  f(t,X)  be_  continuous ,  bounded,  and 
satisfy  a  Lipschitz  condition  in  a  domain  D 
of  the  (n  +  l)-dimensional  (t,X)  space ,  and 
suppose  X  is_  a  solution  of  the  system  on  an 
interval  I :  a  <_  t  £  b .   There  exists  a  6  >  0 
such  that  for  any  (t,  i^r^)    where  a  <  t,  <  b, 
and  d(X-.  ,X(t-|  ) )  <  6  there  exists  a  unique 
solution  Y  on^  I  with  Y(tp,t^ ,X^)  =  Xq  where 
a  <  tp  <  b.   Moreover,  Y  is_  continuous  in  D. 

(Proofs  of  the  continuity  of  solutions  may  be  found  in 
Kreider,  Kuller,  Ostberg,  p.  390  and  Coddington,  Levinson, 

p.  22.) 
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What  happens  to  these  solutions  over  an  extended  period 
of  time  involves  the  theory  of  stability  which  is  the  primary 
focus  of  this  thesis. 

We  will  be  interested  mainly  in  first-order  autonomous 
systems,  that  is,  "time  independent"  systems.   Such  a  system 
may  be  written  in  the  form 


X  =  f(X)   v/here  X  denotes 


xct^)  =  Xq. 


dX 
dt 


The  follovjing  examples  provide  some  insight  into  such 
autonomous  systems. 

EXAMPLE  2.1   Consider  a  simple  pendulum  of  mass  m  swinging 
without  friction  on  a  weightless  rod  of  length  1.   By 
Newton's  second  law,  the  equation  of  motion  of  the  pendulum 
is  given  by 


-mgsin  0  =  ml 


dt2 


or 


6  =  -g/1  sin  e 

2 

By  setting  g/1  =  k  and  6  =  oj  we 

obtain  the  following  system: 
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-mgs'iAQ 


t 


ftii'fion 


e  =  cj 


CO  =  -k   sin  6. 


It  can  be  shown  [Kreider  ,  Kuller,  Ostberg,  p.  397]  that 
the  general  solution  to  this  system  is 


2     2 

CO  =  2k   cos  0  +  c 


where  c  >_  -2k  is  a  constant  depending  on  the  initial 
conditions. 

Graphing  several  of  these  solutions  in  the  6to-plane  we 
have  the  following  diagram. 


FIGURE  2.1 


To  see  how  this  diagram  was  obtained,  we  see  that  w  is 

2 

the  angular  velocity.   Suppose  c  =  -2k  .   Then  when  9=0, 
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(jo  =  0  and  hence  6  and  co  are  both  zero.   Thus  the  system  is 
not  changing  with  time.   The  point  (^)  is  the  complete 
trajectory  for  this  case.   Physically  this  shov;s  that  if 
the  mass  is  started  at  the  equilibrium  position  v/ith  no 

angular  velocity,  the  mass  will  remain  at  the  equilibrium 

2 
position  indefinitely.   Suppose  c  =  2k   and  the  system  is 

started  at  6  =  tt  .   Then  again  to  =  0  and  hence  6  and  w  are 

both  unchanging  with  time  and  the  point  (^)    describes  the 

complete  trajectory.   Physically  this  means  that  if  the 

mass  is  balanced  vertically  above  the  equilibrium  point 

with  no  angular  velocity,  it  will  also  remain  in  that  position 

indefinitely.   Clearly  the  same  analysis  can  be  made  for 

9  =  tniT  which  is  exactly  what  we  would  expect.   If  this 

system  is  started  at  any  other  angle  G,  it  will  have  an 

initial  velocity  which  will  cause  the  mass  to  approach  the 

angle  9  =it   or  9  =  -tt  and  it  will  get  arbitrarily  close 

but  will  never  actually  reach  that  point.   Neither  will  this 

mass  swing  back  away  from  this  point  but  rather  it  will 

continue  to  approach  the  angle  9  =  it  or  9  =  -n  as  time  goes 

to  infinity.   This  is  the  physical  meaning  of  the  trajectories 

"•IT  TT 

between  the  points(^  )  and  (r^)' 

2  2 

Now  suppose  -2k   <  c  <  2k  .   From  the  previous  remarks 

it  is  clear  that  the  mass  will  have  a  velocity  at  9  =  0  but 
its  velocity  is  not  great  enough  to  permit  it  to  approach 
the  angle  9  =  tt  or  9  =  -tt.   We  would  expect  in  this  case 
that  the  mass  would  swing  back  and  forth  across  the  equili- 
brium point  always  reaching  the  same  angle  on  either  side 

2l\ 


of  the  equilibrium  point.   This  is  the  physical  interpre- 
tation of  the  elliptical  trajectories  shown  in  the  diagram. 

2         2 
Finally,  suppose  c  >  2k  .   Then  w  is  always  positive 

which  means  that  w  is  always  positive  or  always  negative. 
If  (I)  is  always  positive,  then  6  is  always  positive  and 
hence  6  is  always  increasing.   A  similar  statement  follows 
if  CO  is  always  negative.   Physically  this  is  the  case  when 
the  mass  has  been  given  enough  initial  velocity  to  swing 
completely  over  the  top  position.   In  a  frictionless  envir- 
onment we  would  expect  this  behavior  to  continue  indefinitely. 
Such  is  the  physical  interpretation  of  the  top  and  bottom 
trajectories  in  the  diagram. 

To  see  the  direction  of  a  particle  along  any  trajectory, 

« 

pick  an  angle  8.   For  0  <  9  <  tt,  we  see  that  to  <  0  so  to  is 
decreasing.   For  -tt  <  6  <  0,  we  see  that  to  >  0  so  to  is 
increasing.   This  completes  our  analysis  of  the  diagram. 

It  is  noted  that  in  this  example  the  solution  is  not 
given  as  a  function  of  time.   To  obtain  such  would  mean  the 
difficult  task  of  solving  a  non-linear  differential  equation. 
This  is  not  any  more  informative  since  the  movement  of  a 
point  along  any  trajectory  with  increasing  time  can  easily 
be  determined  from  the  system  describing  the  motion  as  we 
have  just  observed. 

A  more  elementary  approach  to  the  problem  is  by  observing 
that  for  small  angles  0,  sin  6  is  approximated  by  6.   Thus 
our  system  would  be  given  by 
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9  =  0) 


0)  =  -k  e. 


It  can  be  shown  that  the  solution  to  this  system  is 


e  =  Asln  kt  +  Bcos  kt 


0)  =  Akcos  kt  -  Bksin  kt 

where  A  and  B  are  constants  depending  on  the  initial  conditions 
[Spiegel,  p.  210].   It  is  seen  that  this  solution  describes 

simple  harmonic  motion  and  closely  approximates  our  trajec- 

2  2 

torles  for  small  c  where  -2k   <  c  <  2k  . 

EXAMPLE  2.2  We  consider  the  following  vibrating  spring 
problem. 

yy////////////// 


I 


Equilibrium              o 
position- , Z? 


W 


J 


The  equation  of  motion  is  given  by: 


g 


+  by  +  ky  =  0 
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where  W  is  the  weight  of  the  mass,  g  is  the  force  of  gravity, 

b  is  the  damping  constant  and  k  is  the  proportionality 

constant  of  the  spring.   Suppose  W  =  6  lbs,  g  =  32  lbs, 

k  =  12  and  b  >_  0  is  left  as  an  undetermined  constant  for 

now. 

Then  it  can  be  shown  [Spiegel,  p.  195]  that  the  general 

solution  for  y  is 

8, 


rbt 

3  ^-^       '  '3 


y  =  e  ^   (Acos  (|  J  9  -  b^  t)  +  Bsin  (|J9  -  b^  t)) 


where  the  constants  A  and  B  depend  on  the  initial  conditions. 

If  V7e  specify  the  initial  conditions  to  be  y(0)  =  1/3 
and  y(0)  =  0,  v/e  have  the  solution, 

_8^^^         

y  =  {he    ^   (3  cos(|  J9  -  b^  t)  +   |,.^-  ,,sin(|j9-b^t)) 
^  ^  J  9-b'^    ^ 

To  make  this  example  more  explicit  we  will  consider  two 
cases  for  b.   First  an  undamped  system  v;ith  b  =  0  and  then 
a  damped  system  with  b  =  1.5.   With  b  =  0,  we  have 


y  =  5"  cos   8t 


With  b  =  1.5,  we  have 


y 


1   -H 


=  ^  e"^^(3cos  4  Jl  t  +   n  sin  4  JT  t) 


By  use  of  a  trigonometric  identity,  this  last  equation  may 
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be  written 


y  =  Lpe-^t  sin  (^Jl  t  +  1) 
^  3 


Now  that  we  have  the  solutions,  we  change  this  second- 
order  problem  into  a  system  of  first-order  equations.  Let 
X  =  (^)    and  y  =  X.   Then  for  b  =  0 

V 

w  -^ 
X  =  f(X)  =  (        ) 

X 

X(0)  =  (°/3) 


For  notation,  let  the  initial  point  (-i/o)  ~  ^n*  ^^^^   the 
solution  to  this  system  is 

-|sin  8t 
X(XQ,t)  =  (  ^      ). 

^os  8t 


For  b  =  1.5>  we  have 


-bgx    -kgy 
X  =  f (X)  =  /   W   "   W 


X(0)  =  (J/^) 


and  the  solution  to  this  system  is 
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.  |e-^^cos(iint+^)  -  |^3e"^^sin(^^F3t+^) 
X(X^,t)  =3  3    3  3 


It  is  clear  that  as  t  -f  «»,  X(X^,t)  ->   (^)  which  we  would 
intuitively  expect  from  a  damped  spring. 

We  sketch  these  solutions  on  the  x,y-plane. 

y  .  y 


FIGURE  2.2 

For  the  undamped  problem  we  see  that  its  trajectory 
describes  simple  harmonic  motion.   It  continues  to  oscillate 
indefinitely  with  the  same  maximum  displacement  from  the 
equilibrium  position.   On  the  other  hand,  the  displacement 
and  its  rate  of  change  both  decrease  with  increasing  time 
for  the  damped  problem. 

Since  we  have  two  different  systems,  it  is  necessary  to 
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plot  their  solutions  on  two  different  planes.   This  helps 
to  avoid  any  misunderstanding  in  our  future  work. 

EXAMPLE  2.3   Consider  the  second-order  differential  equation 
y  +  y  =  0,  where  y  is  considered  as  a  function  of  time. 
By  setting  x  =  y,  this  equation  may  be  converted  to  the 
system 

X  =  -y 
y  =  X 

Now  considering  X  =  (  )  with  the  initial  condition  that  at 
t  =  0,  the  solution  passes  through  the  point  (q)j  we  have 
the  following  initial  value  problem 

X  =  f(X) 
X(0)  =  (J). 

where  f  is  the  linear  function  defined  by 

f    (t,x,y)  _ 

f(t,x,y)    =    (^    /,  n)    =    (      ).      A   simple    check   shows   that 

ipVi^jX,yj  X 

this  system  has  the  solution 


x(t)  =  (llll). 


If  we  let  X^  denote  the  point  (q),  then  X-^  denotes  the 
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Initial  condition.   Further  since  the  solution  depends  on 
this  initial  condition,  the  solution  X(t)  is  more  properly 
stated  X(XQ,t)  =  (11^  I). 

There  are  three  observations  we  wish  to  make  about  this 

solution.   First,  we  see  that  X(X„,0)  =  X^ .   Next, 

X(X(XQ,t2),t^)  =  X(XQ,t-L  +  t^). 

cos  tp 
This  is  true  since  X(X(X^,tp) ,t  )  =  X( (  .   .  ),t^). 

The  term  on  the  right  is  the  solution  of  the  system  which 

cos  tp 

has  the  initial  condition  Y„  =  (  .   4_  ) .   Its  solution  is 

0    sin  tp 

cos  tp  cos  t  -  sin  tp  sin  t 

X(Y^,t)  =  (    ^    .   4.  ^   .   4.      ^).   Hence 
0'      cos  tp  sm  t  +  sm  tp  cos  t 

cos  tp  cos  t,  -  sin  tp  sin  t^ 
^^^O'h^  =  ^cos  t2  sin  t^  +  sin  t^  cos  t.^^*   ^^^^^  a 

trigonometric  identity  we  have 


cos  (t^  +  tp) 
^^^0>h)  =  ^sin  (tj  +  t^)^-   T^-^ 


X(X(XQ,t2),t^;  =  X(YQ,t^)  =  X(XQ,t-^  +  t^).   Finally  X(XQ,t) 
is  continuous  because  both  of  its  components  are  continuous. 

It  is  clear  that  the  graph  of  the  solution  is  a  circle 
in  the  x,y-plane,  and  if  we  consider  the  graph  to  be  the 
trajectory  of  the  point  X_^ ,  that  point  moves  counter-clockwise 
with  increasing  time.   This  can  be  seen  by  choosing  some 
positive  value  for  x.   Then  y  is  positive  so  y  is  increasing. 
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FIGURE  2.3 

If  we  alter  the  initial  condition  to  X^  =  (,  ),  the 

a  cos  t  -  b  sin  t 
solution  is  X(X^,t)  =  (    •   +-  .  h     +- ^  ^^^   ^^^   trajectory 

is  a  circle  passing  through  the  point  X_. 

Now  referring  to  Example  2 . 3  we  see  that  the  solution 

2 
X(X_,t)  is  a  map  from  the  product  space  R  xR  into  the  space 

2  2 

R  .   We  note  also  that  R  is  a  metric  space.   The  three 

observations  made  about  the  solution  are  precisely  the 
properties  defining  a  dynamical  system.   The  fact  that  the 
system  in  Example  2.3  does  define  a  dynamical  system  will 
be  made  clear  by  the  following  definition  where  in  this  case, 

the  map  it  is  defined  by  7T(x,t)  =  X(x,t)  where  x  denotes  a 

2 
point  in  R  .   In  all  further  work,  X  denotes  an  arbitrary 

metric  space. 

DEFINITION  2.^   A  dynamical  system  on  X  is  the  triplet 
(XjR,!!),  where  tt  is  a  map  from  the  product  space  XxR  into 
the  space  X  satisfying  the  following  axioms : 
1.   For  every  x  in  X,  ttCxjO)  =  x. 
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2.  For  every  x  in  X  and  t, jt^  in  R, 
TT(Tr(x,t  )  ,t2)  =  irCxjt^  +  tp),  and 

3.  The  map  tt  is  continuous. 

The  space  X  and  the  map  tt  are  called  the  phase  space  and 
the  phase  map,  respectively. 

From  this  point  on,  whenever  no  confusion  will  result, 
the  symbol  tt  will  be  suppressed.   Thus  the  image  under  ir 
of  the  point  (x,t)  will  be  written  more  simply  as  xt .   The 
first  two  axioms  then  become 

1.  For  every  x  in  X,  xO  =  x, 

2.  For  every  x  in  X  and  t,  ,t_  in  R, 
xt^Ct^)  =  x(t^  +  t^). 

Similarly,  if  M  d  X  and  A <C  R,  MA  is  the  set 
{xt:  X  e  M  and  t  e  A}. 

The  phase  map  determines  two  other  maps  when  either  x  or 
t  is  fixed.   For  fixed  t  in  R  the  map  tt  :X  ->•  X  defined  by 
TT  (x)  =  xt  is  called  a  transition;   for  fixed  x  in  X  the 
map  IT  :R  ->  X  defined  by  tt  (t)  =  xt  is  called  a  motion 
through  X. 

THEOREM  2. 5  For  each  t  iii  R,  it   is_  a  homeomorphism  on  X  onto 
itself. 

PROOF  For  any  t  in  R,  it   is  clearly  continuous  because  tt 
is  continuous.   tt   is  injective  since  xt  =  yt  implies  x  =  y 
because  x  =  xO  =  x(t  -  t)  =  xt(-t)  =  yt(-t)  =  y(t  -  t)  =  yO 
=y.   To  see  that  tt   is  surjective,  let  y  e  X  be  arbitrary. 
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Then  it  (x)  =  y  where  x  =  y(-t).   Finally  to  show  that  tt 
has  a  continuous  inverse,  it  is  sufficient  to  show  that 

IT    is  the  inverse  of  tt  .   To  see  this  we  note  that  if  it 

s  t   ^ 

and  7T   are  two  transitions,  the  composition  it  ©tt   is  the 

t  +  s 
transition  tt     :  for  any  x  in  X, 


IT  OTT^(x)  =  TT  (tt  (x))  =  TT  (xs)  =  XS(t)  =  x(s+t)  =  TT^    (x) 


It  is  also  clear  that  tt   is  the  identity  transition  because 

for  any  x  in  X,  tt  (x)  =  xO  =  x.   Now  since 

tt"  OTT   =  TT   ~   -   "^    i   the  transition  tt"   is  the  inverse  of 
t 

TT  . 


REMARK  2.6   For  any  x  in  X  and  [a,b]  C   R,  the  set  x[a,b] 
is  compact  and  connected. 

PROOF  The  set  {x}  is  compact  and  connected  in  X  and  [a,b] 
is  compact  and  connected  in  R.   Thus  {x}x[a,b]  is  compact 
and  connected.   Then  since  x[a,b]  is  the  continuous  image 
of  this  set,  x[a,b]  is  compact  and  connected.  // 

REMARK  2.7  Observe  that  the  transitions  tt  ,  t  in  R,  form  a 
commutative  group  with  the  group  operation  being  the  compo- 
sition of  transitions.  ^  • 

PROOF   The  closure  property  is  trivial  for  if  t,s  are  in  R, 

t   s    t  +  s 
t  +  s  is  in  R  and  tt  ott   =  tt      as  we  have  already  observed 
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A  J    ^j  t     ,    S       Vv  t       s+v  t+s+v 

Associative:      ir   0(77   ott    )    =   it   o-n-  =77 


t+S        V  /     t        Sn        V 

=    IT         Ott      =    Ctt    ott    ;o7t 


t      s  t      s  t 

Commutative:      it   077    (x)    =   tt    (tt    (x))    =   -n    (xs)    =    (xs)t 


x(s+t)    =    x(t  +  s)    =    (xt)s   =   7T^(xt) 


IT      (it      (x))     =     tt     ott      (x) 


Identity:     The  identity  is  tt   since  tt  ott   =  tt 

t    0+t    0   t 

=  TT   =  TT      =  TT  OTT 


Inverses:  The   inverse   of  tt      is   tt~      since 


t      -t  0 

TT     OTT  =     TT      . 


We  now  give  a  precise  definition  of  what  we  have  been 
calling  trajectories. 

2.7  DEFINITION  The  maps  y(x),  y    (x),  and  y~(x)  from  X  into 

2   are  defined,  for  any  x  in  X, 


y(x)  =  {xt:  t  e  R}, 


y"*"(x)  ='{xt:  t  e  r"^}. 


y"(x)  =  {xt:  t  e  R  }. 
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For  any  x  in  X,  the  sets  y(x),  y  (x),  and  y~ (x)    are  called 
respectively  the  trajectory ,   the  positive  semi-trajectory 
and  the  negative  semi-traj ect ory .   Note  that  for  any  x  in  X, 
Y(x)  =  xR. 

If  two  trajectories  y-i  ^^^   Yp  have  a  point  in  common, 
then  the  trajectories  are  identical.   This  fact  is  guaranteed 
by  the  existence  and  uniqueness  theorem.   For  if  x  e  y-i 
and  X  e  Yp >  this  theorem  states  that  there  exists  a  unique 
solution  passing  through  the  point  x.   Hence  y-,  ~  Yp*   This 
means  that  two  different  trajectories  can  never  cross  each 
other.   This  is  an  important  fact  and  is  very  useful  in 
diagramming  dynamical  systems. 

Another  example  of  a  dynamical  system  which  will  prove 
to  be  useful  is  a  dynamical  system  defined  on  a  torus. 
Consider  a  differential  system  in  the  plane 


f(X) 


where  f  satisfies  the  hypotheses  of  the  existence  and  unique- 

X  ^1 

ness  theorem.   Set  X  =  (  )  and  f  =  ( -,  )  and  assume  that  the 

y        12 

functions  f,  and  fp  are  periodic  with  period  1  in  each  of 
the  variables  x  and  y.   Then 


f^(x,y)  =  f^(x+l,y)  =  f^(x,y+l),   i  =  1,2 


We  take  the  square  {(x,y):  0  _<  x  <  1 ,  0  <_  y  <  1}  and 
identify  the  sides  x  =  0  and  x  =  1  and  the  sides  y  =  0  and 
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y  =  1.   Then  the  points  (0,0),  (0,1),  (1 ,0)  and  (1 ,1)  are 
identified.   Then  we  project  the  plane  onto  the  torus  by 
identifying  any  point  (x,y)  in  the  plane  with  a  point  (x,y) 
of  the  torus  where  x  =  x  (mod  1)  and  y  =  y  (mod  1).   Then 
the  trajectories  in  the  plane  projected  onto  the  torus  yield 
trajectories  of  a  dynamical  system  on  the  torus. 


rf}-6yvNiJ 


It  is  clear  from  Example  2.3  that  if  we  pick  any  point 
X  on  the  trajectory,  xt  is  also  on  the  trajectory  for  any  t 
in  R.   In  Example  2.2  if  we  pick  a  point  x  on  any  trajectory 
and  consider  the  positive  semi-trajectory ,  then  xt  belongs 
to  y''"(x)  for  any  t  e  R"^.   The  same  results  can  be  seen  for 
the  negative  semi-trajectory.   This  observation  leads  us 
naturally  into  the  following  concept  of  invariance. 

2.8   DEFINITION  A  set  M  C  X  is  called  invariant  whenever 

xt  e  M  for  all  x  e  M  and  t  e  R. 
It  is  called  positively  invariant  whenever  this  holds  with 
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R  replaced  by  R  and  is  called  negatively  invariant  if  it 
holds  with  R  replaced  by  R~. 


From  this  definition  we  see  that  a  set  M  C  X  is  invari- 
ant if  and  only  if  it  is  both  positively  and  negatively 
Invariant  for  if  M  is  invariant,  we  let  x  be  an  arbitrary 
element  of  M.   Then  xt  e  M  for  all  t  e  R  and  in  particular 
for  t  e  R  ,  xt  G  M  and  for  t  e  R~,  xt  e  M.   Thus  M  is  both 
positively  and  negatively  invariant.   If  M  is  both  positively 
and  negatively  invariant,  then  for  any  t  e  R  ,  xt  e  M  and 
for  any  t  e  R~,  xt  e  M.   Thus  xt  e  M  for  all  t  e  R  and  M 
is  invariant. 

This  concept  leads  us  to  some  interesting  and  important 
properties  associated  with  such  sets. 


THEOREM 


2 . 9   Let  {M. }  be  a  collection  of  positively  invari- 


ant ,  negatively  invariant  or  Invariant  subsets  of  X.   Then 
their  intersection  and  their  union  have  the  same  property. 
A  set  M  C  X  ls_  positively  invariant  if  and  only  if  the  set 
X-M  is_  negatively  invariant .   M  i_s  invariant  if  and  only 
If  X-M  ls_  invariant. 

PROOF  Let  the  sets  M.  be  positively  invariant.   Let  A  =VjM. 
and  B  =  Am  .   For  any  x  e  A,  x  e  M.  for  some  i.   Thus 
xt  e  M.  for  all  t  e  R  since  M.  is  positively  invariant. 
Then  xt  e  A  for  all  t  e  R  .   Hence  A  is  positively  invariant. 
Let  X  e  B.   Then  x  e  M.  for  all  i.   Then  since  each  M  is 
positively  invariant,  xt  e  M.  for  all  t  e  R  and  all  i. 
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Thus  xt  e  B  for  all  t  e  R  and  B  is  positively  invariant. 
The  proofs  for  the  negatively  invariant  and  invariant  cases 
are  entirely  analogous.   Let  M  be  positively  invariant. 
Let  X  e  X-M  and  t  e  R~  and  suppose  xt  e  M.   Since  -t  e  R  , 
xt(-t)  e  M.   But  xt(-t)  =  xO  =  X  which  implies  x  e  M  and 
we  have  a  contradiction.   Thus  xt  e  X-M  and  X-M  is  negatively 
invariant.   The  proofs  of  the  converse  and  the  third  part 
are  entirely  similar.  // 

THEOREM  2.10   Let  M  C  X  be  positively  invariant  3  negatively 
invariant  or  invariant.   Then  the  follov/ing  statements  are 
true. 

1.  The  closure  of  M,  cl  M,  and  its  interior,  int  M, 
have  the  same  property^ 

2.  Each  of  its  components  has  the  same  property. 
If_  M  is_  invariant ,  then  so  is  its  boundary  bdy  M. 

PROOF  1.   Consider  the  case  of  invariance.   Let  x  e  cl  M 

and  t  e  R.   Then  there  is  a  sequence  {x  }  in  M  such  that 

X  -*■  X.   By  the  invariance  of  M,  x  t  e  M  for  all  t  e  R 
n        ^  '   n 

and  each  n.   Then  by  continuity  we  have  that  x  t  ->-  xt .   Thus 
xt  e  cl  M  and  cl  M  is  positively  invariant.   The  proofs  for 
the  negatively  invariant  and  positively  invariant  cases 
are  entirely  analogous. 

Assume  M  is  positively  invariant.   Then  X-M  and 
hence  cl  (X-M)  are  negatively  invariant.   Consequently 
int  M  =  X  -  cl  (X-M)  is  positively  invariant.   The  proof  for 
M  negatively  invariant  is  entirely  similar. 
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2.   Let  M  be  invariant.   Now  M  =UM.  where  the  M.'s 
are  disjoint  components  of  M.   Suppose  M.  is  not  invariant 
for  some  i.   Then  there  exists  x  e  M.  and  t  e  R  such  that 
xt  g^  M  .   Now  M.  is  a  maximal  connected  subset  of  M  and 
x[0,t]  is  a  connected  set.   Thus  if  xt  e  M   for  some  j, 

J 

M.  U  M.  is  a  connected  set  and  M.  and  M.  are  not  disjoint. 
Thus  xt  i   M   for  all  j  v/hich  implies  that  xt  ^  M.   This 
contradicts  the  invariance  of  M  and  hence  each  component  of 
M  is  invariant.   For  the  converse  let  each  component  of  M 
be  invariant.   Then  since  M  =UM.  and  each  component  is 
invariant,  it  follows  that  M  is  invariant.   The  proofs  of 
the  other  cases  are  similar. 

Suppose  M  is  invariant,  then  so  is  X-M.   Consequently 
el  M  and  cl  (X-M)  are  invariant.   Hence  bdy  M  =  cl  M  H  cl  (X-M) 
is  invariant  and  int  M  =  X  -  cl  (X-M)  is  invariant.         // 

For  the  following  theorem  we  need  to  define  y(M),  y  (M) , 
and  y"'(M)  where  M  C  X.   We  define  these  as  we  would  expect; 
that  is  y(M)  =  {Mt:  t  e  R},  y"*"(M)  =  (Mt:  t  e  r"*"}  and 
similarly  for  y"*(M). 

THEOREM  2.11   For  any  x  e  X,  the  sets  y(x),  y'''(x),  and 
Y  (x)  are ,  respectively ,  invariant ,  positively  invariant 
and  negatively  invariant .   A  set  M  d  X  is_  invariant ,  posi- 
tively invariant ,  or  negatively  invariant  if  and  only  if, 
respectively,  y(M)  =  M,  y  (M)  =  M,  or  Y~(M)  =  M. 
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PROOF  Let  xt^  e  Y~(x)  for  t^  e  R"  and  let  t^  e  R~.   Then 
xt^(t2)  =  x(t^  +  t^)  =  xt^  for  t^  =  t^  +  tp.   Since 
t^jt^  e  R  ,  we  have  t-,  e  R~  and  xt-.  c   y~{x).      Hence  Y~(x) 
is  negatively  invariant.   The  proof  of  the  other  cases  is 
entirely  analogous.   For  the  second  part,  let  M  be  negatively 
invariant.   Now  Y~(M)  =  {Mt :  t  e  R~}.   In  particular  for 
t  =  0,  Mt  =  M  and  MC  Y~(M).   Suppose  Y~(M)  ^  M.   Then  there 
exists  at,  e  R~  such  that  Mt,  <X   M.   Then  there  exists 
x  e  M  such  that  xt   ^  M  which  contradicts  the  negative 
invariance  of  M.   Thus  Y~(M)  C  M  and  Y~(M)  =  M.   Conversely, 
suppose  y"(M)  =  M.   Then  (Mt :  t  e  R")  =  M  and  Mt  C  M  for 
all  t  G  R~.   Thus  M  is  negatively  invariant.   The  proofs  of 
the  other  cases  are  similar. 

The  following  remarks  are  useful  in  understanding  this 
material  and  prove  some  useful  results. 

REMARKS  2.12 

1.  The  sets  X  and  0  are  invariant. 

PROOF  Let  X  e  X.   Then  xt  e  X  for  all  t  e  R  since  xt  ^  0. 
Thus  X  is  invariant.   Since  0  =  X-X,  it  follows  that  0  is 
invariant . 

2.  A  set  M  C  X  is  invariant,  positively  invariant, 
or  negatively  invariant  if  and  only  if  for  each  x  e  M, 
respectively,  y{x)   C  M,  y  (x)  CM,  or  y  (x)  C  M. 
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PROOF  Let  M  be  Invariant.  Then  by  Theorem  2.11  y{M)  =  M. 
Thus  {Mt:  t  e  R}  =  M.  Hence  for  any  x  e  M,  xt  e  M  for  all 
t  e  R.  But  {xt:  t  e  R}  =  y(x).  Thus  y(x)  C  M.  Conversely, 
suppose  y(x)  c  M  for  each  x  e  M.  Then  {xt:  t  e  R)  d  M  and 
xt  e  M  for  all  t  e  R  and  all  x  e  M.  Hence  M  Is  invariant. 
The  proofs  of  the  other  cases  are  similar. 

Consider  once  again  Example  2.1,   For  the  points  („  ) 

2 
n  =  ±1,  ±2,  ...  we  see  that  6  =  o)  =  0  and  w  =  -k   sin  6  =  0. 

This  shov7S  that  each  x  =  (^  )  does  not  change  position  with 

increasing  or  decreasing  time.   Such  points  are  called 

critical  or  rest  points.   VJe  now  state  the  formal  definition. 

DEFINITION  2.13   A  point  x  e  X  is  said  to  be  a  critical 
point  if  X  =  xt  for  all  t  e  R. 

REMARK  2.14   If  x  is  critical  and  y  ?^  x,  then  x  does  not 
belong  to  the  trajectory  Y(y). 

PROOF  Suppose  x  is  critical  and  y  5^  x.   Then  {x)  =  y(x) 

and  since  y(x)  is  invariant,  {x)  is  invariant.   Thus 

X-{x}  is  invariant.   Since  y  7^  x,  y  e  X-{x}  so  Y(y)  ^  X-{x}. 
Thus  X  ^  Y(y)-  // 

With  this  remark.  Figure  2.1  should  more  properly  be 
shown  as  follows : 
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The  next  theorem  contains  several  characterizations  of 
critical  points  but  first   the  follov^ing  lemma  will  be  found 
useful. 


LEMMA  2.15   If  X  e  X  and  x  =  xt  for  some  t  e  R,  then 
X  =  x(nt)  for  all' integers  n. 

PROOF   If  X  =  xt,  then  x(-t)  =  xt(-t)  =  x(t  -  t)  =  x  so  we 
need  only  consider  positive  integers.   We  complete  the 
proof  by  induction.   Clearly  if  x  =  xt  for  some  t  e  R  , 
then  X  =  x(lt).   Assume  x  =  x(nt).   Then 
X  =  xt(nt)  =  x(t  +  nt)  =  x((n  +  l)t). 


// 


THEOREM  2.16   Let  x  e  X.   Then  the  following  are  equivalent 

1.  The  point  x  is_  critical 

2.  The  singleton  {x}  is_  y(x) 
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3.  The  singleton  {x}  ±s_  y    (x) 

4.  The  singleton  {x}  is_  Y~(x) 

5.  The  singleton  (x)  is_  x[a,b]  for  some  a  <  b 

6.  There  Is  a  sequence  {t  },  t   >  0,  t   ^0  with 
—  _  — :i n  '   n    '   n     

X  =  xt   for  each  n. 
n 

PROOF   The  equivalence  of  statements  1.  with  2.,  3-,  and  H. 

are  trivial.   We  next  prove  the  equivalence  of  statemnts  1. 

and  6.   If  x  is  critical  then  statement  6  follovjs  trivially 

since  x  =  xt  for  all  t  e  R.   Conversely,  assume  there  is 

a  sequence  {t},t   >0,t   -^0  with  x  =  xt   for  each  n 
n  '   n    '   n  n 

and  let  t  e  R.   If  t  =  kt   for  some  integers  k  and  n,  then 

n  ^  ' 

xt  =  x(kt  )  =  X  by  the  lemma.   Otherwise  since  t   >  0  for 
n       ''  n 

all  integers  n,  t/t   e  R  and  is  located  between  two  integers 

k  and  k  +1.   Thus  k   <  t/t   <  k  +1  and 
n      n  n      n    n 

k  t   <  t  <  (k  +  l)t  .   Let  6  =  min  {(k  +  l)t   -  t,  t  -  k  t^}. 
nn         nn  n      n'      nn 

Since  t   ->■  0 ,  there  exists  an  integer  m  >  n  such  that 
n     * 

t   <  6/2.   As  above  there  exists  integers  k  and  k  +  1 
m  "^     m      m 

such  that  k  t   <  t  <  (k  +  l)t  .   Since  this  interval  has 
mm         m      m 

length  (k  +  l)t   -  k  t   =  t   <  6/2,  this  interval  must  be 
°     m      m    m  m    m 

contained  inside  the  previous  interval  since 

max    {(k     +l)t      -t,    t-kt}<   6/2.      Thus   we   have 
m  m  *  mm 


k  t      <  k  t      <   t    <    (k     +   l)t      <    (k     +   l)t    . 
nnmm  m  m  n  n 


Continuing  in  this  manner,  we  have  constructed  a  sequence 
{k  t  }  with  the  property  that  k  t   ->  t .   Now  by  the 
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continuity  axiom,  x(k  t  )  -^^  xt  and  since  x(k  t  )  =  x  for  each 

n,  by  the  lemma,  we  have  x  =  xt .   Thus  x  is  critical. 

To  prove  the  equivalence  of  statements  1.  and  5.  we 

first  note  that  if  x  is  critical,  then  5.  follows  trivially. 

On  the  other  hand,  suppose  x  =  x[a,b]  for  some  a  <  b  and 

let  t  e  R.   If  t  =  kt^  for  some  integer  k  and  t^  in  [a,b], 

then  xt  =  x(kt  )  =  x  by  the  lemma.   Otherwise,  since  [a,b] 

is  closed,  there  is  a  sequence  {t  }  in  [a.b]  with  t   ->  a. 

n       *         n 

Consider  the  sequence  {t   -  a}.   Clearly  (t   -  a)  >  0, 

n  n 

(t   -  a)  -^  0  and  x  =  xa  =  x(-la)  by  the  lemma.   Since 

t   e  [a,b],  X  =  xt  .   Combining  these  last  two  statements 
n     '   '       n  ^ 

we  have  x  =  xa  =  x(-a)  =  xt  (-a)  =  x(t   -  a).   Thus  the 

n         n 

hypotheses  of  statement  6.  are  satisfied  and  x  is  critical.  // 

From  this  theorem  we  know  that  if  we  follow  a  point  of 
any  trajectory  in  a  dynamical  system  and  discover  that  that 
point  is  at  rest  for  any  period  of  time,  then  the  point 
remains  at  rest  indefinitely. 

Continuing  to  seek  motivation  from  our  examples,  we 
again  refer. to  Example  2.3.   We  see  that  if  x  is  any  point 
on  the  trajectory,  x  =  X(X-,t)  for  some  t  e  R  but  also 
X  =  X(X  ,t  +  2n7T)  for  n  =  ±1,  ±2,  ...   Transfering  this  to 
our  notation  of  a  dynamical  system,  we  have  that 
X  =  X_t  =  Xj-^(t  +  2n7T).   As  we  would  expect,  any  point  with 
this  property  is  called  a  periodic  point.   Formally,  we 
state: 
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DEFINITION  2.17   A  point  x  e  X  is  said  to  be  periodic  if 
there  is  a  T  5^  0  such  that 

xt  =  x(t  +  T)   for  all  t  e  R. 

A  number  T  for  which  this  holds  is  called  a  period  of  x. 
If  X  is  a  periodic  point  then  both  the  motion  and  the 
trajectory  are  said  to  be  periodic. 

The  restriction  that  T  7^  0  is  necessary  since  it  is 
true  that  xt  =  x(t  +  0)  for  all  t  but  the  point  x  may  not 
be  periodic.   It  is  also  noted  that  any  critical  point  x 
is  periodic  since  every  T  e  R  is  a  period  for  x.   The 
following  theorems  and  remarks  prove  some  of  the  expected 
properties  of  these  points. 


THEOREM  2.18   If  {x  }  is_  a  sequence  of  periodic  points  with 

positive  periods  T  -^0  and  x  -^  x,  then  x  is  critical. 
L. t: n     n     ' 

PROOF  As  in  the  proof  of  Theorem  2.16,  given  any  t  e  R, 

there  are  integers  k  and  k  +  1  such  that 
^     n      n 


k  T   <  t  <  (k   +  1)T   =  k  T  +  T  . 
nn—       n      n    nn    n 


Since  T  ->■  0 ,  k  T   ->  t .   Then  x   =  x  T   since  x   is  periodic 

n*nn  nnn  n 

and  X  T  =  X  (k  T  )  by  lemma  2.15.   Combining  these  last 
n  n    n  n  n   "^ 

two  statements  we  have  x  =  x  T  =  x  (k^T  )  ->  x  t .   Now 

n    n  n    nnn  n 


US 


since  x   -»-  x  we  have  x  =  xt ,   As  t  e  R  was  arbitrary, 

X  =  xt  for  all  t  e  R.   Hence  x  is  critical.  // 

From  the  definition  we  can  see  that  a  point  x  c  X  is 
periodic  if  and  only  if  there  is_  a  T  7^  0  with  x  =  xT.  For 
if  X  is  periodic,  we  know  there  exists  a  T  7^  0  such  that 
xt  =  x(t  +  T)  for  all  t  e  R.   But  in  particular,  for  t  =  0, 
we  have  x  =  xT.   On  the  other  hand,  if  there  exists  a  T  7^  0 
such  that  X  =  xT  and  if  t  is  any  element  in  R,  then 
xt  =  xT(t)  =  x(T  +  t)  =  x(t  +  T).   Since  t  was  arbitrary, 
we  see  that  x  is  periodic. 

In  Example  2.3  we  noted  that  T  =  27t  was  a  period  of  x 

and  T  =  ^77  was  a  period  of  x.   Clearly  in  this  example, 

X  has  an  infinite  number  of  periods.   We  v;ould  expect  that 

x  would  have  a  smallest  non-negative  period,  which  we  will 

call  its  fundamental  period.   If  there  exists  an  infinite 

number  of  periods  for  a  point  x  in  a  dynamical  system  whose 

fundamental  period  is  0,  then  as  an  immediate  consequence  of 

the  last  theorem, ,we  know  that  x  is  critical.   This  is  true 

since  we  can  certainly  form  a  sequence  {T  }  of  the  periods 

such  that  T  ->-  0  and  since  the  sequence  {x  }  where  x  =  x 
n  ^        n         n 

for  all  n  has  the  property  that  x  -^  x,  the  hypotheses  of 
the  theorem  are  fulfilled.   The  following  theorem  proves 
the  existence  of  a  fundamental  period  and  also  proves  why 
we  need  only  consider  positive  periods. 

THEOREM  2.19   If  X  e  X  is  periodic  but  not  critical,  then 
there  is  a  T  >  0  such  that  T  is_  the  smallest  positive  period 


^7 


of  X.   Further  !£  T'  is  any  other  period  of  x,  then 
T'  =  nT  for  some  Integer  n. 

PROOF  Let  P  =  {t  >  0:  t  is  a  period  of  x}.   P  ?^  0  since  if 

T  is  a  period  of  x,  T  7^  0,  then  x  =  xT.   Thus 

x(-T)  =  xT(-T)  =  x(T-T)  =  X  so  -T  is  also  a  period  of  x. 

Since  either  T  or  -T  is  positive,  ?   ^  0.      Now  set  T  =  inf  P. 

Such  a  T  exists  since  P  is  bounded  below  by  0  and  hence 

has  an  inflmum.  Clearly  T  >_  0  so  suppose  T  =  0.   Then  there 

is  a  sequence  {t  }  in  P  with  t   ->-  0.   Since  x  =  xt   for 

n  n  n 

all  n,  X  must  be  a  critical  point  by  Theorem  2.19.   Hence 

T  >  0.   We  see  that  T  is  also  a  period  of  x  since  there  is 

a  sequence  {t  }  in  P  with  t   -»-  T  because  T  =  inf  P  and 
^        n  n 

X  =  xt   for  all  n.   Thus  by  the  continuity  axiom, 

X  =  xt   ->  xT.   Thus  T  is  a  period  of  x  and  by  definition 

it  is  the  smallest  positive  period  of  x.   Finally  let  t  e  R 

be  any  period  of  x.   If  t  7^  nT  for  any  integer  n,  then 

there  exists  an  integer  n  such  that  nT  <  t  <  (n  +  1)T. 

But  since  T  is  a  period  of  x,  nT  is  a  period  of  x  by 

lemma  2.15.   Thus  we  have  x  =  xt  =  x(nT)  which  gives 

xt(-nT)  =  x(nT)(-nT)  =  x(nT-nT)  =  x.   Then 

x  =  xt(-nT)  =  x(t-nT)  which  shows  that  t-nT  is  also  a  period 

of  X.   But  since  0  <  t-nT  <  T  we  have  a  contradiction 

since  T  v;as  the  smallest  positive  period  of  x.  // 

REMARK  2.20   If  x  e  X  is  periodic,  then  yCx)  and  y  (x) 
are  compact. 


l\Q 


PROOF  Let  X  be  periodic  with  positive  period  T.   Then 
y(x)  =  x[0,T]  and  y  (x)  =  x[0,T]  since  for  any  t  e  [0,T], 
xt  G  x[0,T]  and  for  any  t  with  t  >  T,  t  =  t^  +  kT  for  some 
integer  k  and  t,  e  [0,T].   Clearly  similar  results  hold  for 
t  <  0.   From  Remark  2.6,  x[0,T]  is  compact.   Thus  y(x)  and 
Y  (x)  are  compact.  // 


Prom  this  analysis  it  is  now  clear  that  l_f  any  trajectory 
in  a  dynamical  system  ever  crosses  itself,  then  the 
trajectory  must  be  periodic .   For  suppose  xt,  =  xt^  for 
t^  <  t^.      Then  x  =  x(tp  -  t  )  =  xT  where  T  =  t   -  t   >  0. 
But  this  shows  that  x  must  be  periodic. 
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III.   PROLONGATIONS  AND  LIMIT  SETS 


Referring  again  to  Example  2.2,  It  vias   noted  that  In  the 

damped  spring  case  the  solution  converges  to  the  point  (q) 

as  t  -»-  CO.   This  solution  converges  In  the  sense  that  there 

is  a  sequence  {t  }  with  t   £  R  for  all  n,  t   ->  »,  and 

n        n  '   n    ' 

xt   -V  (n)«   This  point  is  called  a  positive  limit  point, 
n    0 

Formally,  the  definition  is  stated: 

DEFINITION  3.1   Define  maps  A**",  A"  from  X  into  2^  by  setting 
for  each  x  e  X, 

1.  the  set  A  (x)  =  {y  e  X:  there  is  a  sequence 

{t  }  in  R  with  t   ->  +oo  and  xt   -^  y), 
n  n  n   '^  ' 

2.  the  set  A~(x)  =  {y  e  X:  there  is  a  sequence 

{t  }  in  R  with  t  -»■  -oo  and  xt  ^  y}. 
n  n  ^ 

For  any  x  e  X,  the  set  A  (x)  is  called  its  positive  limit 
set,  and  the  set  A~(x)  is  called  its  negative  limit  set . 

REMARK  3.2   If  x  e  X  is  periodic,  then  a'^(x)  =  A"(x)  =  y(x). 

PROOF  Let  X  G  X  be  periodic.   By  Remark  2.20,  y{x)    is 
compact.   Let  y  e  y(x).   Then  y  =  xt  for  some  t  e  R.   Since 
X  is  periodic,  y  =  x(t  +  T)  =  xT(t  +  T)  =  x(t  +  2T).   Clearly 
this  idea  can  be  extended  inductively  so  consider  the 
sequence  {t  +  nT}.   Since  T  >  0,  t  +  nT  -^  +«>  and  y  =  x(t  +  nT) 

for  all  n.   Thus  x(t  +  nT)  ->-  y  and  y  e  A  (x).   Hence 

+  + 

Y(x)  C  a  (x).   Now  let  y  e  A  (x).   Then  there  exists  a 
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sequence  {t^}  In  R  with  t   ->-  +°o  and  xt  -^   y.   For  all  n, 
xt^  e  y(x).   Hence  ixt^}  <Z  y{x) .      Since  this  sequence 
converges,  it  must  converge  to  a  point  in  y(x)    because 
Y(x)  is  compact.   Thus  y  e  y(x)  and  A  (x)  O  y{x).      Hence 
A  (x)  =  y(x).   The  proof  for  A~(x)  =  y(x)  is  entirely 
analogous.  // 

An  immediate  consequence  of  this  remark  is  that  i_f  a 
point  X  is_  a  critical  point ,  then  A  (x)  =  A~(x)  =  x.   From 
these  two  observations,  it  is  an  easy  task  to  determine  the 
positive  and  negative  limit  sets  of  Examples  2.1,  2.2,  and 
2.3.,  because  in  Example  2.1  we  have  critical  points, 
periodic  trajectories,  and  trajectories  which  clearly  do 
not  converge.   In  Example  2.2,  we  have  either  a  periodic 
trajectory  for  the  undamped  case  or  a  trajectory  whose 
positive  limit  set  is  the  origin  and  negative  limit  set  is 
empty.   Finally  in  Example  2 . 3  we  have  a  periodic  orbit. 
However,  a  point  may  not  be  periodic  and  yet  have  more  than 
one  point  in  its  positive  or  negative  limit  set.   The 
following  examples  help  to  illustrate  this  idea. 

EXAMPLE  3-2   Consider  the  differential  system  defined  in 

2 

R  by  the  following  equations  in  polar  coordinates: 


dt  =  ^^^  -  ^^ 


d0_  ^  -, 
dt 
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An  extensive  analysis  is  not  necessary  here  because  all  the 
facts  needed  can  be  derived  directly  from  the  equations. 
This  system  does  in  fact  define  a  dynamical  system  [Bhatia, 
Szego,  p.  20].   Since  the  change  in  6  is  constant,  only  the 

change  in  r  affects  the  system.   Clearly  for  r  =  0, 

dr 

-TTT  =  0  SO  the  origin  is  a  critical  point.   For  r  =  1, 

dr 

•^   =  0  so  r  is  again  constant  with  time  and  the  solution 

describes  a  periodic  trajectory.   Suppose  r  <  1.   Then 

dr 
dt 


is  positive  so  r  is  increasing  with  time.   If  r  >  1, 
dr 


then  -^-rr   is  negative  and  r  is  decreasing  with  time.   Such 
dt      ^  ^ 

trajectories  are  shown  in  the  diagram  below.   It  is  noted 

dr 
that  although  -^   is  always  positive  for  r  <  1,  r  will  never 


reach  r  =  1,  for  as  r 


dr 
dt 


0. 


FIGURE  3.2 

As  this  diagram  shows,  for  points  p  with  0  <  r  <  1,  A  (p) 
is  the  unit  circle  and  A~(p)  is  the  origin.  For  points  p 
with  r  >  1,  A"''(p)  is  the  unit  circle  and  A~(p)  =  0.   These 
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facts  become  clear  when  it  is  realized  that  a  sequence 

{t  }  can  be  found  with  t   ->  «  such  that  xt   approaches  any 
n  n  n 

given  point  on  the  unit  circle.   Since  the  unit  circle  is 

periodic  and  the  origin  is  a  rest  point,  their  positive 

and  negative  limit  sets  are  obvious.   The  direction  of 

the  trajectories  is  counter-clockwise  because  6  is  increasing 

with  time. 

EXAMPLE  3.3  Consider  the  differential  system  defined  in 

2 

R  by  the  differential  equations 

X  =  f(x,y) 
y  =  g(x,y) 

where  the  functions  f  and  g  are  defined  by 


0  if  X   >  1 


f(x,y) 


y(l  -  x^) 


(1  +  y^)(l  -  p(x)q(y)) 


if  Ixl  <  1, 


Here  the  functions  p(x)  and  q(x)  are  any  continuously 
differentiable  functions  satisfying  the  following  conditions: 
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0  <  p(x)  <  I  if  0  <  X  <  1, 

p(x)  =  0  if  X  <_  0, 

p(x)  =1  if  X  >  1, 

0  <  q(y)  <  i-  if  y  <  0, 

q(y)  =0  if  y  >  0. 

This  system  defines  a  dynamical  system  [Bhatia,  Szego  p.  21] 
and  it  is  analyzed  in  the  same  manner  as  the  preceding 
example.   For  x  =  y  =  0,  f(x,y)  =  0  and  g(x,y)  =  0  so  x  and 
y  are  unchanging  with  time.   Consequently,  the  point  (^)  is 
a  critical  point.   For  x  >_  1,  f(x,y)  =  0  and  g(x,y)  =  1. 
The  trajectories  of  such  points  are  straight  lines  parallel 
to  the  y-axis  with  direction  in  the  +y  direction.   A  similar 
analysis  shows  that  for  x  <_  -1 ,  the  trajectories  are  straight 
lines  parallel  to  the  y-axis  with  direction  in  the  -y 
direction.   Now  suppose  0  <  x  <  1.   Then  g(x,y)  >  0  so  y 
is  increasing  with  time.   If  y  <  0,  f(x,y)  >  0  so  x  is 
increasing.   If  y  >  0,  then  f(x,y)  <  0  so  x  is  decreasing. 
A  similar  analysis  holds  for  -1  <  x  <  0 .   Thus  the  solutions 
describe  spiralling  trajectories.   The  question  is  do  they 
spiral  outward  or  toward  the  origin  or  are  they  periodic 
trajectories.   Such  an  analysis  from  the  equations  is  very 
involved  but  the  trajectories  do  in  fact  spiral  outward 
[Bhatia,  Szego,  p.  22].   The  trajectories  are  sketched  in 
the  diagram  below. 

5^ 


IJ 


FIGURE  3.3 

From  the  diagram  one  can  see  that  for  any  point  except  the 
origin  in  the  strip  -1  <  x  <  1,  its  positive  limit  set  is 
the  two  straight  lines  x  =  ±1  and  its  negative  limit  set 

is  the  origin.   Since  the  origin  is  critical,  its  positive 

2 
limit  set  is  itself  and  for  any  other  point  in  R  ,  its 

positive  and  negative  limit  sets  are  empty. 

Proceeding  now  with  the  concept  of  limit  sets,  some  very 
important  theorems  are  useful. 


THEOREM  3.4   For  an^  x  e  X, 

1.  the  sets  A  (x)  and  A"(x)  are  closed  and  invariant » 

2.  the  closure  cl  (y  (x))  =  y    (x)  U  A  (x)  and  the 
closure  cl  (y~(x))  =  Y~(x)  VJ  A~(x). 

PROOF   1.   Consider  the  case  of  A  (x).   Let  {y  }  be  a 

'^n 


sequence  in  A  (x)  with  y 


n 


If  it  can  be  shown  that 
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+  + 

y  e  A  (x),  then  A  (x)  is  closed.   Proceeding  In  this  manner, 

it  is  noted  that  for  each  positive  integer  k,  there  is  a 

sequence  {t  }  in  R  with  t   -»■  +0°  and  xt  -^   y,  since  each 
n  n  n   ''k 

y,  e  A  (x).   We  may  assume  without  loss  of  generality  that 

k  k 

d(y,  ,xt  )  <  1/k  and  t   >  k  for  n  >  k  where  d  is  the  metric 
''k'   n  n  —         — 

for  the  space  X.   These  assumptions  may  be  made  since  the 

xt   get  arbitrarily  close  to  y,  and  hence  any  point  in  the 

sequence  at  distance  greater  than  1/k  from  y,  may  be  removed 

with  the  resultant  subsequence  converging  to  y,  .   The  same 

reasoning  applies  to  t   ^  k.   Consider  nov;  the  sequence 

{t  }in  R  with  t   =  t^.   The  significance  of  this  "diagonal" 
n  n    n         ^  ^ 

sequence  is  in  consideration  of  (xt  }  and  for  purposes  of 
clarity,  we  sketch  the  latter  sequence  below:- 


^c^[x<:t^i:-? 


Then  t   -»■+«>  and  the  claim  is  that  xt   -^  y 
n  n 

observe  that 


To  see  this. 


d(y,xt^)  <  d(y,yj^)  +  d(y^,xt^)  <  d(y,y^)  +  1/n. 


Since  1/n  and  d(y,y  )  tend  to  zero,  this  implies  that 
d(y,xt  )  -^  0.   Consequently,  xt   ->  y  and  y  e  A  (x).   Thus 


n 


n 


+  + 

A  (x)  is  closed.   To  see  that  A  (x)  is  invariant,  let 
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y  e  A  (x)  and  t  e  R  be  arbitrary.   There  Is  a  sequence 

{t^}  in  R  with  t^  ->  +00  and  xt  -*-  y.   Then  by  the  continuity 

axiom,  xt  (t)  -y   y  t .   Since  xt  (t)  =  x(t   +  t)  and 

(t^  +  t)  ->-  +00  this  implies  that  yt  z   a"*'(x).   Thus  a''"(x)  is 

invariant . 

To  prove  the  second  part,  it  is  recalled  that  y  (x)  =  xR  . 

By  the  definition  of  a"^(x)  it  follows  that 

cl  (y"*'(x))  3  y'^(x)  U  a'''(x).   Now  let  y  e  cl  (y"*"(x)).   Then 

there  is  a  sequence  {y  }  in  y  (x)  such  that  y   ^-  y .   Since 

this  sequence  is  in  y  (x)  it  follows  that  y   =  xt   for  some 
^  '  ''n     n 

t   e  R.   Either  the  sequence  {t  }  has  the  property  that 

t  -^  +<»,  in  which  case  y  cA  (x),  or  by  the  Bolzano-Weierstrass 

theorem  there  is  a  subsequence  {t   }  with  t   -»■  t  £  R  . 

"k        "k 

But  then  xt   ^  xt  which  is  in  v  (x),  and  since  also 
"k 

xt   -»■  y,  it  follows  that  y  =  xt  e  y  (x)*   Thus 
"k 

+        +       + 

cl  (y  (x))C  y  (x)  U  A  (x).   The  proofs  for  the  negative 

case  are  entirely  similar.  // 


The  properties  of  limit  sets  just  proved  will  be  found 
to  be  very  useful  in  further  work. 

THEOREM  3.5   For  x,y  e  X,  if  y  e  y'*"(x),  then  A"^(y)  =  A  (x) 

PROOF  Let  y  e  y  (x).   Then  there  is  a  t  e  R  such  that 

y  =  xt .   Now  let  u  e  A  (y).   Then  there  exists  a  sequence 

{t  }in  R  with  t  ->  +»  such  that  yt  ^  u.   Consider  the 
m  m  m 
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sequence  {t  +  t  }  in  R.   Clearly  t  +  t   ->  +«>.   Then 

in  m 

yt   ->-  u  Implies  yt^  =  xt  (t  )  =  x(t  +  t  )  -*-  u.   Thus 
m  m       m  m 

+  +       + 

u  e  A  (x)  and  hence  A  (y) d  A  (x). 

On  the  other  hand,  let  z  e  A  (x).   Then  there  is  a 

sequence  {t^}  in  R  with  t   ->  +00  and  xt   ->  z .   Consider  the 

ni  mm 

sequence    (t      -  t}    in  R.      Clearly    (t      -  t)    -»■  <»   since   t      -►  ». 

Then  xt      ->   z   implies   xt      =   x(t   +   t      -  t)    =   xt(t      -   t)    = 
n  '^  n  m  m 

y(t^  -  t)   ->   z.      Thus   z   e   A'^(y)    and   a"*'(x)  C  A'^(y). 

REMARK  3.6   Let  x,y  e  X  such  that  x  is  not  periodic.   Then 
the  setOCy  (y):  y  e  y  (x)}  =  0. 

+  + 

PROOF   Since  y  e  y  (^) >    there  is  a  t  e  R   such  that  y  =  xt . 

+  +  + 

Now  suppose  s  e  A{y  (y):  y  e  y  (x)}.   Then  s  e  y  (y)  which 

implies  that  s  =  yt^  for  some  t   e  R  .   But  since  y  =  xt , 

it  follows  that  s  =  xt(t-.)  =  x(t  +  t,  ) .   Let  tp  e  R  such 

that  tp  >  t  +  t,.   Then  for  y  =  xtp,  s  li   y''"(y).   Thus 

s  ^  f\{y    (y):  y  £  y  (x)}and  we  have  a  contradiction.   Hence 

n^y'^'Cy):  y  e  y"^(x)}=  0.     ^  .  // 

REMARK  3.7   For  any  x  e  X,  a''"(x)  =A{c1  (y"^(y)):  y  e  y'^(x)}  = 
/\{y    (xn):  n  is  an  integer}. 

PROOF  If  X  is  periodic,  the  proof  is  trivial  since 
cl  (y  (y))  =  A  (x)  for  all  y.   So  we  assume  x  is  not  periodic 
Let  {t  }  be  any  sequence  in  R  such  that  t  -^   +°°.   Without 
loss  of  generality  it  may  be  assumed  that  t   >  0  for  all  n 
since  any  points  in  the  sequence  for  which  this  does  not 
hold  may  be  removed  from  the  sequence  without  affecting 
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t^  -*■  +~.   Since  the  intersection  in  the  statement  is  taken 

over  all  possible  y's  in  y  (x),  set  y  =  xt   where  each 

n     n 

^n  ^  ^^n^-   ^^^  ^'^'^^n^  "  ^^n^'  ^  c  R"^)  =  {x(t^  +  t):  t  e  R^}. 
The  claim  is  that  H  {cl  (Y"^(y))}  =  ^  (cl  (y'^iy^^)))^ 

By  Theorem  3.^,  fl  (cl  (y^Cy))}  =  n{y'^(y)U  A"^(y)}  = 
f\{y'^(y)}V    f^ih'^iy)}    and 

n{cl  (y^(y^))}  =  n{y"*"(y^)  U  A^(y^)}  =n{y''(y^)}U  n{A-^(y^)} 
By  Remark  3.6,  /0{y'*"(y)}  =  J3  =  A^Y"^(y^)}  and  by  Theorem  3-5, 
A'^(y)  =  A"^(x)  =  A'^(y^).   Thus 
n  {cl  (y'^(y))}  =  nih^iy)}   =  A'^(x)  and 
n  {cl  (y'^(y^))}  =  n{A'^(y^)}  =  A'^(x).   Then,  as  claimed, 

n  {cl  (y  (y))}  =  A  {cl  (y  (y  ))}  and  as  shown  above, 

+  +        + 

n {cl  (y  (y)):  y  e  y  (x)}  =  A  (x).   To  prove  the  second 

part,  let  the  sequence  {t  }  e  (n).   Then  the  proof  is 
identical.  // 

In  the  following  proof  and  elsewhere  in  this  thesis, 
the  notation  S(M,£)  and  S[M,e]  will  be  used  for  the  sets 
{x:  d(x,M)  <  e},  and  {x:  d(x,M)  <_  e}  respectively  where 
M  C  X,  e  e  R  and  d  is  the  metric  for  the  space  X.   In  the 
case  where  M  =  {x},  the  notation  will  be  shortened  to 
S(x,e)  and  S[x,e]. 

THEOREM  3.8   Let  x  e  X.   If  cl  (y'*"(x))  is  compact,  then 
A  (x)  is_  a  non-empty  compact  and  connected  set , 

PROOF  Let  cl  (y  (x))  be  compact  and  let  {t  }  be  any 

sequence  in  R  such  that  t   ->■  +».   Now  consider  the  sequence 
^  n 
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{xt^}.   Since  xt^^  e  cl  (y"^(x))  for  all  n  and  cl  (y'*'(x)) 
is  compact,  this  sequence  has  a  convergent  subsequence 
{xtj^}which  converges  to  a  point  y  e  cl  (y^Cx)).   Thus  since 

^k  ^  "*"°°  ^^^  ^^k  "^  ^'  ^^  have   A  (x)  7^  0  since  y  e  a"^(x). 
Now  cl  (y  (x))  =  y  (x)  U  a"^(x).   By  Theorem  3.4,  a"''(x)  is 
a  closed  set  and  since  it  is  a  closed  subset  of  a  compact 
set,  A  (x)  is  itself  compact. 

The  proof  of  connectedness  is  by  contradiction.   Suppose 

+  + 

A  (x)  is  compact  but  not  connected.   Then  A  (x)  =  P  U  Q 

where  P,Q  are  non-empty,  closed  disjoint  sets.   Since 

cl  (y  (x))  is  compact,  it  is  locally  compact  and  hence  there 

exists  an  e  >  0  such  that  S[P,e]  and  S[Q,e]  are  compact  and 

disjoint.   To  see  that  S[P,e]  is  compact,  it  is  noted  that 

in  a  locally  compact  space,  every  point  has  at  least  one 

compact  neighborhood.   Taking  the  union  of  all  these  compact 

neighborhoods,  yields  a  compact  neighborhood  of  the  set  P. 

Then  there  is  an  e  >  0  such  that  S[P,e]  is  a  closed  subset 

of  this  neighborhood  and  hence  S[P,g]  is  compact.   The 

same  reasoning  applies  to  S[Q,e].   Since  a  metric  space 

is  normal,  there  exist  disjoint  open  sets  containing  P  and 

Q  respectively.   Hence  for  a  proper  choice  of  e,  the  sets 

S[P,e]  and  S[Q,e]  are  disjoint  and  compact. 

Now  let  y  e  P  and  z  e  Q.   Then  there  exist  sequences 

{t  }  and  {t  }  with  t  -*■+«>   and  t   ->■  +0°  such  that  xt^^  y 
n        n        n  n  n   '^ 

and  XT   ->  z.   Without  loss  of  generality  it  may  be  assumed 
that  xt   G  S(P,e),  XT   e  S(Q,e)  and  t^^  ~  ^n  ^  °  ^°^  ^^■'"  ^* 
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Since  the  trajectory  segments  x[t  ,t  ] ,  n  =  1,  2,  . . .  are 

compact  connected  sets,  they  clearly  Intersect  S(P,£)  and 

S(Q,e).   Thus  in  particular,  there  is  a  sequence  {T  } 

with  t^  <  T^  <  T^  such  that  xT   e  H(P,g)  where 

H(P,e)  =  {x:  d(x,P)  =  e}.   Since  H(P,e)  is  a  closed  subset 

of  S[P,e],  H(P,e)  is  compact.   Since  {xT  }  is  a  sequence 

n 

in  H(P,e),  it  has  a  convergent  subsequence  {xT  }  which 

^         m 

converges  to  a  point  y'  and  because  T  ->-  +«>,  it  follows 

m     ' 

that  y'  e  A  (x).   But  y'  t^  ?  0   Q   v;hich  contradicts  the 
assumption  that  A  (x)  =  P  U  Q.   Hence  A  (x)  is  connected.  // 

Since  only  the  local  compactness  of  cl  (y  (x))  was  used 
in  the  proof,  an  immediate  corollary  to  this  theorem  is 
that  i_f  the  space  X  is_  locally  compact ,  then  A  (x)  is_  connected 
whenever  it  is  compact . 

The  converse  of  the  previous  theorem  holds  if  the  space 
X  is  locally  compact.   However,  the  proof  relies  on  the  one- 
point  compactification  of  X.   The  fact  that  the  one-point 
compactification  is  an  extension  of  a  dynamical  system  to 
another  dynamical  system  is  proved  in  the  following  lemma. 

LEMMA   Given  a  dynamical  system  (X,R,7t)  on  a  locally  compact 
space  X,  let  X'  =  X  U  {w}  be  the  one-point  compactification 
of  X  and  define  tt  '  :  X'xR  ->  X'  by 


7r(x,t)    if  X  e  X,   t  e  R 
TT'(x,t)  = 

w     if  X  =  CO,   t  G  R 
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Then  (X',R,7t')  Is  a  dynamical  system  on  X' 

PROOF  First,  by  definition  of  tt', 

7t(x,0)  =  X   if   X  e  X 
TT  '  ( X  ,  0 )  = 

0)        if  X  =  0) 

Thus  7t'(x,0)  =  X  for  all  x  e  X'. 

Secondly,  if  x  e  X,  then  tt  '  (tt  '  (x,t^)',t2)  =  TT'(x,t,  +  tp) 

follows  trivially  so  suppose  x  =  w  .   Then 

TT' (tt' (to,t^)  ,t2)  =  TT'(a3,t2)  =  '^  =  TT'(w,t^  +  t  ^)  .      Thus 

tt' (tt' (x,t^),t2)  =  TT'(x,t-j_  +  t^)  for  all  x  e  X'.   It  remains 

to  show  that  tt'  is  continuous.   To  see  this,  we  will  utilize 

the  characterization  that  if  the  inverse  image  of  an  open 

set  is  open,  then  the  map  is  continuous. 

So  let  U  be  any  open  set  in  X' .   If  w  ^  U,  then  U  is 
open  in  X  and  thus  (tt')"^(U)  =  tt~-^(U).   Now  tt-^(U)  is  open 
because  tt  is  continuous  and  hence  (tt')~  (U)  is  open.   If 
to  e  U,  then  X-U  is  compact  and  since  w  ^  X'-U,  X'-U  =  X-U 
and  hence  (tt' )~^(X'-U)  =  (tt ')""'- (X-U)  =  tt"-^(X-U).   Now  since 
tt  is  continuous,  tt~  (X-U)  is  closed  and  hence  (tt')~  (x'-U) 
is  closed  so  (tt')~  (U)  is  open.   Thus  the  inverse  image  of 
every  open  set  is  open  and  hence  tt  '  is  continuous.        // 

THEOREM  3.9   If  the  space  X  is_  locally  compact ,  then  for 
any  x  e  X,  c1(y  (x))  is_  com.pact  whenever  A  (x)  is_  non-empty 
and  com.pact . 
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PROOF  Let  A  (x)  be  non-empty  and  compact.   Consider  the 

extension  of  the  dynamical  system  to  the  one-point  compactl- 

flcatlon  X'  =  X  U  {oj}.   Certainly,  for  any  x  e  X, 

A'^(x)  =  A  (x)  U  {a)}or   A'^(x)  =  A'''(x)  where  A'^(x)  Is  the 

positive  limit  set  of  x  in  X'.   However  since  A"'"(x)  is 

compact,  it  follows  that  A  (x)  =  A'*"(x)  holds  because 

A  (x)=cl^,(A  (x))  =  cl^(A  (x))  =  A  (x).   Also  we  have  that 

for  any  x  e  X,  y  (x)  =  y  (x)  since  7r'(x,t)  =  7T(x,t).   Thus 

we  have 

cl^.C-^/Cx))  =  y'^(x)  U  A^(x)  =  y'^(x)  U  A'^'Cx)  =  cl^Cy'^Cx)) 


and  since  cl„, (y  (x))  is  a  closed  subset  of  a  compact  space, 
it  is  compact.   Thus  cl(y  (x))  is  compact.  // 

THEOREM  3.10   If  cl(y'''(x))  is_  compact,  then  d(xt ,  a'*'(x)  )  -^  0 
as  t  ->-  +00. 

PROOF  Let  cl(y  (x))  be  compact  and  suppose  d(xt,A  (x))  /  0 

as  t  -»-  +00.   Then  there  exists  a  sequence  {t  }  in  R  such  that 

M        n 

t  -►  +00  and  d(xt  ^A  (x))  >  e  for  all  n  >  N  where  N  is  an 

integer  and  e  >  0.   But  since  cl  (y  (x))  is  compact,  the 

sequence  {xt  }  has  a  convergent  subsequence  (xt,  }  which 

converges  to  a  point  p  e  cl  (y  (x)).   Now  since  t   -*■  +00^ 

it  follows  that  t,  ->■  +°°  and  since  xt, -^  p,  we  have  that 

p  e  A  (x).   Also  since  xt,  ■+   p,  this  sequence  gets  arbitrarily 

close  to  p  and  hence  d(xt  ,p)  <  e  for  k  sufficiently  large. 
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Hence  we  have  d(xt^,  A  (x))  <  g  for  k  sufficiently  large 
but  this  is  a  contradiction  and  the  theorem  is  proved.    // 

The  study  of  limit  sets  can  be  very  useful  in  predicting 
.the  behavior  of  a  system.   If  the  system  describes  a  physical 
problem  and  the  trajectories  are  sketched,  it  is  a  relatively 
simple  task  to  discover  the  limit  sets  and  hence  predict 
the  future  behavior  of  the  system.   For  example,  consider 
again  the  pendulum  problem.  Example  2.1.   If  it  is  known 
that  the  pendulum  is  started  at  the  equilibrium  position  and 
given  an  initial  angular  velocity  of  w  =  2k,  then  one  would 
predict  that  the  pendulum  V70uld  approach  the  position  where 
9  =  IT  and  would  never  change  direction  to  swing  back.   This 
result  would  be  known  without  having  someone  watch  the 
pendulum  to  actually  see  what  happens.   It  is  true  that 
this  is  a  very  elementary  problem,  but  the  ideas  involved 
have  broader  applications  to  problems  of  much  greater 
difficulty. 

It  is  important  that  the  reader  have  a  firm  understanding 
and  familiarity  with  the  symbols  and  sets  thus  far  defined 
because  the  following  additional  sets  and  symbols  become 
quite  important  in  the  next  two  sections  of  this  thesis. 

DEFINITION  3.11  For  each  x  e  X,  define  the  following 
subsets  D  ,  D~,  J  ,  J~  of  X, 

1.   the  set  D  (x)  =  {y  e  X:  there  is  a  sequence  (x  } 

in  X  and  a  sequence  {t  }  in  R  such  that  x^  -»■  x  and  x  t   -^  y), 

n  n  n  n 

en 


2.  the  set  D  (x)  =  {y  e  X:   there  is  a  sequence  {x  } 

n 

in  X  and  a  sequence  {t  }  in  R~  such  that  x  ->  x  and  x  t   -»-  y } , 

n  n  n  n   '^  ' 

3.  the  set  J  (x)  =  {y  e  X:   there  is  a  sequence  {x  } 

n 

in  X  and  a  sequence  {t  }in  R  such  that  x  ->  x ,  t   ->  +°°, 

n  n    '   n     ' 

and  x  t  ->■  y), 
n  n   "^  * 

^.   the  set  J  (x)  =  (y  e  X:   there  is  a  sequence  {x  } 

^        n 

in  X  and  a  sequence  (t  }in  R"  such  that  x  ->  x,  t  ->  -co . 

n  n    '   n     ' 

and  X  t  -^  y}. 

For  any  x  e  X,  the  set  D  (x)  is  the  first  positive  prolongation 

of  X  and  the  set  D~(x)  is  the  first  negative  prolongation 

of  x. 

The  sets  J  (x)  and  J~(x)  are  called,  respectively,  the 

first  positive  and  the  first  negative  prolongational  limit 

set  of  X. 


These  sets  are  called  "first"  because  higher-order 
prolongations  can  be  defined.   However,  in  this  thesis, 
such  higher-order  prolongations  will  not  be  discussed. 
Consequently,  the  word  "first"  will  be  suppressed. 

From  the  definitions  it  is  clear  that  any  x  e  X, 
D"^(x)  D  y"^(x),  D"(x)  3>  y~(x),  J'^(x)  Z)  A"^(x),  and  J"(x)  Z>  A"(x) 
To  achieve  an  intuitive  feeling  for  these  sets  and  to  see 
that  the  inclusions  just  noted  may  be  proper,  we  present 
the  following  example. 

EXAMPLE  3-12  Consider  a  predator-prey  problem  where  it 
is  assumed  that  a  single  predator  population  feeds  on  a 
single  prey  population.   Let  x  denote  the  size  of  the  prey 
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population  and  let  y  denote  the  size  of  the  predator 
population.   Then  this  problem  can  be  expressed  by  the 
system 


^  =  ax  -  bxy,   x  >  0, 


^  =  cxy  -  dy,   y  >  0, 


where  a,b,c,  and  d  are  all  non-negative  constants. 

These  equations  support  our  intuition  since  we  would 
expect  that  any  change  in  the  size  of  the  predator  population 
depends  not  only  on  the  size  of  the  prey  population  but  also 
on  the  size  of  the  predator  population.   For  example  if  the 
prey  population  is  small  and  the  predator  population  is 
large,  we  would  expect  that  the  predator  population  would 
decrease.   This  expectation  is  confirmed  by  the  equations. 

Rather  than  solving  this  system  in  terms  of  t,  it  will 
be  analyzed  directly  from  the  differential  equations. 
Clearly  for  x=0  =  y,  -T^=0  =  -r^  and  the  point  (q)  is 
a  critical  point.   One  can  also  see  from  the  equations 
that  the  point  (  ,?^)    is  a  critical  point.   For  any  point 


p,   we  have  ||  >  0  and  ||-  = 
trajectory  lies  on  the  x-axis.   Physically,  this  is  the 


p  =  (^),  we  have  ^  >  0  and  §t-  =  0-   Thus  its  associated 
^    0  *         dt         dt 


case  when  there  is  no  predator  population  so  the  prey 
population  increases.   For  any  point  p  =  (  ),  dt"  ~  "^  ^^^ 
-r^  <  0  so  its  trajectory  lies  on  the  y-axis  with  direction 
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in  the  negative  y-direction.   This  is  the  case  when  there  is 
no  prey  population  so  the  size  of  the  predator  population 
decreases.   For  small  x  and  large  y,  ;qj  <  0  and  ;r^  <  0. 


dt 


dt 


Thus  both  X  and  y  are  decreasing  and  x  continues  to  do  so 

dx 
until  y  =  a/b.   At  this  point  ^  =  0  and  as  y  continues  to 

decrease,  x  begins  to  increase.   When  x  reaches  a  value  of 

d/c,  -rp  =  0  and  as  x  continues  to  increase,  y  begins  to 

increase.   The  physical  situation  described  by  this  case 

is  quite  clear.   The  trajectories  for  this  system  are 

sketched  below. . 

r 


a/b- 


FIGURE  3.12 

For  any  point  p  =  (  )  we  have  A  (p)  =  (q)  but 
j'*'(p)  =  {(^):x  >  0}.   The  fact  that  A"*'(p)  =  (q)  means  that 

if  there  is  no  prey  population,  the  predator  population  will 

+       x 
eventually  decrease  to  nothing.   The  fact  that  J  (p)  =  (q) 

means  that  even  though  the  prey  population  may  decrease 

considerably  so  that  the  predator  population  decreases  to 

nothing,  there  will  be  enough  of  the  prey  population  left 
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to  eventually  increase  its  size  to  any  amount.   We  also 

notice  from  this  diagram  that  for  any  point  p  =  (°), 

A'^(p)  =  (^)  but  J'^(p)  =  {(^):x  >  0}.   Also 

D'^(p)  =  Y^(p)  U  {(q)},  D"(p)  =  y"(p)  and  J"(p)  =  0.   For 

any  point  p  =  (^),  D'^(p)  =  y'^Cp),  j'^(p)  =  0,  D"(p)  =  Y~(p)U 

{(  )},  and  J"(p)  =  {(..)}.   For  any  other  point  p, 
<y  y 

d'*"(p)  =  y'^(p),  D"(p)  =  Y~(p),  J'^(p)  =  0  =  J~(p).   It  may 
be  argued  here  that  we  cannot  form  infinite  sequences  since 
"size"  is  not  a  continuous  function.   However,  in  an  actual 

experiment,  the  following  values  for  the  constants  a,b,c,d 

—  R 
were  observed:   a  ~  1,  b  ~  .01,  c  ~  3  x  10   ,  and  d  ~  0.^5. 

These  values  were  obtained  by  observing  paramecia  feeding 

dx 
on  yeast  cells.   The  value  of  y  for  which  ^  changes  sign 

is  then  at  y  -  100  and  the  value  of  x  for  which  -^   changes 

sign  is  at  x  -  1.5  x  10  .   These  values  are  large  enough 

to  assume  continuity.   Also  any  mathematical  model  idealizes 

the  physical  situation:   the  question  is  how  "close"  [Bailey, 

p.  2^5,  246]. 

Now  we  consider  some  basic  properties  and  important 

theorems  about  these  limit  sets. 


THEOREM  3.13   For  any  x  e  X, 

1.  the  set  D  (x)  is_  closed  and  positively  invariant , 

2.  the  set  J  (x)  ls_  closed  and  invariant , 

3.  the  equality  d'^(x)  =  y'^(x)U  J''"(x)  holds. 
Analogous  results  hold  for  D~(x)  and  J~(x). 
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PROOF  Statement  1  will  follow  from  statements  2  and  3  so 
statement  2  is  proved  first.   This  proof  parallels  the  proof 
of  Theorem  3.^.   Let  {y  }  be  a  sequence  in  j"'"(x)  with 


n 


y  -^  y.   For  each  integer  k,  there  are  sequences{x^}  in  X 
n  ^        n 

and{t   }in  r"^  with  xjj  ^  x,  t^  -^  +~,  and  x^t^  ^  y,  .   Without 
n  n    '   n     '      n  n   •'k 

loss  of  generality,  it  may  be  assumed  that  t^  >  k, 

^^\>^)   1   1A>  and  ^^^^^n'^k^  -  "^^^  ^°^  n  >  k.   Now  consider 


the  sequences  (x'^},  {t^}.   Clearly  x^  ->  x,  t^  ->■  +°°, 

n  '   n  ''   n    '   n     ' 


and 


X  t   ->  y.   This  is  true  since 


^^''n^n'y^  -  ^^^n^n'^n^  ^   d(y^,y)  <  1/n  +  d(y^,y).   Thus 

+         + 
y  e  J  (x)  and  J  (x)  is  closed. 

+  + 

To  see  that  J  (x)  is  invariant,  let  y  e  J  (x)  and  t  e  R. 

There  is  a  sequence  (x  )  in  X  and  a  sequence  (t  }  in  R 

n  ^        n 

such  that  t  -^   +«,  X  -^  X,  and  x  t^  -»■  y.   Now  consider  the 
n     '   n    '      n  n   '' 

sequence  (t   +  t}.   Clearly  t   +  t  ->  +«°,  and 

x  (t   +  t)  =  X  t  (t)  ->  yt.   Since  x  ->-  x,  it  follows  that 
n  n         n  n      "^  n    ' 

yt  e  J  (x).   Since,  t  e  R  v^ras  arbitrary,  it  follows  that  • 

J  (x)  is  invariant. 

Consider  statement  3.   Observe  that 

D  (x)  ^  Y  (x)Vj  J  (x)  always  holds  from  the  definition.   So 

let  y  e  D  (x).   Then  there  is  a  sequence  {x  }  in  X  and  a 

sequence  {t  }  in  R  with  x  ->'  x,  x  t   -*■  y.   It  may  be  assumed 
^        n  n    *   n  n   *^        '^ 

that  either  t   -»-teR  ort   -*■  +°°,  if  necessary  by  taking 
n  n     * 

subsequences.   In  the  first  case  x  t   ->  xt  by  the  continuity 
axiom.   Hence,  xt  =  y  e  y  (x)  as  t  e  R  .   In  the  second  case 
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y  e  J  (x)  by  definition.   Thus  y  e  y"''(x)  U  j"'"(x)  and 
D"^(x)0  y'*'(x)  U  J'*'(x).   Therefore,  d"*"(x)  =  y'*'(x)  U  J'^(x). 

To  see  statement  1,  recall  that  cl  (y''"(x))  =  y'^(x)  \J   a'*"(x) 
But  it  has  been  shown  that  a'^(x)  C  J'*'(x).   Hence 
D"^(x)  =  cl  (y"^(x))  V   J'^(x),  because  j'^(x)  73  a'^(x).   Since 
both  these  sets  are  closed  and  positively  invariant,  it 
follows  that  D  (x)  is  closed  and  positively  invariant.    // 

The  following  theorem  will  also  be  found  useful  in 
subsequent  sections. 

THEOREM  3.1^  Let  X  be  locally  compact.   Then  A'*"(x)  t^   0 
whenever  J  (x)  is_  non-empty  and  compact . 

PROOF   Suppose  a"*'(x)  =  0.   Then 

cl  (y  (x))  =  Y  (x)  U  A  (x)  =  Y  (x)  since  A  (x)  =  0.   Hence 

Y  (x)  is  closed.   It  is  claimed  that  y  (x)  is  disjoint  from 

J  (x)  for  if  Y  (x)  n  J  (x)  7^   0,  then  by  the  invariance  of 

+      +       +  + 

J  (x),  Y  (x) C  J  (x).   Since  J  (x)  is  compact,  for  any 

sequence  {t  }  in  R  such  that  t  ->■  +«>,  the  sequence  {xt  } 

must  have  a  convergent  subsequence.   It  follows  that  the 

point  to  which  that  subsequence  converges  is  a  point  in 

A'''(x)  and  hence  A'^(x)  7^  0.   Thus  y''"(x)  A  j''"(x)  =  0  as 

claimed.   Since  J  (x)  is  non-empty  and  compact, 

d(Y  (x),J  (x))  >  0.   Thus  since  X  is  locally  compact,  there 

is  a  6  >  0  such  that  S[J  (x),6]  is  compact  and  disjoint 

from  Y  (x).   Now  choose  any  y  e  J  (x).   There  is  a  sequence 

4. 
{x  }  in  X  and  a  sequence  {t  }  in  R  such  that  x^  -*■  x, 
n  ^        n  n 
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^n  "^  "*""'*  ^^^  ^n^n  "^  ^*   "^^  '^^^  ^^  assumed  that 


X  <!fS[J  (x),6],  x^t^  eS[J'^(x),6]  for  all  n.   Then  the  tra- 
jectory segments  x^[0,t  ]  intersect  H(j'''(x),6)  where 

n    n  V   s  /  ,  / 

+  + 

H(J  (x),6)  =  {z  e  X:d(J  (x),2)  =  6}.   Therefore  there  is 

a  sequence  (x^^},  0  <  t^  ^  t  ,  such  that  x  t   e  H(j''"(x)  ,6 )  . 

Since  H(J  (x),6)  is  compact,  it  may  be  assumed  that 

x^Tj^  -*■  z   e   H(J  (x),6).   By  taking  subsequences  if  necessary, 

it  may  be  assumed  that  either  x   -»-  t  e  r"^  or  t   -»■  +°°.   If 

n  n 

T^  ->  t  G  R  ,  then  by  the  continuity  axiom  x  t   -»■  xt  =  z 
n         J       ^  "^        n  n 

+  +         + 

so  z  e  Y  (x)  which  contradicts  y  (x)  ^i   S[J  (x),6]  =  0. 

If  T   ->  +~,  then  z  e  J  (x)  ,  but  this  contradicts 

z  e  H(j'*"(x),5)  as  j"''(x)  O  H(j"*'(x),6)  =  0.   Thus  the  original 

assumption  that  A  (x)  =  0  is  untenable  and  hence 

A"^(x)  7^  0.  // 

THEOREM  3.15   Let  X  be_  locally  compact.   Then  J"*'(x)  is_ 
non-empty  and  compact  if  and  only  if  D  ( x )  is_  compact . 

PROOF  Let  J  (x)  be  non-empty  and  compact.   Then  by  the 
preceding  theorem,  A  (x)  is  non-empty  and,  since  it  is  a 
closed  subset  of  J  (x),  it  is  compact.   But  then  cl  (y  (x)) 
is  compact  by  Theorem  3. 9.   Hence 

D"^(x)  =  y'^(x)  U  J'*'(x)  =  cl  (t'*'(x))  \J  J"^(x)  so  that  d"^(x) 
is  compact.   On  the  other  hand,  suppose  D  (x)  is  compact 

and  assume  J  (x)  =  0.   But  J  (x)  =  0  implies  y  (x)  is 

+  +  + 

compact.   Then  y  (x)  =  cl  (y  (x))  so  cl  (y  (x))  is  compact 

but  this  implies  A  (x)  7^  0  by  Theorem  3.8  and  the  assumption 
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+  + 

that  J  (x)  =  0  Is  contradicted.   Thus  J  (x)  5^  0  and  since 
,  J  (x)  is  a  closed  subset  of  d'''(x),  it  follows  that  j''"(x) 
is  compact.  // 

Before  proceeding  to  the  next  theorem,  the  set  d"*"(M) 

must  be  defined  where  M  ci  X.   Thus  d'^(M)  =  {y  e  X:  there  is 

a  sequence  {x^}  in  X  and  a  sequence  {t  }  in  R  such  that 
n  n 

X  ->•  X  e  M,  and  x  t   -^  y}. 
n        '      n  n   "^ 

THEOREM  3.16   For  an^L  compact  set  M  C  X,  the  set  d'*"(M)  is 
closed  and  positively  invariant . 

PROOF  Let  {y  }  be  any  sequence  in  D  (M)  such  that  y   -^  y . 

k  k 

Then  for  each  y,  there  exist  sequences{x  }  in  X  and  {t  } 
•^k  ^        n  n 

in  R  such  that  x  ^  ->-  x,  and  x  t   ->  y,  where  x,  e   M.   Without 

n    k      n  n   ''k        k 

k  k 
loss  of  generality  it  may  be  assumed  that  d(x  t  ,y,)    <   1/k. 

Now  consider  the  sequences  {x  }  in  M  and  {t  }  in  R  .   Then 

x  t  ->■  y   and  since  M  is  compact,  it  may  be  assumed  that 
n  n   "^n  t-        >  j 

x  -*■  X  e  M.   Thus  x  t   -^  y  since 
n  n  n   '^ 

d(x^t[},y)  <  d(x[]t^,y^)  +  d(y^,y)  <  1/n  +  d(y^,y).   Then 

by  definition,  y  e  D  (M)  and  hence  D  (M)  is  closed.   To 

see  that  D  (M)  is  positively  invariant  let  y  e  D  (M)  and 

let  t  e  R"^  be  arbitrary.   Now  y  e  D  (M)  implies  there  exist 

sequences  {x  }  in  X  and  {t  }  in  R  such  that  x  -^^  x  e  M 
^  n  n  n 

and  X  t   -^  y.   Now  consider  the  sequence  {t   +  t}  in  R  . 
n  n  ^        n 

Since  x  t  -*-  y,  it  follows  that  x  t  (t)  =  x  (t^  +  t )  ->  yt. 
nn  nn       nn 

Thus  yt  e  D'''(M).   Since  t  e  R  was  arbitrary  it  follows  that 
D  (M)  is  invariant.  // 
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We  conclude  this  section  with  an  important  theorem. 

THEOREM  3.17   Let  x,y  in  X.   Then  y  e  j''"(x)  i£  and  only 
if  X  e  J"(y). 

PROOF  Let  y  e  J  (x).   Then  there  is  a  sequence  {x  }  in  X 

and  a  sequence  {t  }  in  R  such  that  x  ->■  x,  t   ->■  +0°,  and 
^        n  n    '   n     ' 

X  t  ->■  y.   Set  X  t  =  y  ,  and  t  =  -t  .   Then  y  -^  y, 
n  n  n  n   *^n'      n     n        ''n   "  ' 

T  -*-  -00  and  yx   =xt(T)  =  x(t   -t)  =  x  ->x. 
n     '     '^nn    nnn     nn    n     n 

Consequently  x  e   J~(y).   The  converse  holds  in  the  same 
manner. 


// 
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IV.   RECURSIVE  AND  DISPERSIVE  CONCEPTS 

Using  the  sets  defined  in  the  last  section,  we  continue 
with  our  analysis  of  dynamical  systems  by  studying  recursive 
and  dispersive  concepts.   To  motivate  the  idea  of  recursive- 
ness, consider  again  Example  3.3.   Suppose  v;e  consider  a 
set  A  to  be  the  union  of  two  circles  centered  respectively 
at  the  points  ("r,)  and  (^)   with  radius  r  where  0  <  r  <  1. 
In  this  example,  all  trajectories  in  the  strip  -1  <  x  <  1 
spiral  outward  from  the  origin  with  the  lines  x  =  -1  and 
X  =  1  as  their  positive  limit  set.   Thus  for  any  point  in 
the  strip  -1  <  x  <  1,  its  trajectory  v;ill  eventually 
intersect  A.   The  trajectory  will  pass  out  of  A  again  but 
it  nevertheless  will  return  some  time  later  and  repeat  the 

process.   This  is  shown  in  the  following  diagram. 

7 


The  set  A  is  said  to  be  positively  recursive  with 
respect  to  the  trajectory.   We  now  state  the  formal  definition 
of  recursiveness. 
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DEFINITION  ^.1   A  set  A  C  X  is  said  to  be  positively 
recursive  with  respect  to  a  set  B  C  X  if  for  each  T  e  R 
there  is  a  t  >  T  and  an  x  e  B  such  that  xt  e  A.   Negative 
recursiveness  may  be  defined  by  using  the  inequality  t  <  T. 
A  set  A  is  self  positively  recursive  whenever  it  is  positively 
recursive  with  respect  to  itself. 

The  simplest  example  of  a  self  recursive  set  is  the 
set  (x)  where  x  is  a  periodic  point.   By  placing  restrictions 
on  the  system  by  demanding  that  the  entering  of  points  from 
B  into  the  set  A  happen  with  some  regularity,  other  concepts 
such  as  almost  periodicity  can  be  discussed  although  such 
is  not  our  purpose  here.   We  now  define  some  concepts  and 
establish  properties  in  connection  with  the  concept  of 
recursiveness. 

DEFINITION  4.2   A  point  x  e  X  is  said  to  be  positively 
Poisson  stable  if  every  neighborhood  of  x  is  positively 
recursive  with  respect  to  {x}. 

This  definition  means  that  the  trajectory  through 
the  point  x  eventually  intersects  every  neighborhood  of  x. 
The  following  diagram  illustrates  such  a  situation. 
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To  give  a  better  understanding  and  some  characteriza- 
tions of  positively  Poisson  stable  points,  the  follov;ing 
theorem  is  presented. 

THEOREM  4.3   Let  x  e  X.   Then  the  follov;ing  are  equivalent. 

1.  The  point  x  is_  positively  Poisson  stable . 

2.  Given  a  neighborhood  U  of_  x  and  a  T  >  0, 
we  have  xt  e  U  for  some  t  >  T, 

3.  The  point  x  l£  contained  in  its  positive  limit 
set  A  (x) , 

h.      The  set_  cl  (y''"(x))  =  A'''(x), 

5.  The  set  y(x)  C  A"^(x), 

6.  For  every  e  >  0  there  is_  a  t  >_  1  such  that 
xt  e  S(x,e) . 

PROOF  The  equivalence  of  statements  1-  and  2  follows  from 

the  definition  and  our  previous  discussion.   To  see  the 

equivalence  of  statements  2  and  3,  consider  a  sequence 

(T  }  in  r"*"  such  that  T  -^  +«.   To  each  n  associate  the 
n  n 

neighborhood   S(x  ,1/n)  of  x.   Then  for  each  t^  there 

exists  at   >  T  such  that  xt^  e  S(x,l/n).   Clearly 

t   -*■  +«  and  xt  -^   X.   Thus  x  e  A  (x).   Conversely,  let  U  be 
n  n 

a  neighborhood  of  x.   Since  x  e  A  (x)  there  is  a  sequence 

{t  }  in  R"^  with  t   -^  +«  such  that  xt^  ^  x.   Thus  xt  must 
n  n  n  a 

eventually  be  in  every  neighborhood  of  x  so  that,  in 

particular,  given  any  T  >  0  there  exists  a  t^  >  T  such  that 

xt   e  U.   Next,  we  establish  the  equivalence  of  statements 
n 
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3  and  k.      Let  x  e  A  (x).   Then  since  A"*"(x)  is  positively 

invariant,  we  have  y    (x)  C  A'*'(x).   Since 

cl  (y  (x))  =  y'''(x)  U  A"^(x),  it  follows  that  cl  (y"^(x))  =  A'^(x). 

Conversely,  let  cl  (y  (x))  =  A''"(x).   But  since  we  know  that 

cl  (y  (x))  =  y'^(x)  U  A'*'(x),  it  follows  that  y'^ (x)  Cl   A'*"(x) 

and  hence  that  x  e  A  (x).   To  see  the  equivalence  of 

statements  3  and  5,  let  x  e  A  (x).   Then,  since  A'*'(x)  is 

invariant,  xt  e  A  (x)  for  all  t  £  R.   Thus  y(x)  C  A'''(x). 

Conversely,  let  y(x)  Cl   A  (x).   Then  for  any  t  £  R, 

xt  G  A  (x).   In  particular,  for  t  =  0  we  have  x  e  A  (x). 

Finally,  we  prove  the  equivalence  of  statements  3  and  6. 

To  this  end  let  x  e  A  (x)  and  let  £  >  0.   Then  there  exists 

a  sequence  (t  }  in  R  with  t   ->■  «  such  that  xt   -^  x. 
n  n  n 

Without  loss  of  generality,  we  may  assume  t   ^1  for  all  n. 

Then  since  xt   -^  x,  this  sequence  must  eventually  be  in 

every  neighborhood  of  x.   Hence  there  exists  at   ^1 

such  that  xt   e  S(x,iE:  ).   Conversely,  choose  a  positive 

null  sequence  (£  )  such  that  £   -^0.   Then  for  each  n,  we 
^        n  •  n  ' 

have  at   >  1  such  that  xt   e  S(x,£  )  .   Since  £   -^0, 
n  —  n      *  n  n    ' 

we  have  that  xt   -^  x.   Consider  the  sequence  (t  }  in  R  . 

n  ^        n 

Either  t   ->  +«>  or  it  contains  a  convergent  subsequence 
n  o         ^ 

{t   }  such  that  t   -^  t  £  R"*".   If  t   -^  +°°,  x  £  A"'"(x)  by 
^k  "k 

definition.   If  t   -*■  t ,  then  xt   -»-  xt .   But  since  xt„  ->-  x, 

^k  ""k  ^ 

we  have  x  =  xt  and  therefore  x  is  periodic  with  period  t. 

It  follows  that  X  £  A  (x)  and  the  proof  is  complete.     // 
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From  this  theorem  it  should  be  clear  that  whenever  a 
point  X  is  positively  Poisson  stable  then  so  is  xt  for  every 
t  e  R.   This  is  true  because  if  x  is  positively  Poisson 
stable,  then  x  e  A  (x)  which  implies  xt  e  a'*'(x)  for  all 
t  £  R  by  the  invariance  ofA'^(x).   But  a'*'(x)  =  A'''(xt)  for 
all  t  e  R  by  Theorem  3-5.   Hence  xt  e  A  (xt)  and  xt  is 
positively  Poisson  stable. 

According  to  the  above  theorem,  a  point  x  is  positively 
Poisson  stable  if  x  e  A  (x).   Similarly,  a  point  x  is 
defined  to  be  negatively  Poisson  stable  if  x  e  A~(x).   It 
is  said  to  be  Poisson  stable  if  it  is  both  positively  and 
negatively  Poisson  stable.   If  a  point  x  is  Poisson  stable, 
then  both  its  motion  and  its  trajectory  are  said  to  be 
Poisson  stable. 

From  the  preceding  theorem  observe  that  if 
cl  (y  (x))  =  A  (x),  then  x  is  positively  Poisson  stable. 
Now  consider  the  situation  where  y  (x)  =  A  (x).   In  that 
case  the  point  x  is  positively  Poisson  stable  since 
cl  (y  (x))  =  Y  (x)  \J  A  (x)  =  A  (x).   We  will  now  establish 

that  Y  (x)  =  A  (x)  ±f_   and  only  if  x  is_  a  periodic  point . 

+       + 
For  if  X  is  a  periodic  point,  it  is  clear  that  y  (x)  =  A  (x). 

On  the  other  hand,  if  y  (x)  =  A  (x),  then  x  e  A  (x)  and 

+  + 

by  the  invariance  of  A  (x),  y(x)  =  A  (x).   Then 

xt'  G  Y  (x)  for  t'  <  0  and  therefore  there  exists  a  T  >_  0 
such  that  xt'  =  xT.  Hence  xt'(t)  =  xt(t')  =  xT(t). 


78 


Therefore  xt  =  x(t  +  T  -  t')  for  all  t  e  R,  showing  that 
X  is  periodic  with  period  T  -  t'  >  0. 

With  this  last  fact  it  is  a  natural  question  to  ask 
whether  there  are  Poisson  stable  points  which  are  not 
periodic.   The  following  example  shows  that  such  points  do 
exist, 

EXAMPLE  4.3   Consider  the  dynamical  system  defined  on  a 
torus  by  the  system 


=  f(x,y). 


dt 

^=  af(x,y). 


dt 


where  f(x,y)  is  periodic  in  both  x  and  y  with  period  1. 
That  is  f(x,y)  =  f(x  +  l,y  +  1)  =  f(x  +  l,y)  =  f(x,y  +  1). 
Further,  i-je  have  that  f(x,y)  >  0  for  x  and  y  not  both  0  mod  1 
and  f(0,0)  =  0.   The  constant  a  is  an  irrational  number. 
Clearly  the  point  p  =  (^)    is  a  critical  point  because 

at  this  point  -rp  =  0  =  -r^.   Solving  this  system  simultaneously 

dx 
we  have  -r—   =  1/a  which  implies  y  =  ax  +  c  where  c  is  some 

constant  of  integration.   Thus  the  trajectories  describe 

straight  lines  of  slope  a  in  the  plane.   V/e  may  assume 

a  >  0.   Nov/  no  trajectory  can  contain  p  except  y(p)  =  'fp^ 

since  p  is  a  rest  point  and  since  -rr-  and  -rt-  are  both  positive 

for  all  other  points. 
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We  note  now  that  the  function  y  =  ax  +  c  is  not  defined 
at  the  point  p  since  at  that  point  ^  =  q"'   ^°  ^^  claim 
that  there  is  only  one  trajectory  which  approaches  the 
point  p  as  t  ->  «.   The  following  diagram  will  help  us  to 
determine  this  trajectory.   For  purposes  of  this  diagram, 
we  assume  0  <  a  <  1. 


Consider  the  trajectory  y  =  ax  +  (1  -  a).   As  x  -^  1,  it 
is  clear  that  y  -^  1.   Hence  this  trajectory  approaches 
the  point  (?")  as  t  ->  «.   But  on  the  torus,  the  point 
p  =  (^)  =  (|).   Hence  this  trajectory  approaches  p  as  t  -^^  «>. 
Generalizing  this  approach  to  the  case  of  an  arbitrary 
positive  value  for  a,  the  trajectory  which  approaches  p  as 
t  -»■  «>  is  given  by  the  equation  y  =  ax  +  (1  -  (a  mod  1)). 
To  verify  this  fact,  let  x  -^  1.   Then  y  -^  a  +  (1  -  (a  mod  1) ) 
=  1  +  (a  -  (a  mod  D).   But  a  -  (a  mod  1)  is  just  some 
integer  n  and  hence  as  x  ^  1,  we  have  y  -^  n  +  1 .   Thus 
as  t  -»-  «,  this  trajectory  approaches  the  point  (   ^   )  but 
again,  in  the  torus,  this  point  and  the  point  p  are  identi- 
fied.  Thus  this  is  precisely  the  trajectory  which  approaches 
p  as  t  -»►  «.   Call  this  trajectory  y^'      Similarly  there  is 
only  one  trajectory  which  approaches  the  point  p  as 
t   -^  -  V     It  is  given  precisely  by  y  =  ax  +  0.   It  is  clear 
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that  as  X  ->  0,  we  have  y  -^  0  and  this  verifies  that  this 
trajectory  does  in  fact  approach  p  as  t  ^  -<».  Call  this 
trajectory  Y-,* 

Now  we  claim  that  no  trajectory  in  the  torus  is  periodic 
If  there  is  such  a  trajectory,  then  it  will  pass  through 
the  same  point  twice  and  without  loss  of  generality,  we 
may  assume  that  this  point  lies  on  the  y-axis.   We  may 


make  this  assumption  because  the  trajectory  must  intersect 

0 
c 


the  y-axis.   Assume  that  this  point  is  the  point  (^)  where 


0  <  c  <  1.   Then  if  the  trajectory  given  by  y  =  ax  +  c 
intersects  this  point  twice,  we  have  (^)  =  i ^^   ^  ^)mod  1 
for  some  n  >  0.   But  this  implies  that  c  =  (an  +  c)  mod  1. 
Since  0  <  c  <  1,  we  know  that  c  =  c  mod  1  =  (m  +  c)mod  1 
where  m  is  any  integer.   Hence  it  follows  that 
(an  +  c)mod  1  =  (m  +  c)mod  1.   This  implies  that  an  =  k  for 
some  integer  k.   Now  since  n  7^  0 ,  it  follows  that  a  =  k/n. 
But  this  implies  that  a  is  a  rational  number  which  is  a 
contradiction.   Thus  no  trajectory  can  pass  through  the  . 
same  point  twice  and  hence  no  point  except  p  is  periodic. 

We  also  note  here  that  the  trajectories  Y]_  and  y^  are 
not  the  same  for  if  they  were,  then  the  equations 
y  =  ax  +  (1  -  (a  mod  1))  and  y  =  ax  +  0  must  be  equal. 
But  this  implies  that  a  mod  1  =  1  and  thus  a  must  be  an 
integer  which  contradicts  the  fact  that  a  is  irrational. 

We  now  sketch  several  of  the  trajectories  in  the 
following  diagram. 


81 


To  complete  our  analysis  of  this  example,  we  see  that 


we 


A  (Y2)  =  "tP^  and  that  h~{y^)    =  (p).   For  the  point  p, 
have  A  (p)  =  A  (p)  =  {p}.   For  any  other  trajectory,  we 
have  A  (y)  =  A  (y)  =  the  torus  since  these  trajectories 
are  not  periodic.   Also,  A  (y-,  )  =  A~(Yp)  =  the  torus. 
Hence  we  have  that  all  the  points  on  Y-.  are  positively 
Poisson  stable  since  for  each  x  e  Y-i  >  x  £  A  (x).   Also 
points  on  Yp  are  negatively  Poisson  stable.   All  other 
points  are  Poisson  stable. 

If  we  change  this  example  so  that  f(x,y)  >  0  for  all 
x  and  y,  then  we  will  have  no  rest  points  and  no  periodic 
trajectories.   Hence  every  trajectory  is  dense  in  the  torus 
and  the  positive  and  negative  limit  sets  of  each  point  is 
the  torus. 


We  now  introduce  the  notion  of  a  non-wandering  point. 

DEFINITION  4.4   A  point  x  e  X  is  said  to  be  non-wandering 
if  every  neighborhood  U  of  x  is  self  positively  recursive. 
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It  is  clear  from  this  definition  that  if  a  point  is 
positively  Poisson  stable  then  it  is  non-wandering.   This 
idea  can  also 'be  seen  in  the  following  theorem  v;hich 
gives  some  characterizations  of  non-wandering  points. 

THEOREM  4.5   For  any  x  e  X,  the  following  are  equivalent. 

1.  The  point  x  is_  non-wandering 3 

2.  The  point  x  is_  an  element  of  its  positive 
prolongational  limit  set  J  (x) , 

3.  Every  neighborhood  of  x  l_s  self  negatively 
recursive, 

H.      The  point  x  is_  an  element  of  its  negative 
prolongational  limit  set  J  (x) . 

PROOF  Assume  x  is  non-wandering.   Consider  a  null  sequence 

{E},0<e,e  -^  0 ,    and  a  sequence{t^}  in  R  with  t^  -*-   +». 
n  '      n'   n     '         ^       n  n 

Since  each  S(x,e  )  is  self  positively  recursive,  we  have  an 

X  e  S(x,e  )  and  ax   >  t^  with  x^i^  e  S(x,e^).   Since 

e  -)■  0  we  have  x  -»■  x  and  x  t   ->■  x  and  since  t^  -v  +«>  we  . 
n  n  •         n  n  n 

conclude  that  x  z   J''"(x).   Thus  statement  2  holds.   Assume 

x  e  J'*'(x).   Then  there  is  a  sequence  {x^}  in  X  and  a 

sequence  {t  }  in  R  with  x  ->  x  and  t^  -^   +<=°  such  that 

X  t  -»■  X.   Now  for  any  neighborhood  U  of  x  and  T  >  0 
n  n 

there  is  an  N  such  that  t^  >  T,  x^  e  U,  and  x^t^  e  U  for 
n  _>  N.   Thus  U  is  self  positively  recursive  and,  consequently, 
X  is  non-wandering.   The  equivalence  of  statements  3  and  4 
is  proved  in  the  same  manner.   It  remains  to  show  that 
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statements  2  and  ^  are  equivalent.   This  follov/s  directly 
from  Theorem  3.17.  // 

We  shall  now  prove  some  theorems  which  establish  the 
connection  between  Poisson  stable  points  and  non-wandering 
points. 

THEOREM  k.6      Le^  x  e  X.   Every  y  e  A^(x)  is  non-wandering. 

PROOF  It  is  necessary  to  show  that  if  y  e  A"*'(x)  for  some 

X  e  X,  then  y  e  J  (y).   If  y  e  A  (x),  then  there  is  a 

sequence  {t  }  with  t   ->-  +00  and  xt   -*■  y.   Since  t   ->■  +«>, 
n        n  n  n 

we  may  assume,  if  necessary  by  taking  a  subsequence,  that 

t   ,  -  t   >_  n  for  each  n.   Then  setting 

T   =  t  ,T  -  t   and  xt   =  x  we  have  x  -*-y,x(t,,  -t)  = 
n    n+1    n       n    n  n   *^  *   n  n+1    n 

^^n+1  "^  ^'  ^^^  "^n  ^   ^n+1  ~  "^n  "^  "'"°°*   "^^^^  y  e  J  (y)  and 

y  is  non-wandering.  // 

THEOREM  4.7   Let  P  C  X  be  such  that  every  x  e  P  is_  either 
positively  or  negatively  Poisson  stable .   Then  every 
X  G  cl  P  is_  non -wandering. 


PROOF  Let  {x  }  in  P  and  x  ->■  x.   We  must  prove  that 
n  n  ^ 

+  + 

X  e  J  (x).   For  each  n  we  have  that  either  x   e  A  (x  )  or 

x   e  A^Cx  ).   Thus  by  taking  a  subsequence  if  necessary,  we 

may  a'ssume  that  x   e  A  (x  )  for  all  n  or  that  x_  e  A  (x  ) 
*'  n       n  n       n 

for  all  n.   In  the  first  case,  for  each  n  there  is  a 

t   >  n  with  d(x  ,x  t  )  <  1/n.   Then  clearly 
n  n  n  n 

d(x,x  t  )  <  d(x,x  )  +  d(x  ,x  t  )  <  d(x,x^)  +  1/n.   This 
'nn—    'n       nnn—      n 


8^1 


shows  that  x^t^  ->■  x  and  consequently  x  e  j'*'(x).  In  the 
second  case  similar  considerations  show  that  x  e  J~(x). 
Thus  by  Theorem  ^.5,  every  x  e  cl  P  is  non-wandering.    // 

A  partial  converse  of  this  last  theorem  will  now  be 
presented  without  proof.   The  proof  is  entirely  constructive 
but  not  particularly  instructive.   A  proof  may  be  found 
in  Bhatia  and  Szego. 

THEOREM  4.8   Let  X  be  complete  and  let  every  x  e  X  be 
non-wandering.   Then  the  set  of  Poisson  stable  points  P 
is  dense  in  X. 

To  see  how  several  of  these  concepts  are  related,  we 
present  the  following  summary.   It  is  clear  that  a  critical 
point  is  a  periodic  point.   Also,  all  periodic  points  are 
Poisson  stable  since  the  positive  and  negative  limit  sets 
of  a  periodic  point  are  the  periodic  trajectory.   Certainly 
the  periodic  point  is  contained  in  its  periodic  trajectory. 
Finally,  all  Poisson  stable  points  are  non-wandering  points 

points  for  if  X  is  a  Poisson  stable  point,  then 

+       + 
X  G  A  (x)C  J  (x).   These  ideas  are  summarized  in  the 

following  diagrajn. 
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We  now  define  a  different  type  of  stability.   This 
notion  will  become  important  in  discussing  "dispersive" 
concepts. 

DEFINITION  ^.9   For  any  x  e  X,  the  motion  tt   is  said  to 
X 

be  positively  Lagrange  stable  if  cl  (y  (x))  is  compact. 
Further,  if  cl  (t'Cx))  is  compact,  the  motion  tt   is  called 
negatively  Lagrange  stable.   It  is  said  to  be  Lagrange  stable 
if  cl  (y(x))  is  compact. 

If  X  =  R  ,  then  the  above  statements  are  equivalent  to 
the  condition  of  the  sets  y    (x),  y~(x),  y(x)  being  bounded, 
respectively.   The  following  statements  about  this  notion 
are  easily  proved  using  Theorems  3.8  and  3.9  and  Remark  3.10. 

1.  If  X  is  locally  compact,  then  a  motion  tt   is 
positively  Lagrange  stable  if  and  only  if 

A  (x)  is  a  non-empty  compact  set. 

2.  If  a  motion  tt   is  positively  Lagrange  stable, 
then  A  (x)  is  compact  and  connected. 

3.  If  a  motion  tt   is  positively  Lagrange  stable, 
then  d(xt,A"*'(x))  -^  0  as  t  ->  +<». 

We  are  now  ready  to  discuss  dynamical  systems  which  are 
noted  by  their  lack  of  the  properties  of  recursiveness. 
One  can  easily  draw  a  parallel  between  these  dispersive 
concepts  and  the  recursive  concepts  just  defined.   First, 
we  need  some  definitions. 
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DEFINITION  4.10   Let  x  e  X. 

1.  The  motion  tt^  is  said  to  be  positively  Lapranpie 

II  X  ^--  —  —  —   I-  T  -  — w^   -      -  W  — 

unstable  v/henever  cl  (y  (x))  is  not  compact. 

2.  The  motion  is  said  to  be  negatively  Lagrange 
unstable  whenever  c1(y~(x))  is  not  compact. 

3.  The  motion  is  said  to  be  Lagrange  unstable 

if  it  is  both  positively  and  negatively  Lagrange 
unstable. 
k.      The  point  x  is  called  positively  Poisson  unstable 
whenever  x  ^   A  (x) . 

5.  The  point  x  is  called  negatively  Poisson  unstable 
whenever  x  ^   A~(x). 

6.  The  point  x  is  called  Poisson  unstable  whenever 
it  is  both  positively  and  negatively  Poisson 
unstable. 

7.  The  point  x  is  called  wandering  whenever 
X  ^  J    (x) . 

This  is  a  rather  lengthy  definition  but  the  definitions 
are  what  we  would  expect.   We  are  now  ready  to  define 
dispersive  concepts  on  the  whole  space  X. 

DEFINITION  4.11   The  dynamical  system  (X,R,tt)  is  said  to  be 

1.  Lagrange  unstable  if  for  each  x  e  X,  the  motion 
TT   is  Lagrange  unstable, 

2.  Poisson  unstable  if  each  x  e  X  is  Poisson  unstable, 

3.  completely  unstable  if  every  x  e  X  is  wandering. 
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^.   dispersive  if  for  every  pair  of  points  x,y  e  X, 
not  necessarily  distinct,  there  exist  neighbor- 
hoods U  of  X  and  U   of  y  such  that  U   is  not 
positively  recursive  with  respect  to  U  and  U 
is  not  positively  recursive  with  respect  to  U  . 

One  may  wonder  why  such  concepts  as  positively  and 
negatively  wandering  points  are  not  defined.   This  is 
unnecessary  since  x  e  J  (x)  if  and  only  if  x  e  J"(x).   One 
can  see  from  the  definition  that  a  completely  unstable 
system  must  be  Poisson  unstable  since  if  x  ?f  j''"(x),  then 
x  £f  A  (x).   The  following  remark  proves  the  relationship 
between  all  the  sets  defined. 


REMARK  ^.12   In  definition  4.11,  each  statement  implies 
the  preceding  statement. 

PROOF   Suppose  (X,R,'it)  is  dispersive.   Then  for  every  pair 

of  points  x,y  e  X,  there  exist  neighborhoods  U  of  x  and 

U  of  y  such  that  U  is  not  positively  recursive  with  respect 

to  U  .   Set  X  =  y.   Then  there  exists  a  neighborhood  of  x 

which  is  not  self  positively  recursive.   Hence  x  is  wandering 

so  X  ^  J  (x).   Since  this  is  true  for  all  x  e  X,  we  have 

that  X  is  completely  unstable.   Next  suppose  X  is  completely 

unstable.   Then  x  f^  J    (x)  for  all  x  e  X.   Then  also 

x  ^  J~(x)  for  all  X  e  X.   Hence  x  jzf  A  (x)  and  x  f^   A~(x) 

since  A  (x)  C  J  (x)  and  A~(x)  C  J^Cx).   Thus  X  is  Poisson 

unstable.   Finally  assume  X  is  Poisson  unstable.   Then 
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X  ^  A  (x)  and  x  i   A~(x)  for  all  x  e  X.   Suppose  cl  (y'^(x)) 
is  compact.   Then  by  Theorem  3.8,  A'^(x)  is  non-empty  and 
corapact,  hence  complete.   Then  by  Theorem  ^1.6  every  point 
in  A  (x)  is  non-wandering  and  by  Theorem  4.8,  the  set  P  of 
Poisson  stable  points  is  dense  in  A''"(x).   Thus  P  is  non- 
empty.  But  if  q  is  a  Poisson  stable  point  in  A'^(x),  then 
q  e  A  (q)  and  the  space  is  not  Poisson  unstable.   Hence 
we  have  a  contradiction  so  that  cl  (y'^(x))  is  not  compact. 
A  similar  argument  shows  that  cl  (y~(x))  is  not  compact. 
Hence  X  is  Lagrange  unstable.  // 

The  follov;ing  diagram  provides  a  summary  of  the  ideas 
in  the  preceding  remarks. 


Thus  we  know  that  a  dispersive  system  is  completely 
unstable,  a  completely  unstable  system  is  Poisson  unstable 
and  a  Poisson  unstable  system  is  Lagrange  unstable.   One 
may  wonder  if  the  reverse  implications  hold.   In  general 
they  do  not  as  the  following  examples  illustrate. 

To  show  that  a  Lagrange  unstable  system  may  not  be 
Poisson  unstable  we  want  x  e  A  (x)  but  cl  (y  (x))  not 
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compact.   Since  x  e  A  (x),  we  know  that  y"''(x)  C  A'*"(x)  and 
hence  cl  (y  (x))  =  A  (x).   Thus  we  must  have  an  example 
where  A  (x)  Is  not  compact  and  yet  x  e  A''"(x).   To  this  end, 
consider  Example  4.3  of  the  dynamical  system  defined  on  the 
torus.   We  delete  the  point  p  from  the  torus  so  that  the 
torus  is  no  longer  compact.   Then  for  any  x  in  the  torus 

such  that  X  ^  Y^,  we  have  A  (x)  =  the  torus  and  hence 

+         +  + 

X  e  A  (x)  but  A  (x)  is  not  compact.   Thus  cl  (y  (x)) 

is  not  compact.   Clearly  for  x  e  Y-,  ,  A~(x)  =  0  and  hence 

cl  (y"(x))  is  not  compact.   Similar  results  hold  for  x  e  Yp 

and  cl  (y  (x))  is  not  compact.   Thus  for  all  x  in  the  torus, 

cl  (y  (x))  is  not  compact  and  cl  (y  (x))  is  not  compact. 

Hence  this  system  is  Lagrange  unstable  but  it  is  not  Poisson 
unstable  because  for  any  x  ^   y^i    '^   ^   ^    (x). 

To  see  that  a  Poisson  unstable  system  may  not  be 
completely  unstable,  consider  the  following  phase  portrait 
in  the  x,y-plane.  ^ 


The  unit  circle  contains  a  rest  point  p  and  a  trajectory  y 
such  that  for  each  point  q  e  y  we  have  A  (q)  =  A~(q)  =  {p} 
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All  trajectories  in  the  interior  of  the  unit  circle  have 
the  same  positive  and  negative  limit  sets  as  y.   All  tra- 
jectories in  the  exterior  of  the  unit  circle  spiral  to 
the  unit  circle  as  t  ->•  +«>,  so  that  for  each  point  q  in 
the  exterior  of  the  unit  circle  v;e  have  A  (q)  =  {p}  \J  y, 
and  A  (q)  =  0.   If  we  consider  the  dynamical  system  obtained 
from  this  one  by  deleting  the  rest  point  p  then  the  resul- 
tant system  is  Poisson  unstable.   However  it  is  not 
completely  unstable  because  for  each  q  e  y,  J  (q)  =  y 
so  q  e  J  (q) . 

Finally  to  see  that  a  completely  unstable  system  may 
not  be  dispersive  consider  the  following  dynamical  system 
defined  by  the  differential  equations 

X  =  sin  y 

2 
y  =  cos   y 

This   system   contains   trajectories   y,     given  by 


Y^  =^(x,y):    y   =   kTT   +   tt/2},        k  =   0,    ±1,    ±2, 


These  are  lines  parallel  to  the  x-axis.   Between  any  two 

consecutive  y,  's  the  trajectories  are  given  by 
'k 

Y  =  {(x,y):  X  +  c  =  sec  y},  where  c  is  a  constant  depending 
on  the  trajectory.   The  phase  portrait  between  the  lines 
y  =  -7t/2  and  y  =  tt/2  is  shown  in  the  diagram  below. 
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It  can  be  seen  that  this  system  is  completely  unstable 


since  for 


any  p  e  y_^j  J  (p)  =  Y   and  for  all  other  points 


p  not  belonging  to  any  Yi^j  J  (p)  =  0.   Thus  p  i   J(p)  for 
any  p  and  hence  every  point  is  wandering.   However  if 
p  e  Y_-|  and  q  e  y  j  then  every  neighborhood  of  q  is  recur- 
sive with  respect  to  any  neighborhood  of  p  and  the  system 
is  not  dispersive. 

The  following  result  shows  that  to  prove  Lagrange 
instability,  it  is  sufficient  to  prove  that  the  space 
is  positively  Lagrange  unstable. 

THEOREM  4.13  A  dynamical  system  (X,R,7t)  is_  Lagrange  unstable 
if  and  only  if.  it_  is_  positively  Lagrange  unstable. 

PROOF  The  necessity  is  trivial.   On  the  other  hand,  let 
the  space  be  positively  Lagrange  unstable.   Then  for  each 
X  e  X,  cl  (y'*"(x))  is  not  compact.   Suppose  cl  (y~(x))  is 
compact.   Then  A~(x)  is  non-empty  and  compact.   Since  A  (x) 
is  invariant,  yt  e  A~(x)  for  all  t  e  R  and  y  e  A~(x).   Then 
since  A~(x)  is  closed,  cl  (Y'*'(y))C  A~(x)  and  hence 
cl  (Y"''(y))  is  compact.   This  is  a  contradiction.   Hence 
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cl  (y  (x))  is  not  compact  for  each  x  e  X.   Thus  the  space 

is  also  negatively  Lagrange  unstable,  yielding  the  result.  // 

For  a  dynamical  system  defined  in  the  x,y-plane,  the 
following  theorem  provides  an  important  fact. 

THEOREM  4.14   For  dynamical  systems  defined  by  a  system  of 
differential  equations  in  the  euclidean  plane,  the  concepts 
of  Lagrange  instability  and  complete  instability  are 
equivalent . 

PROOF  Let  this  system  be  given  by  the  differential 
equations 

X  =  f (x,y) , 
y  =  g(x,y) , 

where  f(x,y)  and  g(x,y)  are  continuous  functions. 

Clearly  if  this  system  is  completely  unstable,  then  . 
it  is  Lagrange  unstable.   So  assume  the  system  is  Lagrange 
unstable  and  that  there  is  a  point  p  which  is  non-wandering. 
Then  p  e  J"*'(p).   The  point  p  cannot  be  a  critical  point  for 
then  cl  (y''"(p))  would  be  compact  contrary  to  hypothesis. 
Likewise,  there  exists  a  neighborhood  V  of  non-critical 

points  about  the  point  p. 

P-, 
Let  p  =  (   ).   Since  p  is  not  a  critical  point,  then 
P2 

either  f(p^,P2)  7^  0  or  g(p-,^,P2)  ?^  0-   Assume  fCp^,?^)  i-    0. 
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Since  f(x,y)  is  continuous,  there  is  a  neighborhood  W  about 
p  where  f(x,y)  ^   0  for  all  points  (^)  e  V/.   Now  consider 
the  function  H:R^  ->  R  defined  by 

H(x,y)  =  tan"-''[||^?^j].   Clearly  H(x,y)  is  continuous  in 

W.   Since  H(x,y)  =  tan"   [^]  =  tan"  {-r-)  y   we  may  consider 

X  ^^ 

the  function  H  to  be  a  "slope"  function  since  it  provides 

us  with  the  slope  of  any  trajectory  at  a  given  point. 

Now  by  the  continuity  of  H(x,y)  there  is  a  neighborhood 

U'  such  that  for  all  points  (  )  in  U' , 

|H(x,y)  =  H(p^,Pp)|  <  e,  where  e  >  0  is  chosen  such  that 

U'  C  W.   The  relationship  between  the  sets  V,  W,  and  U' 

is  shown  in  the  following  diagram. 


Now  since  |H(x,y)  -  HCp^jP^)!  <  e  for  all  (^)    in  U', 
we  see  that  the  neighborhood  U'  contains  only  points  whose 
trajectories  have  slopes  that  differ  from  the  slope  of 
Y(p)  at  p  by  less  than  e.   Now  let  6  >  0  be  so  chosen 
such  that  6  <  e  and  6  <  1/10.   Now  let 

U  =  {(^)  e  U':  |H(x,y)  -  H(p^,P2)|  <  <5 } .   The  number  1/10 
was  chosen  so  that  we  are  assured  that  we  can  construct 
the  following  "transversal"  through  all  the  trajectories 
lying  within  U. 


9^ 


Let  L  be  a  line  segment  through  the  point  p  making 
an  angle  of  90°  with  the  motion  through  p.   Such  a  line 
can  be  constructed  since  we  are  working  in  the  plane.   Now 
every  trajectory  in  U  must  cross  this  line  L  by  our  above 
construction  of  U. 


We  claim  that  each  trajectory  crossing  L  in  U  must  do 
so  in  the  same  direction  as  y(p)«   To  see  this,  let  r  e  U 
and  suppose  y(^)  crosses  L  in  the  opposite  direction  of 
y(p).   Let  p'  be  the  point  where  y(^)  nieets  the  line  L. 
Now  let  U-L  =  U,  U  Up.   It  is  clear  that  U   and  Up  are 
disjoint  sets. 


Since  y(i")  and  y(p)  cross  L  in  different  directions, 
we  may  assume  that  for  some  small  time  t'  >  0,  p't'  e  U^ 
and  pt'  e  Up.   Now  let  L'  be  the  segment  of  L  which  joins 
the  points  p  and  p'  and  consider  Tr(L',t').   Since  L'  is  a 
connected  set,  'iT(L',t')  is  a  connected  set.   Also  part 


95 


of  L't'  v/111  be  in  U^  and  part  of  L't'  will  be  in  U 
since  p't'  e  U  and  pt '  e   \]    . 


Now  L't'  is  a  connected  set  so  there  is  at  least  one 
point  q  such  that  q  e  L'H  L't'.   Thus  it  follows  that 
q  =  qt '  since  it  would  be  impossible  that  any  other 
trajectory  passing  through  L'  at  time  t  =  0  would  cross 
the  point  q  in  time  t  =  t'.   For  this  to  happen,  we  would 
have  the  following  situation. 


Clearly  this  is  an  impossibility  in  U  and  hence  we 
have  q  =  qt ' .   But  this  implies  that  q  is  a  periodic  point 
and  hence  cl  (y  (q))  is  compact  and  we  have  a  contradiction. 
Hence  all  trajectories  cross  L  in  the  same  direction. 

Now  p  G  J  (p)  so  there  exist  sequences  {x  }  in  X 

and  {t  }  in  R  such  that  x  -^  p  and  t   ->  +00  with  x  t   -*  p. 
n  n  n  n  n 

If  any  of  the  x  's  are  periodic,  then  cl  (y  (x  ))  is  compact 
and  we  have  a  contradiction.   Hence  no  x  is  periodic  and, 
without  loss  of  generality,  we  may  assume  each  x  lies  on 
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L.   Now  for  n  sufficiently  large,  we  may  assume  that 
^n^n  ^  ^  since  x^t^  ->  p.   Then  it  follows  that  there 
exists  a  T  e  R  such  that  x^t^(T)  lies  on  L.   Call  this 
point  z. 

It  is  clear  that  z  lies  between  the  points  x  and  p. 

n     ^* 

for  if  not  we  have  the  situation  shown  below  and  it  is 
clear  that  in  this  case,  x^t^  -»■  p  is  an  impossibility. 


Thus  we  must  have  the  following  situation. 


Now,  let  L^  be  the  segment  of  L  v/hich  joins  the  points 

X  and  z.   Since  z  =  x  (t   +  T),  we  consider  the  set 
n  n  n     ' 

{x  [0,t   +  T]  U  L-,}  =  A.   This  is  a  Jordan  curve  and  hence 
n   *  n         1 

the  interior  of  A  together  v;ith  A  is  a  compact  subset  of 

2 

R  .   Let  B  denote  this  compact  set.   Now  p  e  B  so  we 

consider  y  (p).   We  claim  that  pt  e  B  for  all  t  e  R  .   If 

this  is  not  the  case,  then  y  (p)  must  cross  the  set  A 

+  + 

at  some  point  since  y  (p)  is  a  connected  set.   But  y  (p) 
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cannot  cross  A  anywhere  along  x  [0,t   +  T]  because  different 


trajectories  cannot  intersect.   Then  the  only  alternative 

1 


is  that  Y  (p)  niust  cross  L,  . 


But  this  too  is  impossible  since  to  do  so  would  mean 
crossing  L  in  the  opposite  direction  from  y(x^).   Hence 
pt  e  B  for  all  t  e  R"*".   Hence  y    (p)  C  B  and  since  B  is 
compact,  we  have  cl  (y'''(p))<=  B.   Therefore  cl  (y  (p)) 
is  compact.   Thus  we  have  our  final  contradiction  and  the 
assumption  that  p  e  j'''(p)  is  untenable.   Thus  the  system 
is  completely  unstable.  // 

(The  above  proof  was  communicated  to  me  through  a  private 
conversation  with  Dr.  William  Chewning. ) 

We  now  continue  with  several  characterizations  of 

dispersive  systems. 

< 

THEOREM  4.15   A  dynamical  system  (X,R,tt)  is_  dispersive  if 

and  only  if  J'^(x)  =  0  for  each  x  e  X. 


PROOF  Let  the  system  be  dispersive.   Assume  there  is  a 
point  X  e  X  such  that  J^(x)  5^  '  0  .   Then  if  y  e  j"*'(x) 
there  are  sequences  {x^}  and  {t^}  such  that  x^  ->  x. 
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t^  -y   +CO  and  x^t^  ->  y.   This  shows  that  for  any  neighborhoods 

U  of  X  and  U   of  y,  U  t  H  U,  ¥  0   as   the  element  x  t 
^  y        X  n    y  n  n 

is  contained  In  this  Intersection.   But  this  contradicts 

the  assumption  that  the  system  Is  dispersive  because  U   Is 

then  positively  recursive  with  respect  to  U  .   Hence 

J  (x)  =  0  for  all  X  e  X.   Conversely,  let  J'*"(x)  =  0   for 

all  X  e  X.   We  claim  that  there  exist  neighborhoods  U  of 

^  X 

X  and  U  of  y  and  a  T  >_  0  such  that  U  t  A  U  =0  for  all 
t  >_  T.   For  if  not  there  exist  sequences  (x  }  {y  }  {t  } 

such  that  x^  "^  ^»  ^n  ^  ^n^n'  ^n  "^  ^'  ^^^  ^n  "^  "'"°°'  ^°  ^^^^ 

y  e  J  (x).   But  this  is  a  contradiction.  Similarly  U 

is  not  positively  recursive  with  respect  to  U  .   Hence  the 
system  is  dispersive.  // 

THEOREM  4.16   For  a  given  (X,R,7t)  the  following  are  equivalent. 

1.  The  space  (X,R,7t)  1s_  dispersive . 

2 .  For  any  two  points  x , y  in_  X  there  are  neighbor- 
hoods U  of  X  and  U  of  y  and  a  constant  T  >  0 

X  y 

such  that  U  H  U  t  =  0  for  all  t.  It  I  >  T. 
X    y ill- 
s'  For  any  two  distinct  points  x,yinX,  y^fJ  (x) 

and  X  ^  J    (y ) . 

PROOF  Let  (X,R,7t)  be  dispersive.   Then  for  any  two  points 

x,y  in  X,  there  exist  neighborhoods  U  and  U   such  that  U 

X      y  X 

is  not  positively  recursive  with  respect  to  U  and  U  is 

not  positively  recursive  with  respect  to  U  .   Hence  for  some 

T  >  0,  U  n  U  t  =  0  for  all  t  >  T  and  U^t  H  U„  =  0  for  all 
'   X '  '  y  X     y 

t  >  T  which  implies  U^  fl  U,  (-t)  =  0.   Thus   U^  H  U  t  =  0  for 

X    y  A     y 

all  t,  |t|  >  T. 
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Assume  statement  2  and  suppose  x  e  J  (y).   Then  there 

exist  sequences  (y  }  and  {t  }  such  that  y^  "*"  y  and  t^  -♦■  +°° 

and  y  t   -»-  x.   For  n  sufficiently  large  v;e  may  assume  that 
n  n 

V   e  U  and  y  t   e  U   and  t   >  T.   Then  we  have 
•''n    y     '^n  n    x      n 

U  n  U  t   7^  0  which  is  a  contradiction.   The  proof  that 
X    y  n 

y  ^  J  (x)  is  similar. 

Finally  assume  that  for  any  tv/o  distinct  points  x,y  in 
X,  y  ^  J'^(x)  and  x  i   J"^(y).   Suppose  for  all  neighborhoods 
U  and  U  and  any  T  >  0  we  have  U^  0  U  t  ^   0  for  some 
t  >  T.   Since  this  is  true  for  all  neighborhoods,  it  will 
be  true  for  a  nested  sequence  of  neighborhoods  about  the 
points  X  and  y.   Hence  consider  for  each  n,  the  neighborhoods 
S(x,l/n)  and  S(y,l/n).   Now  set  T  =  T^.   Then 
S(x,l)  Pi  S(y,l)t^  f   0  for  some  t^  >  T.^  +  1.   Set  "Y^   =  t^. 
Then  there  exists  a  t^  >  T2  +  2  such  that 

S(x,l/2)  r\   S(y,l/2)t2  i-   0-      In  this  manner  we  have  constructed 
sequences{x^},  (y^),  (t^} ,  with  x^  --  x,  t^--^  ~>"  Vn  "  ^n' 
and  y   ->  y.   We  note  that  the  process  in  constructing  the 
sequence  {t  }  cannot  terminate  for  if  it  did,  then  we  have 
found  a  T,  >  0  such  that  S(x,l/k)  Pi  S(y,l/k)t  =  0  for  all 
t  >  T,  contrary  to  our  hypothesis.   Thus  we  have  x^^t^  -^  y 
so  y  e  J'*"(x)  but  this  is  a  contradiction  and  hence  U^ 
cannot  be  positively  recursive  with  respect  to  U^.   A 
similar  argument  holds  for  U^  not  positively  recursive 
with  respect  to  U^.   Hence  (X,R,tt)  is  dispersive.        // 
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We  now  introduce  the  fl-nal  concept  in  this  section. 
It  is  that  of  parallelizable  dynamical  systems.   We  are 
headed  toward  one  particular  theorem  and  consequently  many 
lemmas  and  theorems  will  be  presented  v;ithout  proof  since 
their  proofs  are  constructive  in  nature  or  are  not  particu- 
larly informative.   However,  the  proofs  may  be  found  in 
Bhatia  and  Szego  pages  48  through  53.   Formally  we  begin 
with  the  definition  of  a  parallelizable  dynamical  system. 

DEFINITION  4.17   A  dynamical  system  (X,R,tt)  is  called 
parallelizable  if  there  exists  a  set  S  C  X  and  a  homeo- 
morphism  h:  X  ->  SxR  such  that  SR  =  X  and  h(jct)  =  (x,t) 
for  every  x  e  S  and  t  e  R. 

To  proceed  in  the  study  of  parallelizable  dynamical 
systems  we  need  to  develop  a  theory  of  sections. 

DEFINITION  4.18   A  set  S  C  X  is  called  a  section  of  (X,R,tt) 
if  for  each  x  e  X  there  is  a  unique  time  t(x)  such  that 
xt(x)  c   S. 

The  set  S  in  Definition  4.17  is  in  fact  a  section  of 
(XjRjTt)  for  if  z  is  any  element  in  X,  then  z  =  xt  for  some 
X  e  S  and  t  e  R.   If  this  was  not  the  case,  then  SR  7^  X. 
Then  the  value  (-t)  is  precisely  t(z)  given  in  Definition  4.18 
It  is  clear  that  zt(z)  =  z(-t)  =  x  e  S.   To  see  that  (-t) 
is  unique,  assume  there  exists  (-f)  e  R  such  that 
z(_t')  =  y  e  S.   Then  xt  =  y t ' .   Now  using  the  homeomorphism 


101 


h,  we  have  h(xt)  =  h(yt')  which  Implies  (x,t)  =  (y,t'). 
Thus  we  have  x  =  y  and  t  =  t'.   Thus  t(z)  Is  unique. 

Then  to  summarize,  it  is  clear  that  if  a  dynamical 
system  is  parallelizable  it  must  have  a  subset  S  which 
contains  a  point  of  every  trajectory  so  that  SR  =  X.   More- 
over, the  set  S  can  contain  only  one  point  of  every  trajec- 
tory for  otherwise  t(x)  would  not  be  unique.   To  see  that 
this  last  statement  is  true,  suppose  there  is  a  trajectory 
which  intersects  the  section  twice.   Then  there  are  points 
x,y  in  this  trajectory  such  that  x,y  e  S  and  xt  =  y  for 
some  t  >  0.   Then  t(x)  =  0  since  xo  =  x  e  S  but  also 
xt  =  y  e  S  and  hence  t  =  t(x).   Since  t  >  0,  this  contradicts 
the  uniqueness  of  t(x).  This  implies  in  particular  that 
the  dynamical  system  cannot  have  any  critical  points  or 
periodic  trajectories.   To  see  how  a  parallelizable  system 
compares  with  the  other  dispersive  concepts  defined  in  this 
section,  consider  the  following  example. 

2 
EXAMPLE  ^.19   Consider  a  dynamical  system  defined  in  R  by 

the  differential  equations 


if  =  f(x,y). 


^  =  0 
dt   ^' 


where  f(x,y)  is  continuous,  and  moreover  f(x,y)  =  0  whenever 
the  point  (^)  is  of  the  form  (j^^)   ^'ith  n  a  positive  integer 
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For  simplicity  we  assume  that  f(x,y)  >  0  for  all  other 
points.   This  system  is  sketched  in  the  diagram  below. 


->- 


w 


->- 


w 


fr*r 


Ji^l 


r 


Now  we  form  the  system  in  which  we  are  interested  by 
deleting  the  sets 


I   =  {(^):  X  <  n,  y  =  1/n},   n  =  1,  2,  3,  •• 
n     y     - 


from  the  plane  R  .   This  system  is  sketched  below. 

Y 


in 


(IV 


-^ 


a\< 


i\^-^ 


It  is  clear  that  this  system  is  dispersive  since  the 

only  way  two  neighborhoods  U^  of  x  and  U^  of  y  could  be 

positively  recursive  with  respect  to  each  other  would  be 

if  the  points  X  and  y  lie  on  the  same  trajectory.   Then 

clearly  we  could  find  a  T  >  0  for  which  U^  0  U^t  =  0  or 

U  n  U  t  =  0  for  all  |tl  >  T.   Hence  this  system  is  disper- 
y    X 
sive  but  it  is  not  parallelizable  since  no  section  S  can 

exist  for  which  t(x)  would  be  continuous.   This  fact  can 

be  seen  when  it  is  realized  that  the  only  section  possible 
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is  one  that  looks  like  y  =  1/x  +  c  for  any  constant  c 
together  with  some  section  which  intersects  the  trajectories 
below  the  x-axis.   If  this  "lower"  section  intersects  the 
trajectory  on  the  x-axis,  then  t(x)  could  not  be  continuous. 
If  the  lower  section  looks  like  y  =  -1/x  +  c,  then  t(x) 
would  be  continuous  but  no  section  will  intersect  the  tra- 
jectory on  the  X-axis.   Hence  the  system  is  not  parallelizable 

We  can  get  an  example  of  a  parallelizable  dynamical 
system  from  the  previous  example  by  considering  the 
differential  equations 


dx 
dt 


=  1, 


^  =  0. 
dt 


This  system  is  sketched  below. 


(S) 


-^ 


Clearly  any  line  parallel  to  the  y-axis  is  a  section  of 
the  system. 

These  examples  also  give  an  intuitive  feeling  for  the 
notion  of  a  parallelizable  dynamical  system  in  that  all 
trajectories  are  in  some  sense  parallel. 
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The  function  t(x)  will  be  basic  in  our  development. 
In  general  the  function  t(x)  need  not  be  continuous  however 
for  a  dynamical  system  to  be  parallelizable  there  must 
exist  a  section  where  t(x)  is  continuous  on  X.   This  fact 
is  proved  in  the  following  theorem. 

THEOREM  4.20   A  dynamical  system  (X,R,7t)  is_  parallelizable 
if  and  only  if  it  has  a  section  S  with  t(x)  continuous  on  X. 

PROOF   Suppose  (XjRjir)  has  a  section  S  with  t(x)  continuous 
on  X.   Then  clearly  SR  =  X.   We  now  define  the  homeomorphism 
h:  X  -^  SxR  by  h(x)  =  (xx  (x)  ,-t  (x) ) .   To  see  that  h  is  one- 
to-one  let  h(x)  =  h(y)  for  x  and  y  in  X.   Then  xt(x)  =  yiCy) 
and  -t(x)  =  -xCy).   Hence  xt(x)  =  yT(x)  e  S  which  implies 
that  X  =  y  for  otherwise  T(y)  is  not  unique.   Hence  h  is 
one-to-one.   The  function  h  is  continuous  by  the  continuity 
of  t(x)  and  the  continuity  axiom.   The  inverse 
h""^:  SxR  ^-  X  is  given  by  h"''"(x,t)  =  xt  and  is  clearly 
one-to-one  and  continuous.   Thus  h  is  a  homeomorphism  so 
the  definition  of  a  parallelizable  dynamical  system  is 

satisfied. 

Now  assume  (X,R,tt)  is  parallelizable.   Then  the  set  S 
in  the  definition  is  a  section  of  X.   Since  for  any  x  e  X, 
x  =  yt  for  some  y  £  S  and  t  e  R,  we  set  t(x)  =  -t .   Then 
xt(x)  =  x(-t)  =  y  G  S.   The  continuity  of  t(x)  follows 
from  the  continuity  of  h  in  the  definition.  // 
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To  continue  this  development  we  need  to  define  special 
types  of  sections.   This  Is  done  In  the  following  definition 

DEFINITION  4.21   An  open  set  U  In  X  will  be  called  a  tube 
If  there  exists  a  t  >  0  and  a  subset  S  C  U  such  that 
SI  C  U  and  for  each  x  e  U  there  Is  a  unique  t(x), 

T 

|t(x)|  <  T,  such  that  xt(x)  g  S.   Here  I^  =  (-t,t).   U 
Is  also  called  a  x-tube  with  section  S,  and  S  a  (t  -  U)- 
sectlon  of  the  tube  U.   If  I^  =  R  then  U  Is  an  --tube  and 
S  an  (~  -  U)-sectlon. 

DEFINITION  4.22   Given  an  open  ~-tube  U  with  a  section  S 
and  t(x)  continuous  on  U,  and  given  sets  N  and  K  such  that 
N  C  K  C  S  where  N  Is  open  in  S  and  K  is  compact,  we  shall 
call  KR  the  compactly  based  tube  over  K.   Then  Indeed  t(x) 
restricted  to  KR  is  continuous  on  KR. 

Proceeding  toward  our  final  result,  the  next  theorems 
and  lemma  are  presented  without  proof. 

THEOREM   If  X  is  locally  compact  and  separable,  and  if 
every  x  £  X  is.  a  wandering  point ,  then  there  exists  a 
countable  covering  {K^R}  of  X  by  compactly  based  tubes 

K  R  each  with  t^(x)  continuous  on  K^R. 
n  n 

LEMMA  A  compactly  based  -tube  U  with  section  K  of  a 
dispersive  dynamical  system  (X,R,7t)  is  closed  in  X. 

We  also  need  the  following  important  topological 
theorem. 
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TIETZE  EXTENSION  THEOREM  Let  X  be_  a  normal  topological 
space ,  let  A  be_  a  closed  subset  of  X  and  let  f  be  a 
continuous  function  on_  A  to_  the  closed  interval  [a,b]. 
Then  f  has  a  continuous  extension  g  v/hich  carries  X  into 
[a,b]. 

The  final  lemma  needed  to  prove  our  result  is  quite 
instructive  and  is  completely  presented  belovj. 

LEMMA  Let  U, ,  U2  b£  two  compactly  based  tubes  of  a  disper- 
sive dynamical  system  with  sections  K^,  K2  and  continuous 
functions  t^(x)  and  t^Cx)  respectively.   If  U.^  0  ^^^^ 
then  U  =  U  VJ  Up  is  a  compactly  based  tube  with  a  section 
K  -D  K   and  a  continuous  function  t(x).   Moreover ,  if  the 
time  distance  between  K^  and  K2  along  tra^^ectories  in 
U  n  U2  is  less  than  x  (>  0),  the  time  distance  between  K 
and  K2  along  tra.jectories  in  U  is  also  less  than  t  . 

PROOF  In  this  proof,  the  following  diagram  will  be  useful 


(;.ni/j 


By  the  previous  lemma,  U^  and  U^  are  closed.   They  are 
invariant  because  each  tube  contains  a  section  K.^  and  ^^ 
respectively  and  K^R  C  U^  and  K2RC  U2.   Hence  U^  H  U^  is 
invariant  and  closed.   Further,  K^  H  U^  is  compact  and 
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non-empty,  as  shown  in  the  diagram.   Set  Sp  =  K^  H  U 
and  S   =  K  n  ^2'      ^^^   trajectory  in  U^ 0  U^  intersects 
S,  in  exactly  one  point  and  intersects  Sp  in  exactly  one 
point.   Thus  for  any  x  e  U^  O  U^ ,  t^(x)  =t^(x)    +   T^(xT2(y')). 
This  is  true  because  xt-,(x)  =  xt2(x)  (t-|^(xt2(x)))=  xd^  (x)  + 
T^(xTp(x)))  and  there  are  no  rest  points  or  periodic 
trajectories  in  a  dispersive  system.   The  function  t^  is 
continuous  on  Sp  v/hich  is  compact  and  now  we  apply  the 
Tietze  extension  theorem.   This  is  possible  since  every 
metric  space  is  normal.   Thus  the  continuous  function  t^ 
is  extended  to  a  continuous  function  t  defined  on  Kp  where 
t(x)  =  T^(x)  for  X  e  Sp.   If  t^(x)  g  (-t,t)  for  x  e   Sp , 
we  have  t(x)  e  (-t,t)  for  x  e  Kp. 

Notice  now  that  {xt(x):  x  e  Sp}  =  S^  and  t(x)  being 
continuous,  {xt(x):  xg   Kp}  is  compact  as  Kp  is  compact. 
Set  K  =  K-j_U  {xt(x):  x  e  Kp}  and  define  t*(x)  on 
KR  =  K-jR  VJ  KpR  as  follows: 

Tt^Cx)  for  x  G  K^R 

T*(x)  =  < 

(xpCx)  +  T(xTp(x))   for  X  G  KpR 

T*(x)  is  continuous  on  KR  and  we  need  only  verify  that  if 
X  G  U^n  Up,  then  T^(x)  =  Tp(x)  +  xCxip  (x))  which  has 
already  been  proved  since  t(x)  =  x-^Cx)  for  x  g  Sp.   Thus 
the  lemma  is  proved. 
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Finally,  we  prove  the  theorem  that  is  the  main  focus 
of  this  section  on  parallelizable  dynamical  systems. 

THEOREM  ^.23   A  dynamical  system  (X,R,7t)  on  a  locally 
compact  separable  metric  space  X  1_^  parallelizable  if  and 
only  if  it  is  dispersive. 

PROOF  Let  the  dynamical  system  (X,R,it)  be  dispersive.   By 

Theorem  ^.20  it  is  sufficient  to  prove  that  X  has  a  section 

S  with  t(x)  continuous  on  X.   By  a  previous  theorem,  there 

is  a  countable  covering  {U  }  of  X  by  compactly  based  tubes 

U  with  sections  K  and  continuous  functions  t  (x).   We 
n  n  n 

replace  this  covering  by  a  like  covering  {U  }  of  compactly 

1 
based  tubes  which  we  construct  as  follows.   Set  K^  =  K  , 

and  U,  =  U-"".   Beginning  with  U"'-  and  U^  we  use  the  previous 

1  2 

lemma  to  enlarge  K  to  a  compact  set  K  ,  thus  obtaining 

2     1  2 

the  compactly  based  tube  U  =  U  U  U2  with  t  (x)  continuous 

on  U^.   This  leaves  K,  unaltered.   Having  found  U  in  the 

same  manner,  we  take  it  together  with  U^^-^  and  construct 

U""^^  with  k'''*"^Z>k",  and  t^'^^(x)  continuous  on  u"""*"^.   Now 

set  S  =  U  k",  then  X  =  SR  and  the  function  t(x)  defined  by 

t(x)  =  T^(x)  for  X  e  u"  is  continuous  on  X,  with  the  property 

that  xt(x)  e  S.   Moreover  t(x)  is  unique  for  each  x  e  X. 

Thus  X  has  a  section  S  with  continuous  t(x)  defined  on  X. 

The  system  (X,R,tt)  is  thus  parallelizable. 

Conversely,  suppose  the  system  (X,R,it)  is  parallelizable. 

Then  there  exists  a  section  S  with  t(x)  continuous  on  X. 
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Suppose  J    (x)    ^  0   for  some  x  e  X.   Let  y  e  J''"(x). 

Then  since  y  e  X,  there  exists  y  e  S  and  t   e  R  such  that 

yt"  =  y.   Now  since  J  (x)  is  invariant,  y(-t  )  e  J'*'(x) 

which  implies  that  y  e  J  (x).  -Then  there  exist  sequences 

{x^}  and  {t^}  with  x^  -f  x,  t^  ->-  +<»,  and  x  t   -»■  y .   Now 

for  each  x  ,  there  exists  x  e  S  and  F  £  R  such  that 
n'  n         n 

X  F  =  X  .   Similarly  there  exists  x  e  S  and  F  e  R  such 
n  n    n  x 

that  xt   =  X.   Then  since  x  t  -^  y .  it  follows  that 
X  n  n   "^  * 


X  t  (t  ) 
n  n  n 


y 


Consider  for  each  fixed  k,  k  =  1,  2,  ...,  the  sequence 

{x  F  (t,  )}.   Since  x  t  ->■  xt  ,  we  have  by  the  continuity 
nnk  nn     x  ^ 

axiom  that  x  F  (t,  )  ->  xF  (t,  ) .   Hence  we  may  assume  by 
n  n   K      X   K 

taking  subsequences  if  necessary  that 

d(x  t  (t,  ),xF  (t,  ))  <  1/k  for  k  =  1,  2,  ...   Then  it  follows 
n  n  k  '   X  k   —  '   ' 

that  xF  (t  )  ■>  y  since 
X  n    "^ 


d(y,x  t  (t  ))  +  1/n.  Now  since  xt^(t^)  ->  y,  let 

T   =  t   +  t  .   Since  t   ->-  +°°,  it  follows  that  t^-^  +' 

n    X    n  n      '  '-^ 

Hence  we  have  that  y  e  A  (x). 


Now  for  each  n  set  y^  =  xt^.   Then  y^  "^  y  •   Since  t^ 
we  may  assume  that  t   >  n  for  each  n.   Let  U  be  an 


-^       +00 
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arbitrarily  small  neighborhood  of  y.   Then  since  y  -»■  y. 


n 


^n  ^  ^y  ^^^  ^-^^   n  >  N  where  N  is  a  large  enough  integer. 


We  know  that  t(x)  is  continuous  and  since  T(y)  =  0,  it 

follows  that  T(y^)  is  close  to  0  for  all  n  >  N.   Hence  we 

may  assume  that  (i^(y  )(  <  1/n  for  each  n  >  N.   But  this  is 

impossible  since  T(y  )  is  unique  and  we  have  already 

determined  that  T(y^)  =_t  <  -n  for  all  n.   Thus  we  have 

n     n 

a  contradiction  so  it  is  impossible  that  y  e  J  (x)  and 
hence  J  (x)  =  0  for  all  x  e  X.   Then  by  Theorem  4.15  we 
have  that  (X,R,it)  is  dispersive.  // 


The  relationships  between  the  various  dispersive  concepts 

is  now  conveniently  summarized  below. 


Ill 


V.   STABILITY  THEORY 

In  this  section  we  will  study  stability  theory  first 
from  the  abstract  point  of  vlev/  that  we  have  established 
and  then  from  an  applied  point  of  viev/  that  uses  systems 
of  differential  equations  and  studies  the  eigenvalues  of 
such  a  system.  To  begin  we  need  some  basic  definitions. 
Throughout  this  section >  it  will  be  assumed  that  the  space 
X  is_  locally  compact  and  the  set  M  is_  a  non-empty  compact 
subset  of  X  unless  explicitly  stated  otherwise . 

DEFINITION  5.1   With  a  given  M  C  X  we  have 

1.  the  set  A  (M)  =  {x  e  X:  A"^(x)n  M  i-    0}, 

w 

2.  the  set  A(M)  =  {x  e  X:  A'*"(x)  ¥■   0   and  a'^(x)  C  M} . 

The  set  A  (M)  is  called  the  region  of  weak  attraction  and 
w  — '^ 

the  set  A(M)  is  called  the  region  of  attraction  of  the  set 
M.   Moreover,  any  point  in  A^(M)  or  in  A(M)  is  said  to  be 
weakly  attracted  or  attracted  to  M  respectively. 

Prom  this  definition  it  is  clear  that  a  point  x  is 

weakly  attracted  to  a  set  M  if  and  only  if  there  is  a 

sequence  {t  }  in  R  with  t   ->  +~  and  d(xt^,M)  ->■  0.   Also  a 
n  n  n 

point  X  is  attracted  to  a  set  M  if  and  only  if  d(xt,M)  ->  0 
as  t  ->■  +«. 

The  above  attractor  sets  are  important  in  the  development 
of  stability  theory  and  the  following  theorem  provides 
needed  information  about  them.   The  following  lemma  is  useful 
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LEMMA  5.2   Let_  X  be  an^  metric  space  and  let  x  e  X.   Then 
for  every  t  e  R, 

1.  t2ie  set  h'^(x)    =  A"*'(x)t  =  A'*"(xt), 

2.  the  se^  j"*'(x)  =  j'^(x)t  =  j'*'(xt). 

PROOF  The  first  equality  In  statements  1  and  2  follows 

from  the  invarlance  of  the  sets  a'*'(x)  and  j"''(x)  respectively. 

To  see  the  second  equalities,  consider  the  case 

A  (x)t  =  A  (xt).   Let  z  e  A  (x)t.   Then  there  is  a  y  e  A"'"(x) 

with  z  =  yt  and  a  sequence  {t  }  ,  t  -»■  +«,  with  xt  ->   y. 

Then  by  the  continuity  axiom,  xt  (t)  -^^  yt .   But 

xt^(t)  =  xt(t  )  and  since  t  ->■  +«>,  we  must  have  that 

+  +        + 

yt  e  A  (xt).   Thus  A  (x)t  C  A  (xt).   The  argument  is  clearly 

+        + 
reversible  so  we  also  have  that  A  (xt)C  A  (x)t.   Kence 

statement  1  is  proved.   The  proof  of  the  second  inequality 

in  statement  2  is  entirely  analogous.  // 

THEOREM  5. 3   For  any  given  M,  A  (M)  "D  A(M)  and  the  sets 


w 


A^(M)  and  A(M)  are  invariant . 


PROOF  Let  x  e  A(M).   Then  a'''(x)  t^  0  and  A''"(x)C  M.   Then 

clearly  A  (x)  H  M  7^  0.   Hence  x  e  A^(M).   To  see  the 

invariance  of  the  sets  A  (M)  and  A(M) ,  consider  the  case 

w 

for  A  (M).   Let  x  e  A  (M)  and  let  t  e  R.   Now 
w  w 

+  +       + 

A  (x)  n  M  7^  0  and  from  Lemma  5.2,  A  (x)  =  A  (xt).   Hence 

A'''(xt)n  M  7^  0  which  implies  that  xt  e  A^(M)  .   Since  t  e  R 

was  arbitrary,  we  have  that  A,  (M)  is  invariant.   The  proof 

w 

of  the  invariance  of  A(M)  is  entirely  analogous.         // 
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Using  these  sets,  we  are  now  ready  to  define  the  various 
types  of  stability  to  be  investigated  in  this  section. 

DEFINITION  5.^   A  given  set  M  is  said  to  be 

1.  a  weak  attractor  if  A,  (M)  is  a  neighborhood  of  M, 

2.  an  attractor  if  A(M)  is  a  neighborhood  of  M, 

3.  stable  is  every  neighborhood  U  of  M  contains  a 
positively  invariant  neighborhood  V  of  M, 

11.   asymptotically  stable  if  it  is  stable  and  is 

an  attractor, 
5.   unstable  if  it  is  not  stable. 

Geometrically,  M  is  stable  if  given  any  neighborhood  U 
of  M,  there  is  a  neighborhood  VCU  such  that  any  trajectory 
that  is  within  V  at  time  t  =  0  remains  within  V  for  all 
t  >  0.   If,  in  addition,  all  trajectories  approach  M  as 
t  increases,  then  M  is  asymptotically  stable.   Finally 
M  is  unstable  if  every  neighborhood  of  M  contains  trajec- 
tories which  are  arbitrarily  close  to  M  at  time  t  =  0  but 
which  eventually  leave  the  neighborhood  for  all  time.   These 
three  cases  are  illustrated  in  the  diagram  below. 


Siymp 


n«4aUe 


im 


We  have  seen  specific  instances  of  these  cases  in 
several  of  our  examples.   For  instance  in  the  pendul 


urn 


problem,  (Example  2.1),  the  point  (q)  is  stable;   in  the 
damped  spring  problem  (Example  2.2),  the  point  (^)  is 
asymptotically  stable  and  in  Example  3.2,  if  the  system  is 
restricted  to  the  interior  of  the  unit  circle,  the  point 
(j^)  is  unstable.   In  this  latter  example,  if  x  is  any  point 
on  the  unit  circle,  then  the  point  x  is  a  weak  attractor 
although  it  is  not  stable.   In  order  to  see  an  example  of 
an  attractor  which  is  not  stable,  consider  the  following 
example: 

Let  the  planar  dynamical  system  be  defined  by  the 
differential  equations  (in  polar  coordinates): 

r  =  r(l  -  r) , 
e  =  sin^  (e/2). 

The  trajectories  of  this  system  are  sketched  below.   They 
consist  of  two  critical  points  at  (q)  and  (q),  a  trajectory 
on  the  unit  circle  with  (q)  as  the  positive  and  negative 
limit  set  of  all  points  on  the  unit  circle.   All  other  points 
have  the  point  (J)  as  their  positive  limit  set  and  hence 
the  point  (q)  is  an  attractor.   It  is  not  stable  because 
every  neighborhood  will  contain  part  of  the  unit  circle 
where  6  >  0.   Hence  no  positively  invariant  subset  of  this 
neighborhood  can  exist. 
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We  are  now  ready  to  prove  some  theorems  about  the 
concepts  just  introduced. 

THEOREM  5.5   I_f  M  is_  a  weak  attractor  (respectively,  an 
attractor) ,  then  the  set  A  (M)  (respectively ,  A(M))  i^  an 
open  neighborhood  of  M. 

PROOF  Let  N  denote  either  of  the  sets  A^(M)  or  A(M).   Then 
N  is  an  invariant  neighborhood  of  M.   Consequently,  bdy  N 
is  invariant  by  Theorem  2.10,  is  disjoint  from  M,  and  is 
closed.   Then  for  each  x  £  N,  A  (x)  0  M  5^  0,  whereas  for 
each  X  e  bdy  N,  A  (x) C  bdy  N  since  bdy  N  is  invariant 
and  closed.   Since  bdy  N  H  M  =  0,  we  conclude  that 
N  n  bdy  N  =  0.   Thus  N  is  open.  // 

We  note  that  if  a  set  M  is  stable  then  it  is  the  inter- 
section of  positively  invariant  neighborhoods  and  M  is 
therefore  positively  invariant.   Immediately  from  this  we 
see  that  if  the_  singleton  (x)  is_  stable,  then  {x  }  must  be 
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positively  Invariant  so  consequently ,  x  must  be  a  critical 
point . 

The  next  two  theorems  show  different  characterizations 
of  stable  sets  but  first  we  establish  the  following  lemma 
and  Its  corollary  which  do  not  depend  on  the  local  compact- 
ness of  X. 

LEMMA  5 . 6   Let  X  be  an  arbitrary  metric  space.   Let  x  e  X 

+  +       + 

and  w  G  A  ( x ) .   Then  J  ( x ) C  J  ( w ) . 

PROOF  Given  w  e  A  (x)  and  any  y  e  J  (x),  there  exist 
sequences  ^^^'^  j  "^n  "^  "''°°'  xt^->w,  and  {tM  and  (x^^), 

X  -*■  X .  t '  -»-  +~,  x  t'  -»■  y.   We  may  assume.  If  necessary  by 
n    '   n     '   n  n 

choosing  subsequences,  that  t^!^  -  t^  >  n  for  each  n. 
Consider  for  each  fixed  k,  k  =  1,  2,  ...,  the  sequence 

(x  t'}.   By  the  continuity  axiom  ^^'^{^'^  ^'^{^>   ^   =  ^>    ^»  

We  may,  therefore,  without  loss  of  generality,  assume  that 

for  each  fixed  k,  d(xTj^,x^Tj^)  <  1/k  for  n  >  k.  '  This  shows 

that  X  T'  -»-  w,  because  d(w,x^T^)<_d(w,xT^)  +  d(xTj!^,x^T^)  < 

d(w,XT')  +  1/n.   Now  notice  that  x^t^  =  x^T^(t^  -t^), 

and  x^t^  ->  y,  x^t^  ->  w,  and  t^!^  -  t^  >  n.   Hence 

y  e  J''"(w).   As  y  e  j'*'(x)  was  arbitrary,  we  have 

J'*'(x)  C  J'^(w),  and  the  lemma  Is  proved.  // 

COROLLARY   Given  M  and  x  e  A^(M) ,  then  J  (x)  <^  J    (M) C   D  (M) 

PROOF   Indeed  for  any  x  e  A^(M)  ,  A'^(x)  O  M  /   0.   Let  w 
be  any  element  In  A^(x)  H  M.   Then  by  the  previous  leirjna, 
j"^(x)C  j'^'Cw).   But  j'^(w)C  J"^(M)  and  we  know  that 
J"^(M)  C  d"^(M).   Hence  j'^(x)  C  j'^'CM)  d  d"^  (M) .  // 
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We  now  give  an  important  characterization  of  stability 
of  a  set  M. 

i 

THEOREM   5.7      A   set_  M  is_  stable    if   arid_  only   if  d'''(M)    =   M. 

I 

PROOF  Let  D  (M)  =  M  and  suppose  if  possible  that  M  is  not 

stable.   Then  there  is  an  e  >  0,  a  sequence  {x  },  and  a 

n  * 

sequence  {t  },  with  t   >  0,  d(x^,M)  -^  0,  and  d(x  t  ,M)  >  e. 
n         n—       n       '        n   n  ~ 

We  may  assume  without  loss  of  generality  that  e  >  0  has 
been  chosen  so  small  that  S[M,e]  and  hence  H(M,e)  is  com.pact 
where  H(M,e)  =  {x  e  X:  d(x,M)  =  e}.   This  is  possible  since 
X  is  locally  compact.   Further,  we  may  assume  that 
x  ->  X  c  M.   We  can  now  choose  a  sequence  {x  },  0  £  x   £  t  , 
such  that  x^x^  e  H(M,e),  n  =  1,  2,  Since  H(M,e)  is 

compact,  we  may  assume  that  x  x   ->  y  e  H(M,e).   Then  clearly 

+       + 
y  c  D  (x)C  D  (M),  but  y   f^  ¥[.      This  contradiction  shows  that 

M  is  stable.   Conversely,  assume  that  M  is  stable.   Then 

given  any  neighborhood  U  of  M  there  is  a  positively  invariant 

neighborhood  V  of  M  with  Y  C   \J .      Since  for  any  x  e  M, 

+  + 

D  (x)  C  cl  WR  for  any  neighborhood  W  of  x,  we  get 

D  (x)  C  cl  V  since  V  is  positively  invariant.   Thus 
D  (M)  C  cl  U  for  any  neighborhood  U  of  M.   Hence 
d''"(M)C  r\icl   U:  U  is  a  neighborhood  of  M}  =  M  as  M  is 
compact.   Since  M  CD  (M)  always  holds,  we  have  D  (M)  =  M.  // 

The  following  remark  shows  what  results  when  we  combine 
several  of  the  concepts  under  discussion. 
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REMARK   I£  M  is_  stable  and  is  a  weak  attractor,  then  M  ls_ 
an  attractor  and  consequently  asymptotically  stable . 

PROOF  Since  M  is  a  weak  attractor,  A  (M)  is  a  neighborhood 

w 

of  M.  We  need  to  show  that  A,,(M)C  A(M)  .   Let  x  c  A  CM)  . 

w  w 

Then  a'^'Cx)  7^  0  since  a'^(x)  H  M  7^  0  and  a'^(x)  C  J"^(x)  C  d'^(M) 
by  the  previous  corollary.   Since  M  is  stable  we  have 
D  (M)  =  M  by  Theorem  5.7  and  consequently  A  (x)  C  M.   Thus 
x  c  A(M)  and  A(M)  is  a  neighborhood  of  M.   Thus  M  is  an 
attractor.   Since  M  is  stable  and  is  an  attractor,  M  is 
asymptotically  stable. 

The  various  implications  that  have  been  proved  thus  far 
are  summarized  in  the  following  diagram. 


sla\>l< 
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We  would  now  like  to  confine  our  attention  to  the  study 
of  stability  for  systems  of  differential  equations.   We 
shall  begin  by  considering  autonomous  systems  of  constant 
coefficient  linear  equations;  that  is,  systems  of  the  form 

X  =  AX 
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where  A  Is  an  nxn  real  valued  matrix.   The  origin  in  R^ 
is  a  critical  point  for  such  a  system,  and  v;hen  zero  is 
not  an  eigenvalue  for  A  (so  that  is  v/hen  A  is  non-singular) 
the  origin  is  the  only  critical  point.   The  solutions  of 
this  system  can  be  described  in  terms  of  the  eigenvalues 
of  A  in  the  following  manner. 

Let  X  =  a  +  bi  be  a  complex  eigenvalue  for  the  nxn 
real  matrix  A,  and  let  E  be  an  eigenvector  in  C  belonging 
to  X.   Then  the  functions 


X  (t)  =  e^^CGj^cos  bt  +  H^sin  bt), 
X^(t)    =  e^-^(H^cos  bt  -  G^sin  bt). 


are  linearly  independent  solutions  of  X  =  AX,  where 


^x  ■*"  ^X 
G^  =  -^-2 — —  and 


i(E,  -  E,  ) 
H  = —  .      [Kreider,  Kuller,  Ostberg, 

p.  255] 


X       2 


Now  when  a  ±  bi,  b  >  0,  are  complex  conjugate  eigenvalues 
of  multiplicity  m,  this  system  has  2m  linearly  independent 
solutions  constructed  from  functions  of  the  form 


G,e^^t^cos  bt 


and 
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H.e  t  sin  bt. 


where  k  is  an  Integer  such  that  0  <_  k  £  m  -  1.   Then 
assuming  t  >  0, 

llG^e^S^cos  bt||  <  IIGJI  e^S^, 

and  since  I  I G^ I  I   is  a  positive  constant, 

lim  llG.e^^t^cos  bt I  I  =  0   if  a  <  0. 

On  the  other  hand  ||G^e^^t  cos  bt | |  is  unbounded  if  a  >  0. 
Similarly, 

lim  ||H  e^^t^sin  bt | |  =  0   if  a  <  0, 

but  llH^e^^t^sin  bt | |  is  unbounded  if  a  >  0. 

Thus  every  trajectory  of  this  system  arising  from  a 
pair  of  complex  conjugate  eigenvalues  for  A  will  approach 
the  origin  of  R  as  t  -»-  +~  i_f  the  real  part  of  these  eigen- 
values is  negative ;  it  will  depart  arbitrarily  far  from 
the  origin,  if  the  real  part  of  these  eigenvalues  is  positive 
In  the  first  case,  the  trajectories  exhibit  the  properties 
required  for  asymptotic  stability  whereas  in  the  second 
case  they  exhibit  instability.   This  reasoning  clearly 
applies  equally  well  to  trajectories  arising  from  real 
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eigenvalues  for  A.   It  remains  to  consider  the  case  where 
A  admits  a  pair  of  pure  imaginary  eigenvalues  ±  bi.   If 
these  eigenvalues  have  multiplicity  one,  the  corresponding 
solutions  of  the  system  are  constructed  from  the  functions 


G^cos  bt   and   H  sin  bt 


and  the  fact  that  cos  bt  and  sin  bt  are  bounded  as 
t  -»-  +00  yet  do  not  tend  to  zero  implies  that  the  resulting 
trajectories  have  the  properties  required  for  stability, 
but  not  asymptotic  stability.   On  the  other  hand,  if  the 
multiplicity  of  ±bi  is  greater  than  one,  the  system  has 
solutions  involving  functions  of  the  form 


k  k 

G^t  cos  bt   and   H  t  sin  bt 


with  k  >_  1.   When  this  happens,  the  origin  is  unstable  since 

|t  cos  bt I    and    |t  sin  bt | 

are  unbounded.   Thus  we  have  proved  the  following  theorem. 

THEOREM  5.8   LfX=AXis_an  nxn  linear  autonomous  system 
whose  coefficient  matrix  is  nonsingular ,  then  the  origin 


in  R^  is 


1.   asymptotically  stable  if  the  real  parts  of  all 
of  the  eigenvalues  of  A  are  negative; 
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2.  stable,  but  not  asymptotically  stable ,  if  A 
has  at  least  one  pair  of  pure  lrriap;lnary  eigen- 
values of  multiplicity  one ,  no  pure  Imag-'inary 
eigenvalues  of  multiplicity  exceeding  one ,  and 
no  eigenvalues  with  positive  real  parts ; 

3.  unstable  otherwise . 

Consider  the  following  plane  autonomous  systems: 


1.   X  =  y 


2.   X  =  y 


3.   X  =  -X 


y  =  X 


y  =  -X 


y  =  -y 


The  coefficient  matrix  for  the  first  system  is 


A  = 


with  eigenvalues  X  =  ±1.   Hence  the  origin  is  unstable 
since  one  of  the  eigenvalues  has  a  positive  real  part. 
For  the  second  system. 


A  = 


-1 


This  system  is  precisely  that  of  Example  2.3  and  has 
eigenvalues  X  =  ±i.   Thus  the  origin  is  stable  but  not 
asymptotically  stable.   This  is  exactly  as  we  would  expect 
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by  recalling  the  geometry  of  Example  2.3 
For  the  third  system. 


A  = 


-1 


-1 


This  system  has  X  =  -1  as  its  only  eigenvalue  and  thus 
the  origin  is  asymptotically  stable. 

The  above  results  can  be  used  to  obtain  a  complete 
description  of  the  trajectories  of  a  linear  plane  autonomous 
system  of  the  form 

X  =  ax  +  by 
y  =  ex  +  dy 


where  the  matrix 


A  = 


is  nonsingular.   This  description  depends  upon  the  nature 
of  the  eigenvalues  ^,^2  °^  ^  which  are  the  roots  of  the 
characteristic  equation 

X^  -  (a  +  d)X  +  (ad  -  be)  =  0. 
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since  the  discriminant  of  this  equation  is  given  by 

A  =  (a  +  d)^  -  4(ad  -  be)  =  (a  -  d)^  +  ilbc, 

it  follows  that  X  and  A  will  be 

1.  real  and  distinct  if  A  >  0, 

2.  real  and  equal  if  A  =  0, 

3.  complex  conjugates  if  A  <  0. 

Moreover,  when  X^   and  X^   are  the  roots  of  this  system,  we 
have  (X  -  X^)(X  -  X2)  =  0.   Thus  we  have 


X^  -  (X^  +  X2)X  +  X^X^  =  0, 


so 


X^  +  Xp  =  a  +  d,  ^1^2   ="  Q'd  -  be  =  det  A 

Thus  when  A  >_  0,  X,  and  Xp  will  have  the  same  sign  if  and 
only  if  ad  -  be  >  0,  and  will  then  be  positive  or  negative 
according  as  a  +  d  is  positive  or  negative.   On  the  other 
hand,  if  A  <  0,  X,  =  a  +  6i,  X^   =   ol   -   3i,  B  >  0,  and 
a  +  d  =  2a. 

In  this  case  X^  and  Xp  will  be  pure  imaginary  if  and  only 
if  a  +  d  =  0.   Otherwise  they  will  have  a  nonzero  real  part 
which  agrees  in  sign  v/ith  a  +  d. 

Now,  combining  these  observations  with  Theorem  5.8,  we 
have  the  following  description  of  the  origin  for  a  plane 
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autonomous  linear  system  with  constant  coefficients. 


1. 


2.   A  >  0 


asymptotically  stable 

unstable 

stable 


ad  -  be  >  0   and  a  +  d  <  0   asymptotically 

stable 
ad  -  be  >  0   and   a  +  d  >  0  unstable 


3.   A  =  0 


ad  -  be  <  0 
"ad  -  be  >  0 
ad  -  be  <  0 


and  a  +  d  <  0 


unstable 

asymptotically 

stable 

unstable 


Thus  far  we  have  not  discussed  the  cases  where  A  is 
singular.  These  cases  are  rather  uninteresting  but  for 
completeness,  they  are  presented  here. 

Suppose  A  is  the  2x2  zero  matrix.   Then  X  =  0  for  all 

2  2 

X  e  R  which  implies  that  each  point  in  R  is  a  critical 

point  and  consequently  each  point  is  stable  since  every 
point  within  some  neighborhood  of  a  point,  remains  within 
that  neighborhood  as  t  -^  +<». 

Suppose. A  is  a  2x2  matrix  with  zero  as  its  only  eigen- 
value.  Then  A  must  be  of  the  form 


A  = 


-a 


-a 


where  a  is  some  real  number. 
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To  find  the  eigenvectors  X^  associated  v;ith  A,  v;e  set 
AXq  =  0.   Then  all  eigenvectors  are  critical  points  and 
are  given  by  the  equation 


»=[; 


-a 
-a 


ax,  ■"  ax,^ 


which  implies  that  x,  =  Xp .   Thus  all  points  on  the  line 
through  the  origin  with  slope  =  1,  are  critical  points. 
Now 


X,  —  ax-|  —  ax.-j 


Xq  —  ax-|  ~  aXp 


so  we  have  x^  =  Xp  which  implies  that  x^  =  x^  +  c  where  c 
is  any  constant.   Thus  the  remaining  trajectories  also  have 
slope  =  1,  and  intersect  the  Xp-axis  depending  on  c. 
Their  directions  are  shown  in  the  sketch  below.   It  is  clear 
that  all  critical  points  are  unstable. 
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Thus  far  v/e  have  studied  only  linear  systems.   To 
study  non-linear  systems  we  need  the  following  idea  of  a 
linear  approximation.   Consider  the  system 

X  =  F(X) 

where  X  is  a  vector  in  r".   We  shall  show  that  v/henever  F 
is  a  function  of  class  C   (that  is,  P  has  continuous  first 
partial  derivatives)  in  a  region  Q   of  R^  containing  the 
origin,  then  we  may  replace  F  by  a  function  of  the  form 

JX  +  G(X), 

where  J  is  an  nxn  matrix  and  G(X)  is  "small"  in  comparison 
with  JX  when  ||x||  is  small.   Precisely,  we  state 


DEFINITION  5.9   Let  F  be  as  above  and  let  Xq  be  a  point  in 


fi.   Then  a  function  L  =  L(X)  is  said  to  be  a  linear 

^0 


approximation  to  F  at  X-,  if  L  is  linear  on  R  and 


[F(X)  -  F(X^)]  -  [L(X)  -  LCX^)] 
lim —  =  0 


IIx-XqII-o  ||x  -  XqII 


Moreover,  the  existence  and  uniqueness  of  the  linear 
approximation  can  be  guaranteed  if  F  is  of  class  C  . 
[Ostberg,  Kreider,  Kuller  pp.  4l7,  ^l8].   This  linear 
approximation  is  in  fact  the  nxn  Jacobian  matrix  of  F  at 
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every  point  X  .   Specifically, 


J  = 


8P. 

9x: 


8P 


n 


ax. 


8F. 

9r" 


n 


9P 


n 


dx 


n 


X  =  X, 


where 


F(X)  = 


F^(X) 


F  (X) 
n 


and  F^:R"  ->  R  for  all  i. 


To  gain  a  better  understanding  of  the  linear  approxima- 
tion, we  give  some  examples. 

We  find  the  linear  approximation  about  („)  for  each  of 


the  following  systems : 
1.   X  =  P^(x,y) 


y  =  P2(x,y) 


where  P,  and  P^  are  polynomials  with  P^(0,0)  =  P2(0,0)  =  0 
Then  F^(x,y)  =  P^(x,y)  and  F2(x,y)  =  P2(x,y).   Now  if 
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P^(x,y)  -  a^x  +  b^y  +  c^xy  +  a^x  +  b^y  +  higher-order  terms. 


P^Cxjy)  =  d^x  +  e^y  +  f^xy  +  d^x   +  e^y   +  higher-order  terms, 


we  have 


9P. 


— ^1  "^  "^1^  "^  2apX  +  higher-order  terms   and 


8F^ 

=  t)-,  +  c^  X  +  2b^y  +  higher-order  terms 

dy     11      d 


Now  evaluating  these  partials  at  the  point  („) ,  we  have 


9F  aF 

TT —  =  a.    and  - —  =  b^ 
8x     1        9y     1 


Similar  results  hold  for  the  partials  of  F^.   Thus  we 
replace  F(X)  by  JX  +  G(X)  where 


J  = 


^1   ^ 


^1   ^ 


,  and  G(X)  = 


2       2 
c  xy  +  a^x  +  b2y   + 


2      2 
f^xy  +  d^x   +  e^y   +  . 


Hence  the  linear  approximation  is  given  by  the  following 
system: 

X  =  JX. 
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2.   X  =  sin  X  +  e^  - 


y  =  xy 


Now,  sin  X  and  e^  have  the  following  Taylor  series 
expansions : 


x^        x5 
sin  X  =   X  -  ~  +  ^ 


7! 


+   - 


e^=l+y+|y+    .. 


Hence   we   have   x  =   x  +  y   - 


3  2 

X  V 

57-  +  |y  +  higher-order  terms 


y  =  xy 

and  applying  the  procedure  used  in  the  previous  example,  we 
have  that  the  approximating  system  is  given  by 

X  =  JX 


where 


J  = 


1    1 


0     0 


To  analyze  this  system  as  to  stability,  we  note  that  the 
matrix  J  is  singular  and  hence  our  theorem  does  not  apply 
Thus  we  must  analyze  J  using  the  same  method  we  used  in 
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the  two  previous  cases  of  singular  matrices.   We  see  that 
J  has  eigenvalues  X  =  0,1  and  hence  does  not  fit  into  the 
category  of  the  other  two  examples.   However,  it  is  an 
easy  task  to  determine  that  all  eigenvectors  lie  on  the 
line  of  slope  -1  through  the  origin  and  are  all  critical 
points.   From  the  approximating  system,  it  is  easily  seen 
that  y  =  0  and  hence  all  trajectories  not  on  the  line  of 
critical  points  are  parallel  to  the  x-axis.   Their  direc- 
tions are  easily  determined  from  the  equations.   They  are 


shown  in  the  following  diagram, 

y 


-<r 


-v 


It  is  clear  that  each  critical  point  is  unstable. 

One  might  question  now  whether  the  complete  theory  we 
developed  for  linear  systems  also  applies  to  the  approxima- 
ting systems  for  non-linear  systems.   Such  is  the  case  for 
asymptotically  stable  and  unstable  systems.   That  is,  if 
the  linear  approximating  system  is  unstable,  then  the  non- 
linear system  is  unstable.   However,  no  definite  conclusion 
can  be  made  about  stable  systems.   It  is  recalled  that  for 
a  system  to  be  stable,  it  could  have  no  eigenvalues  with 
positive  real  parts.   A  non-linear  system  may  have 
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eigenvalues  with  small  positive  real  parts  v/hich  do  not 
appear  in  its  approximating  system.   These  facts  are 
proved  in  the  next  section  and  are  presented  here  Just 
for  information. 
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VI.   LIAPUNOV  STABILITY 

The  basic  feature  of  stability  theory  a  la  Liapunov  is 
that  one  seeks  to  characterize  stability  and  asymptotic 
stability  of  a  given  set  in  terms  of  a  non-negative  real- 
valued  function  defined  on  a  neighborhood  of  the  given  set. 
This  technique  was  devised  by  the  Russian  mathematician 
A.  A.  Liapunov,  and  is  known  as  Liapunov 's  direct  or  second 
method.   It  is  based  on  the  well-known  fact  that  a  physical 
system  loses  potential  energy  in  a  neighborhood  of  a  point 
of  stable  equilibrium.   More  precisely,  a  point  of  stable 
equilibrium  for  a  physical  system  is  a  point  at  v/hich  the 
potential  energy  of  the  system  has  a  local  minimum.   This 
fact  is  known  in  physics  as  Lagrange's  theorem.   We  proceed 
with  the  following  lemma. 

LEMMA  6. 1   Let  the  phase  space  X  be_  arbitrary  and  let  K  C  X, 
Let  ({)  be  any  continuous  real-valued  function  defined  on  K 
such  that  (j)(xt)  <_  cj)(x)  whenever  x[0,t]  C  K,  t  ^0.   Then 
if  for  some  x,  cl  (y  (x))  C  K,  we_  have  <t>(y)    =   (piz)    for 
every  y ,z  e  A  (x) . 

PROOF  Now  there  are  sequences  {t^}  and  {t^}  in  R  such 

that  t   -»-  +0O,  X   -»-+»,  and  xt   -^  y,  xx   ->  z.   We  may  assume 

by  taking  a  subsequence  that  x   >  t^  for  each  n.   Then 

clearly  <^(xt^)    >  (})(xx^)  since  xx^  =  xt^(x^  -  t^), 

X   -  t   >  0,  and  x  t  [0,x   -  t  ]  C  K.   Thus  proceeding  to 
n    n    '      n  n   '  n    n 

13^ 


the  limit  we  have  by  the  continuity  of  ({),  {|)(y)  >_  i>(z)  . 

A  similar  argument  shows  that  <\){z)    >  (t)(y).   Thus  the  lemma 

is  proved.  // 

The  best  known  result  on  asymptotic  stability  is  the 
following. 

THEOREM  6.2   A  compact  set  M  C  X  is_  asymptotically  stable 
if  and  only  if  there  exists  a  continuous  real-valued 
function  E  defined  on  a  neighborhood  N  of  M  such  that 

1.  if  X  e  M,  we_  have  E(x)  =  0  and  if  x  ^   M,  we 
have  E(x)  >  0, 

2.  for  x  ^f  M,  t  >  0  and  x[0,t]  C  N,  we_  have 
E  ( xt )  <  E  ( X ) . 

PROOF   Assume  that  a  function  E  as  required  is  given. 

Choose  a  >  0  such  that  S[M,a]  C  N  and  is  compact.   Let 

m  =  min{E(x):x  e  H(M,a)}  where  H(M,a)  =  {x  e  X:d(M,x)  =  a}. 

Such  a  minimum  exists  because  E  is  continuous  and  H(M,a) 

is  a  closed  and  bounded  set.   By  statement  1  and  the 

continuity  of  E  we  have  m  >  0.   Set 

K  =  {x  G  S[M,a]:E(x)  <  m}.   Then  K  is  compact  and  because 

of  statement  2,  K  is  positively  invariant.   This  establishes 

that  M  is  stable  as  K  is  a  positively  invariant  neighborhood 

of  M.   To  see  that  M  is  an  attractor,  choose  any  compact 

positively  invariant  neighborhood  K  of  M  with  K  C  N.   Then 

for  any  x  e  K,  y"^(x)  C  K  and  since  K  is  compact, 

cl  (y^(x))C  K.   Hence  cl  (y'^(x))  is  compact  and  we  have 
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0  ?^  A  (x)  C  K,  and  Lemma  6.1  shows  that  E  is  constant  on 
A  (x).   But  this  shows  by  statement  2  that  a'^(x)C  M.   Thus 
M  is  an  attractor  and,  consequently,  asymptotically  stable. 

Conversely,  let  M  be  asymptotically  stable  and  A(M) 
its  region  of  attraction.   For  each  x  e  A(M)  define 

<{)(x)  =  sup{d(xt,M):  t  >_  0}. 

Indeed  (|)(x)  is  defined  for  each  x  e  A(M)  because  if 
d(x,M)  =  a,  then  there  is  a  T  >  0  with  x[T,+oo)  C  S(M,a). 
This  is  because  x  e  A(M).   Thus 

(J)(x)  E  sup{d(xt,M):  0  <_  t  ;^  T}. 

Since  d(xt,M)  is  a  continuous  function  of  t,  (|)(x)  is  defined 
Now  (j)(x)  has  the  following  properties:   (})(x)  =  0  for  x  c  K, 
<t)(x)  >  0  for  X  2f  M,  and  {j)(xt)  <_  (|)(x)  for  t  >  0.   This  is 
clear  when  we  remember  that  M  is  stable  hence  positively 
invariant,  and  that  A(M)  is  invariant.   Thus  if  4)(x)  is 
defined  for  any  x  e   A(M)  it  is  defined  for  all  xt  with 
t  e  R  .   We  further  claim  that  this  (|)(x)  is  continuous  in 
A(M).   Stability  of  M  implies  continuity  of  (|)(x)  on  M. 
For  x  ^  M,  this  proof  requires  a  knowledge  of  uniform 
attractors  and  it  will  not  be  presented  here.   However, 
this  function  is  continuous  in  A(M)  [Bhatia,  Szego  p.  6?]. 
This  function  may  not  be  strictly  decreasing  along 
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trajectories  in  A(M)  which  are  not  in  M  and  so  may  not 
satisfy  statement  2  of  the  theorem.   However  we  can  obtain 
the  desired  function  by  setting 


E(x)  =   /  (|)(xT)exp(-T)  dx. 
0 


Clearly  E(x)  is  continuous  and  satisfies  statement  1  in 
A(M).   To  see  that  E(x)  satisfies  statement  2,  let  x  iff  M  and 
t  >  0.   Then  E(xt)  <  E(x)  holds  because  <t>{xt)    <  <^(x)    holds. 
To  rule  out  E(xt)  =  E(x),  observe  that  in  this  case  we 
must  have  (|)(x(t  +   t))    =   (})(xt)  for  all  t  >  0.   Thus  in 
particular,  letting  x  =  0,  t,  2t,  ...  we  get  <^{x)    =  (|)(x(nt)), 
n  =  1,  2,  ...   But  asymptotic  stability  of  M  implies  that 
for  X  e  A(M),  d(xt,M)  -»-  0  as  t  h-  +«.   Thus 
<t)(x(nt))  ->-  0  as  n  ^-  «,  as  (t)(x)  is  continuous.   This  shows 
that  (j)(x)  =  0.   But  as  x  jef  M,  we  must  have  4)(x)  >  0  and 
we  have  a  contradiction.   Hence  E(xt)  <  E(x)  for  x   gi  M 
and  t  >  0.   Thus  statements  1  and  2  are  satisfied.       // 

This  theorem  says  nothing  about  the  size  of  the  region 
of  attraction  of  M.   Thus  if  a  function  E(x)  is  known  to 
exist  in  a  neighborhood  N  of  M,  we  need  not  have  either 
N  a  A(M)  or  A(M) C   N.   This  means  that  the  above  theorem 
cannot  immediately  be  stated  as  a  theorem  on  global 
asymptotic  stability.   The  following  example  helps  to 
illustrate  this  point. 
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EXAMPLE  6.3   Consider  a  dynamical  system  defined  In  the 
x,y-plane  by  the  differential  equations 

X  =  f(x,y), 

y  =  -y, 

where 


f(x,y)  = 


2  2 
X  Is  X  y   >  1 


?  ?         2  2 
2x^y   -  X  If  X  y   <  1 


The  trajectories  of  this  system  are  sketched  below, 


The  origin  (q)  is  asymptotically  stable,  with  the  set 
{(^):x^y^  <  1}  as  its  region  of  attraction.  Consider  now 
the  function 
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E(x,y)  =  y^  + 


2    -    -2 


X 


1  +  x^ 


This  function  satisfies  the  conditions  of  Theorem  6 . 2  in 
the  whole  plane.   To  see  this  we  note 

dE  _  8E  dx  .  9E  dy   ^^  ^  ,        *        :^v  ot? 

dt  -  87  dF  "■  87  dt   ^^i^h  implies  E(x,y)  =  ||f(x,y)  -  y  ||. 

It  can  be  verified  that  E(x,y)  <  0  for  all  ('')  7^  (°),  which 
implies  that  E  satisfies  the  conditions  of  the  theorem  in 
every  neighborhood  of  the  origin.   But  not  every  neighbor- 
hood of  the  origin  is  contained  in  the  region  of  attraction, 
nor  does  it  contain  the  region  of  attraction. 

In  view  of  this  example,  we  wish  to  establish  conditions 
necessary  to  ensure  global  asymptotic  stability.   This  is 
done  by  starting  with  the  following  definition. 

DEFINITION  S.h      A  continuous  real-valued  function  E(x) 
defined  on  a  set  N  C  X  will  be  said  to  be  uniformly  unbounded 
on  N  if  given  any  a  >  0  there  is  a  compact  set  K  C  N, 
K  7^  N,  such  that  E(x)  >_  a  for  x  ^  K. 

The  following  three  theorems  will  be  stated  v;ithout 
proof  since  the  proofs  are  very  laborious.   The  proofs  may 
be  found  in  Bhatia  and  Szego,  pp.  70  -  73. 

THEOREM  6 . 5   Let  M  C  X  be_  a  compact  asymptotically  stable 
set  and  let  M  be_  invariant .   Then  there  exists  a  continuous 
uniformly  unbounded  function  E(x)  defined  on   A(M)  such 
that 
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1.  for  X  e  M,  E(x)  =  0  and  for  x  £f  M,  E(x)  >  0, 

2.  for  all  X  e  A(M)  and  t  e  R,  E(xt)  =  e"'^E(x). 

THEOREM  6.6   If  M  C  X  is_  any  compact  asymptotically  stable 
set_,  there  exists  a  continuous,  uniformly  unbounded  function 
E(x)  on  A(M)  such  that 

1.  for  X  e  M,  E(x)  =  0  and  for  x  i   H,  E(x)  >  0, 

2.  for  X  ;?f  M  and_  t  >  0,  E(xt)  <  E(x). 

THEOREM  6.7  Let  M  C.  X  be_  com.pact  and  let  there  exist  a 
continuous  uniformly  unbounded  function  E(x)  defined  on 
an  open  neighborhood  N  of  M  such  that 

1.  for  X  e  M,  E(x)  =  0  and  for  x  i   M,  E(x)  >  0, 

2.  for  X  jzf  M,  t  >  0  ajid  x[0,t]C  N,  E(xt)  <  E(x). 
Then  M  is_  asymptotically  stable  and  N  C  A(M) .   If  in 

addition i  any  condition  guaranteeing  the  invariance  of  N 
holds,  then  N  =  A(M). 

These  three  theorems  are  very  similar  in  their  conclu- 
sions but  we  can  see  that  the  hypotheses  are  becoming  more 
general.   The  next  results  on  global  asymptotic  stability 
follow  from  these  theorems. 

THEOREM  6.8   A  compact  invariant  set  M  C  X  is  globally 
asymptotically  stable  if  and  only  if  there  exists  a  continu- 
ous uniformly  unbounded  function  E(x)  defined  on  X  such  that 

1.  for  X  e  M,  E(x)  =  0  and  for  x  /zf  M,  E(x)  >  0, 

2.  for  all  X  e  X  and  t  e  R,  E(xt)  =  e~  E(x). 
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PROOF   The  sufficiency  follows  from  Theorem  6.7  since  X  is 
an  invariant  neighborhood  of  M.   The  necessity  follows  from 
Theorem  6.5.  // 

THEOREM  6.9   A  compact  set  M  C  X  is_  globally  asymptotically 
stable  if  and  only  if  there  exists  a  continuous  uniformly 
unbounded  function  E(x)  defined  on  X  such  that 

1.  for  X  e  M,  E(x)  =  0  and  for  x  i   M,  E(x)  >  0, 

2.  for  X  i!f  M  and  t  >  0,  E(xt)  <  E(x). 

PROOF  The  sufficiency  follows  from  Theorem  6.7  and  the 
necessity  follows  from  Theorem  6.6.  // 

These  theorems  allow  us  to  determine  whether  a  set  is 
asymptotically  stable  or  not  without  having  to  solve  the 
system.   However,  if  a  set  is  asymptotically  stable,  it 
may  not  be  a  trivial  problem  to  display  the  continuous 
uniformly  unbounded  Liapunov  function.   As  in  the  last 
section,  we  now  focus  our  attention  directly  on  systems 
of  differential  equations  in  order  to  apply  some  of  the 
above  theory. 

Thus  we  return  to  the  general  case  of  an  autonomous 
system 

X  =  F(X) 


defined  in  a  region  fi  of  r'^  with  a  critical  point  at  Xq 
in  9..      Without  loss  of  generality,  we  may  assume  that  the 


lill 


critical  point  is  the  origin  0  of  r"  for  if  it  is  not,  a 
simple  change  of  variable  x«  =  X  -  X^  will  translate  the 
origin  to  the  point  Xq  v/ithout  changing  the  form  of  the 
system. 

Again  we  seek  a  real-valued  function  E(X)  of  class  C''" 
defined  on  some  region  Q.   of  the  origin  such  that  E(X)  >  0 
and  E(X)  =  0  if  and  only  if  X  =  0.   V/e  will  refer  to  E  as 
the  energy  function  because  of  our  previous  physical 
interpretation.   If  X  =  X(X  ,t)  is  the  solution  of 

X  =  F(X) 


X(0)  =  Xq, 


9E 
the  time  rate  of  change  of  energy,  ^  along  the  trajectory 

defined  by  this  solution  is 


VE  •  X  =  VE-F  =   Z  1^  F, 

l=l^^i   ^ 


Thus  if  this  expression  is  negative  in  Q,  a  particle  moving 
along  any  trajectory  defined  by  the  system  will  be  dissipa- 
ting energy  as  it  intersects  the  surface  E(X)  =  c.   Conse- 
quently a  particle  which  enters  the  region  enclosed  by 
such  a  surface  can  never  gather  enough  energy  to  escape, 
and  its  trajectory  must  remain  in  that  region.   This  implies 
that  the  origin  will  be  asymptotically  stable,  for  the 
given  system. 
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These  remarks  motivate  the  following  definition. 

DEFINITION  6.10   Let  ^   he   a   region  of  R^  containing  the 
origin,  and  let  E  =  E(X)  be  a  real-valued  function  of 
class  C   in  U.      Then  E  is  said  to  be  positive  definite  if 

1.  for  all  X  in  fi,  E(X)  >_  0  and 

2.  the  equality  E(X)  =  0  holds  if  and  only  if 

X  =  0. 

If  in  addition  we  have  VE-F  =   E  | —  F,  <  0  everywhere 

i=l  ^^i   ^  " 

in  Uy   then  E  is  said  to  be  the  Liapunov  function  for  the 

autonomous  system  X  =  F(X). 


In  terms  of  this  definition,  we  can  now  form  some 
conclusions  about  the  stability  of  the  origin. 

THEOREM  6.11   The  origin  is  a  point  of  stable  equilibrium 
for 

X  =  F(X), 


F(Xq)  =  0 


provided  there  exists  a  Liapunov  function  E  for  the  system, 

PROOF  Assume  the  existence  of  the  Liapunov  function  for 
the  system  and  consider  the  set 

S[0,r]  =  {X  e  r'^:  d(0,X)  <  r} .   Let  r  be  so  chosen  that 
S[0,r]  lies  in  Q.      Since  S[0,r]  is  compact  and  E  is 
continuous  and  positive  definite  in  Q,    it  follows  that  E 
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achieves  its  minimum  value  on  S[0,r].   Let  m  be  this 
minimum  value.    Clearly  m  is  positive. 

Again  since  E  is  continuous,  there  exists  an  e  >  0  such 
that  I  |E(X)  -  E(0) I  I  <  m  whenever  |  | X  -  0 1  |  <  e ,  or  simply 
stated,  E(X)  <  m  whenever  ||X||  <  e. 
Consider  the  set  S(0,e).   Let  X^^ 
be  any  point  in  S(0,e)  and  let 
X(XQ,t)  be  the  unique  solution 
of  X  =  F(X)  whose  trajectory 
passes  through  X^  at  time 
t  =  0. 

Now  we  claim  that  the  trajectory 
must  remain  within  S[0,r]  for  all  t  >  0.   Assume  that  this 
is  not  the  case.   Then  ||X(XQ,t  )||  =  r  for  some  t^>   0, 
since  the  trajectory  is  a  connected  set.   Thus 
E(X(XQ,t  ))  >_  m  but  this  implies  that  E  is  increasing  with 
time  and  that  ^  >  0.   Now  15  =  ^ 


dx 


at 


at 


3x.  dt 


=  VE-F.'  Hence 


VE'F  >  0.   But  this  contradicts  the  assumption  that 
VE«F  <_  0  in  fi,  since  this  inequality  implies  that  E  is 
a  non-increasing  function  of  t  along  the  trajectory  in 
question.   Thus  the  trajectory  remains  v/ithin  S[0,r]  for 
all  t  >  0  and  hence  the  origin  is  stable  since  S[0,r]  is 
a  positively  invariant  neighborhood  of  the  origin.       // 

This  argument  suggests  that  if  E  actually  decreases 
along  trajectories  in  a  neighborhood  of  the  origin,  then 


Ikk 


the  origin  would  be  asymptotically  stable.   The  next  theorem 
uses  this  fact. 

THEOREM  6.12   I£  E  is_  a  Liapunov  function  for  X  =  F(X)  with 
the  property  that  -VE*F  is_  positive  definite  in  fi,  then 
the  origin  is  asymptotically  stable. 

PROOF  We  have  seen  that  there  exists  an  e  >  0  with  the 
property  that  E  is  a  non-increasing  function  along  any 
trajectory  X(X^,t)  with  | | Xq | |  <  e.   Since  E  is  positive 
definite  in  fi,  we  claim  that  E  approaches  a  non-negative 
limit  E^   along  this  trajectory  as  t  -»■  °°.   To  see  this,  we 
form  the  sequence  {E(X(Xq  ,t^) ) }  where  t^  ->  «'. 

Now  E(X(X  ,t  ))  assumes  values  in  [0,E(Xq)].   Since 
[0,E(Xp^)]  is  a  compact  set,  every  sequence  in  [0,E(Xq)] 
has  a  subsequence  which  converges  to  a  point  Eq  in  [0,E(Xq)] 
Moreover  E^  is  non-negative. 

We  will  be  done  if  we  show  that  the  limit  is  zero  when 
-VE-F  is  positive  .definite  in  ^.      Then  the  fact  that  E  is 
zero  only  at  the  origin  will  imply  that  the  trajectory 
approaches  the  origin  with  increasing  time. 

On  the  contrary,  suppose  E^  >  0.   Then  since  E  is 
continuous  there  exists  a  6  >  0  such 
that  ||E(X)  -  E(0)||  <  Eq  when 
||X||  <  6.   Consider  the  set 
S['o,6].   Now  the  trajectory 
X(XQ,t)  cannot  enter  the  sphere 
S[0,6]  for  if  it  did,  then 

1^5 


s[o,f;i 


E(X(X  ,t  ))  <  E   for  some  k  >  0.   But  since  E  is  a  non- 
Increasing  function,  it  could  never  leave  S[0,6]  and  thus 
could  not  approach  E„  as  a  limit. 

Now  since  F  is  continuous  and  E  is  of  class  C   in  Q, 
-VE'F  is  continuous  and  hence  assumes  its  minimum  value  on 
the  "annular"  region  6  <_  |  |X|  |  <  r  where  r  is  as  in  the 
proof  of  the  previous  theorem.   This  fact  is  true  because 
this  "annular"  region  is  compact.   Let  m  be  this  minimum 
value  of  -VE'F.   Since  -VE«F  is  positive  definite,  it 
follows  that  m  >  0.   Thus 


aE(X(X„,t)) 

-rr^ <  -m   for  all  t  >  0 

o  t 


Now 


t  8E(X(X  ,t)) 
E(X(XQ,t))  -  E(X(Xq,0))  =   /   ^ dT 


by  the  fundamental  theorem  of  the  calculus.   But 


t  8E(X(X  ,t))        t 
/   —i^ di  <   /   -mdx  =  -mt . 

I     ^-^       -  0 


Thus  we  have  E(X(XQ,t))  <  E(Xq)  -  mt . 

Since  the  right  hand  side  approaches  -»  as  t  ->  °°,  this 
implies  that  E(X(XQ,t))  <  0  for  all  t  >  t^  where 

t   =  —     but  this  contradicts  the  fact  that  E  is  positive 

m 


1^6 


definite  in  ^.      Thus  E^  =  0  and  the  origin  is  asymptotically 
stable.  // 

Thus  far  we  have  presented  theorems  only  on  stability 
or  asymptotic  stability  of  the  origin.   However  it  is  also 
important  to  know  under  what  conditions  the  origin  is 
unstable.   Consequently  we  present  the  following  instability 
theorem  known  as  Liapunov's  First  Instability  Theorem. 

THEOREM  6.13   Let  E  be_  a  real-valued  function  of  class  C 
in  ^   and  suppose  that  E(0)  =  0,  E  ls_  positive  definite  and 
E  is_  able  to  assume  positive  values  arbitrarily  near  the 
origin.   Then  the  origin  is  unstable. 

PROOF  Let  S[0,r]  and  S[0,g]  be  spheres  about  the  origin 
such  that  e  <  r  and  S[0,r]C  ^.   Since  E(X)  and  the  first 
partials  of  E(X)  are  continuous  in  f^ ,  we  have  that  E(X) 
is  bounded  in  S[0,r].   Hence  E(X)  <  k  for  all  X  e  S[0,r]. 
Choose  any  Xq  e  S[0,e]  such  that  E(Xq)  >  0  and  let  X(XQ,t) 
be  the  unique  solution  of  X  =  F(X)  whose  trajectory  passes 
through  Xq  at  time  t  =  0. 

Now  X(XQ,t)  must  eventually  cross  bdy  S[0,r]  for  if  not 
then  E(X(XQ,t))  approaches  a  limit  p  in  [E(XQ),k].   Since 
E(X(XQ,t))  ->  k  as  t  -^  0°,  it  follows  that  E(X(XQ,t))  ->  0. 
Now  E(X)  =  VE-F  and  since  E  and  F  are  continuous,  we  have 
that  E(X)  is  continuous.   Then  since  E(X(XQ,t))  ->■  0,  it 
follows  that  X(XQ,t)  -^  0  since  E  is  positive  definite. 
But  this  is  a  contradiction  since  if  X(XQ,t)  ->  0  as  t  ->  ~, 
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it  must  follow  that  E(X(X^,t))  ->  0  because  E  is  continuous. 
Thus  since  E(Xq)  >  0  we  have  E(X(XQ,t))  <  0  but  this  is 
a  contradiction.   Hence  X(X„,t)  must  eventually  cross  the 
boundary  of  S[0,r]  and  thus  the  origin  is  unstable.      // 


The  important  feature  of  these  theorems  is  that  they 
enable  us  to  determine  the  stability  or  instability  of 
critical  points  without  actually  solving  the  system.   This 
is  particularly  important  in  cases  where  the  solutions  are 
impossible  to  obtain  in  closed  form  or  are  difficult  to 
analyze.   We  can  still  gain  valuable  information  on  their 
stability  provided  that  v/e  can  construct  a  Liapunov  function 
for  the  system.   Although  there  is  no  general  method  for 
constructing  these  functions,  v/e  shall  see  that  this  is 
fairly  easy  to  do  in  many  cases.   To  solidify  our  under- 
standing of  this  notion,  we  consider  the  following  example. 

•  2     2 

EXAMPLE  The  function  E(x,y)  =  x  +  y   is  a  Liapunov 

function  for  the  system 


X  =  -X  +  X  y. 


2 
y  =  -y  +  xy  . 


Clearly,  E  is  positive  definite  and  of  class  C   in  the 
entire  x,y-plane.   Moreover,  since 

'E  =  H  ^1  ^  If  ^2  =  2xe^  +  2ye2  and 


lil8 


F(x,y)  =  (-X  +  x^y)e-^  +  (-y  +  xy^)e. 


we  have 


VE-F  =  -2(x^  +  y^)  +  2x^y  +  2xy^  =  2(x^  +  y^)(xy  -  1) 

Hence  -VE'F  is  positive  definite  in  the  region  xy  <  1 ,  and 
E  satisfies  all  the  requirements  of  Theorem  6.12  in  this 
region.   Thus  the  origin  is  an  asymptotically  stable 
critical  point  for  this  system. 

We  shall  now  describe  a  method  of  constructing  a 
Liapunov  function.   To  do  so  we  begin  by  constructing  a 
Liapunov  function  for  the  constant  coefficient  linear 
autonomous  system 


AX 


when  all  the  eigenvalues  of  A  have  negative  real  parts. 

It  is  noted  that  we  previously  determined  that  the  origin 

is  asymptotically  stable  for  this  system.   Thus  we  will 

not  gain  anything  new  from  the  construction  but  it  will 

allow  us  to  extend  some  of  our  earlier  results  to  non-linear 

systems.   The  construction  goes  as  follows. 

Let  e.,  ,  .  .  .  ,  e   denote  the  standard  basis  vectors  in 
1        n 

R^  and,  for  each  i,  1  <_  i  <_  n,  let 


ihg 


X^(t)  =  X^(e^,t)  = 


x^^(t) 


^n^^) 


denote  the  solution  of  this  system  which  satisfies  the 


initial  condition  X. (0)  =  e . .   Then  if 


n 


Y  =  y-|_e^  +  .  .  .  +  y^e^  is  any  vector  in  R  ,  the  funct 


ion 


X(Y,t)  =  y^X^(t)  +  ...  +  y^X^(t) 


is  the  solution  of  the  system  which  satisfies  the  initial 
condition  X(0)  =  Y.   We  now  set 


E(Y)  =   /  I |X(Y,t)|  p  dt, 
0 


and  note  that  if  this  integral  converges,  then.E(Y)  is 
positive  definite.   It  is  clear  that  E(Y)  >_  0.   Also  if 
Y  =  0,  then  for  all  i,  1  <_  i  <_  n,  we  have  y.  =  0  and  hence 
X(0,t)  =  0.   Thus  E(0)  =  0.   Now  suppose  E(0)  =  0.   This 
implies  that  X(Y,t)  =  0.   Hence  y.  =  0  for  all  i,  1  1  i  1  n 
since  otherwise  X. (t)  must  be  zero  for  some  i.   But  then 


X.(0)  i-   e..   Thus  we  must  have  that  Y  =  0. 


.n 


The  proof  that  E(Y)  is  defined  for  all  Y  in  R  can  be 
found  in  Ostberg,  Kreider,  Kuller,  p.  4l5. 

Finally  we  show  that  E  is  in  fact  a  Liapunov  function 
for  this  system.   We  compute  the  value  of  E(X(XQ,t))  where 
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X(XQ,t)  Is  the  solution  of  X  =  AX  satisfying  X(X  ,0)  =  X  . 


Hence  we  have 


E(X(XQ,t))    =      ;    I  |X[X    ,t),s]|  pds    =      /    ||X(X      s    +   t)||^   di 


=      /    I |X(X^,u)|  p    du, 
t  ^ 


Hence 


^9E(X(XQ,t))    ^   l_     ;"||x(Xq,u)||^    du   =    -||X(XQ,t) 


Zt 


dt 


as  required  since  VE-F  =  8E/9t  along  the  trajectories  of 

2 

X  =  AX.   Moreover,  since  the  real-valued  function  | | X| | 

is  positive  definite  on  R  ,  E  satisfies  all  the  hypotheses 
of  Theorem  6.12  and  again  we  have  proved  that  the  origin 
is  asymptotically  stable. 

We  now  use  this  method  to  construct  a  Liapunov  function 
for  the  following  linear  system: 

X  =  -X 

y  =  -2y. 
This  system  has  the  coefficient  matrix 
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with  eigenvalues  -1,  -2.   In  this  case  the  solution  X.(t) 
and  X^Ct)  such  that  X^(0)  =  e^,  X^CO)  =  e^  are  given  by 


r  -t 


x^(t)  = 


X2(t)  = 


-2t 


Thus  if  Y  =  xe   +  ye  , 


E(Y)  =   /  I  |xe"^e,  +  ye'^'^eJ  |^  dt 
0        -^        '^ 


,  ,  2  -2t  _^  2    -H.    ^^ 
=  /(xe    +ye    )dt 


2     2 

^   +  y 


The  fact  that  this  function  satisfies  the  hypotheses  of 
Theorem  6.12  is  easily  verified  and  thus  the  origin  is 
asymptotically  stable. 

We  now  wish  to  extend  these  results  to  non-linear 
systems.   For  this,  we  recall  the  notion  of  the  linear 
approximation  which  was  developed  in  the  last  section. 
We  recall  that  in  the  system 

X  =  F(X), 

the  function  P  may  be  replaced  by  a  function  of  the  form 
JX  +  G(X),  v:here  J  is  the  n^n  Jacobian  matrix  of  F  and 
G(X)  is  small  in  comparison  with  JX  when  ||x||  is  small. 
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Now  let  F  be  a  function  of  class  C   in  a  region  fi  of 
R  containing  the  origin,  and  suppose  that  P(0)  =  0  and 
the  real  parts  of  all  of  the  eigenvalues  of  the  Jacobian 
matrix  J  of  F  evaluated  at  the  origin  are  negative.   We 
now  show  that  under  these  conditions  the  Liapunov  function 
E  =  E(X)  for  the  constant  coefficient  linear  system 
X  =  JX  that  has  been  previously  constructed  is  also  a 
Liapunov  function  for  the  non-linear  system 

X  =  F(X). 


2 

Now  by  our  earlier  results,  we  know  that  VE« (JX)  <_  -||x|| 

Hence 


VE(X)-F(X)  =  VE(X)-[JX  +  G(X)] 

=  VE(X)-JX  +  VE(X) -GCX) 


<  -||X| 1^  +  I |VE(X)| I  I |G(X) 


This  last  statement  is  a  consequence  of  the  Cauchy-Schwarz 
Inequality.   Now  we  also  know  that  ||VE(X)||  <_2k||x||, 
since  VE  is  a  linear  transformation  on  a  finite-dimensional 
space.   Thus  we  have 

VE(X)-F(X)  <_   -I  |X|  1^  +  2k|  |X|  I  I  |G(X)|  |. 
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We  now  use  the  fact  from  the  definition  of  a  linear 

approximation  that     lim     i  i^if^^  =  0.   But  then  it 

lUlho         "XII 

G(X) 
follows   that        lim        i i^i  r      =   0.      Hence   we    can   find   a 

Ijxlho  ll^ll 

sphere  S[0,6]  about  the  origin  in  r"  such  that 


l|G(X)||  <  4  IIXJI 


for  all  X  e   S[0,6].   Then  within  this  sphere  we  have 

VE(X).F(X)  1  -||X||^  +  (2k||x||)(^||x||)  =  4M^II^ 

and  indeed  E  is  a  Llapunov  function  for  the  system 
X  =  F(X).   Moreover,  since  this  inequality  also  implies 
that  -VE'F  is  positive  definite  in  S[0,5],  the  hypotheses 
of  Theorem  6.12  are  satisfied  and  we  have  established  the 
following  result. 

THEOREM  6.14   Let_  F  =  F(X)  be_  of  class  C"^  in_  a  region  of 
R   containing  the  origin ,    and  suppose  that  F(0)  =  0.   Then 
the  origin  is  asymptotically  stable  for  the  system  X  =  F(X) 
whenever  all  of  the  eigenvalues  of  the  Jacobian  matrix  of 
F  evaluated  at  the  origin  have  negative  real  parts . 

As  before,  there  is  a  companion  result  when  all  of  the 
eigenvalues  of  the  Jacobian  matrix  of  F  have  positive  real 
parts.   In  that  case  the  origin  is  a  point  of  complete 
instability  for  X  =  F(X). 
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These  statements  verify  the  comments  in  Section  V 
about  the  conclusions  that  may  be  dravm  from  a  linear 
approximation.   Finally  we  show  by  means  of  an  example 
that  no  definite  conclusions  can  be  made  from  a  system 
whose  linear  approximation  shows  that  the  origin  is  a 
stable  critical  point. 

EXAMPLE  Consider  the  plane  autonomous  system 

X  =  y  -  xf(x,y) 
y  =  -X  -  yf(x,y) 

where  f  has  a  convergent  power-series  expansion  in  a 

2 

neighborhood  of  the  origin  in  R  and  f(0,0)  =  0. 

Taking  the  linear  approximation,  we  have 
X  =  JX 
where 


J  = 


0 

-1 


In  this  case  J  has  eigenvalues  ±i  so  by  Theorem  5.8,  the 
origin  is  stable  for  the  system  X  =  JX. 
Now  let  E(x,y)  =  x^  +  y  .   Then 
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^^  ^    ^2^"  \(-x  -  yf(x,y))j     2^^  ^  y  )f(x,y). 


Hence  for  this  system  the  origin  is 

1.  stable  when  f(x,y)  _>  0  in  some  neighborhood 
of  (°); 

2.  asymptotically  stable  when  f(x,y)  is  positive 
definite  in  a  neighborhood  of  (q); 

3.  unstable  when  f(x,y)  <  0  in  every  neighborhood 
of  (°). 

This  completes  our  study  of  Liapunov  functions  and 
stability  theory  in  general.   At  this  point,  the  reader 
should  have  sufficient  background  in  this  theory  to  be 
able  to  read  and  understand  further  results  in  this  area. 
Specifically  recommended  for  further  reading  is  Chapter 
VIII,  C   -Liapunov  Functions  for  Ordinary  Differential 
Equations,  in  Bhatia  and  Szego. 
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